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^ ' Abstract: Splitting amplitudes are universal functions governing the collinear behavior 

of scattering amplitudes for massless particles. We compute the two-loop g ^ gg split- 
ting amplitudes in QCD, M = 1, and = 4 super- Yang-Mills theories, which describe the 
limits of two-loop n-point amplitudes where two gluon momenta become parallel. They 
5^ I also represent an ingredient in a direct x-space computation of DGLAP evolution kernels 

at next-to-next-to-leading order. To obtain the splitting amplitudes, we use the unitar- 
ity sewing method. In contrast to the usual light-cone gauge treatment, our calculation 
does not rely on the principal- value or Mandelstam-Leibbrandt prescriptions, even though 
the loop integrals contain some of the denominators typically encountered in light-cone 
gauge. We reduce the integrals to a set of 13 master integrals using integration-by-parts 
and Lorentz invariance identities. The master integrals are computed with the aid of 
differential equations in the splitting momentum fraction z. The e-poles of the splitting 
amplitudes are consistent with a formula due to Catani for the infrared singularities of 
two-loop scattering amplitudes. This consistency essentially provides an inductive proof 
of Catani's formula, as well as an ansatz for previously-unknown 1/e pole terms having 
non-trivial color structure. Finite terms in the splitting amplitudes determine the collinear 
behavior of finite remainders in this formula. 



Keywords: 



QCD, NNLO Computations, Jets, Collider Physics. 



1. Introduction 



Gauge theories form the backbone of the standard SU(3) x SU(2) x U(l) model of particle 
interactions. The computation of perturbative corrections in gauge theories is thus central 
to testing the standard model at high-energy colliders. Such computations are technically 
complicated, so a general understanding of properties of the results is very useful. 

In the past decades, a number of new approaches have been developed to cope with 
this complexity, including helicity methods [0, color decompositions ||2|, |^, ^, |6|, recursion 
relations 0, ideas based on string theory [^, and the unitarity-based method [10, 11 



12, O. The latter technique has been applied to numerous calculations, most recently the 



two- loop calculation of all helicity amplitudes for gluon-gluon scattering [|14| , 15]. 

The subject of two- loop calculations has seen tremendous technical progress in the last 
five years. Much of the progress has been facilitated by new techniques for performing loop 



integrals. Smirnov [|l^] and Tausk [17| gave closed- form expressions for the all-massless 
planar and non-planar double box integrals, respectively. Smirnov and Veretin |l8| and 
Anastasiou et al. |[l9[| provided algorithms for reducing the corresponding tensor integrals. 



More general reduction and evaluation techniques for integrals have followed as well po|. 



21, 22, uM. Using these techniques, several groups have computed the basic two-loop 



QCD (and also QED) amplitudes for four external partons p|, 26, 27, 28, 25, 3C, 31, 32 



15, pSl 34, l3H|, and for three partons and one external vector boson |36, 37, 22|. These 



amplitudes and matrix elements constitute one of the building blocks for next-to-next-to- 
leading order (NNLO) computations in perturbative QCD, in particular of the cornerstone 
processes e^e~ — > 3 jets and — > 2 jets. 

The two- loop amplitudes contain both ultraviolet and infrared singularities. Both 
are conventionally regulated using dimensional regularization. The infrared singularities 
correspond to the circulation of gluons in the loops that are nearly on shell and either soft 
or collinear with one of the external momenta (or to the circulation of other particles in 
the loops that are nearly on shell and collinear with one of the external momenta). The 
form and functional dependence of these singularities was predicted by Catani |^]. The 
a priori knowledge of these infrared singularities has been of great value in the explicit 
computations of two-loop amplitudes cited above. 

The singularities in the virtual corrections have a counterpart in the infrared singulari- 
ties of phase-space integrals of real-emission amplitudes. In the sum over all relevant (phys- 
ically indistinguishable) amplitudes, unitarity in the form of the Kinoshita-Lee-Nauenberg 
theorem dictates that the singularities must cancel. The universality of the singularities 
in the integrals over real-emission amplitudes in turn reflects the universality of factor- 
ization of real-emission amplitudes in soft or collinear limits. For next-to-next-to-leading 
order (NNLO) calculations, the double-emission limits of tree amplitudes, and the single- 
emission limits of one- loop amplitudes, allow one to organize the singular phase-space 
integrations in a process-independent way. The universal functions governing these limits 
are aU known ||, |o|, |4|, ||, ||, ||. 
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In this paper, we will compute the analogous functions, splitting amplitudes, governing 
the universal behavior of two-loop amplitudes as two gluon momenta become collinear. 
For example, in the time-like case where momenta ka and kb are both outgoing, we let 
ka zkp, kf) {1 — z)kp, where kp = ka + kf, is the momentum of the nearly on-shell 
intermediate gluon P, and z is the longitudinal (or light-cone) momentum fraction carried 
by gluon a. Splitting amplitudes will be useful in checking two-loop computations beyond 
four external legs. A putative result for a higher-point amplitude must satisfy non-trivial 
constraints as momenta become collinear. They can also play a role in an alternative 
method of computing the NNLO corrections to the Altarelli-Parisi kernel governing the 
evolution of parton distributions and fragmentation functions |45]. The computation of 
this kernel is of great importance to a program of precision extraction of parton distribution 
functions from experimental data. It has been the object of an ongoing effort by Moch, 
Vermaseren and Vogt p6| , p7[| , just recently completed [48|. The Mellin moments of the 
Altarelli-Parisi kernel are anomalous dimensions of leading-twist operators whose matrix 
elements give rise to parton distributions. One can compute them in a traditional manner, 
by computing ultraviolet divergences of loop corrections. Factorization implies, however, 
that one could compute the kernels in an infrared approach, directly in x-space [49, 50 1. As 
in the computation of differential cross sections, there are both 'virtual' and 'real-emission' 
contributions to the kernel. The splitting amplitudes we compute in the present paper 
provide the doubly-virtual contributions to the NNLO kernel for evolution of the gluon 
distribution, Pgg {x) (at x / 1). Because the gluons have definite helicity, evolution of 
polarized distributions is equally accessible. 

The computation of the Altarelli-Parisi kernel in A/" = 4 supersymmetric gauge theories 
is also of interest, since its Mellin moments are the anomalous dimensions of classes of 
operators. The study of such anomalous dimensions is important to investigations of the 
anti-de Sitter /conformal field theory duality The NNLO anomalous dimensions were 
very recently extracted ||5^ from the corresponding QCD calculation |4^], by adjusting 
the fermion content to match = 4 super- Yang-Mills theory, and then keeping only those 
terms that are leading in "transcendentality." (In x-space, a logarithm or a factor of 
vr is assigned one unit of transcendentality, and an n**^ degree polylogarithm Li„(x) is 
assigned n units. The leading- in-transcendentality terms at two loops carry 4 units of 
transcendentality.) Our results confirm that the leading-transcendentality procedure of 
ref. does correctly extract the AA = 4 result from that of QCD, for the fully virtual 
NNLO contributions to the Altarelli-Parisi kernel, which arise from the two-loop splitting 
amplitude (and from the square of the one-loop splitting amplitude) . 

We gave the result for the splitting amplitude in AA = 4 supersymmetric gauge theory 
in an earlier Letter We will document that calculation in the present paper. The 

computation revealed an unexpected relation between splitting amplitudes at different 
loop orders: the two-loop splitting amplitude can be expressed algebraically in terms of 
the one- loop and tree splitting amplitudes, through 0{e^), where e is the parameter of 
dimensional regularization, D = A — 2e. In the planar (large-A^c) limit, the four-point 
two-loop amplitude can be expressed in a similar 'iterative' form. This is very surprising 
because in a general massless field theory, the analytic structure of a two-loop amplitude 
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can be considerably more complicated than that of one-loop amplitudes. Thus the two- loop 
amplitude in this theory is much simpler than expected. The splitting amplitude relation 
also led to a conjecture of a similar relation between one- and two-loop amplitudes with 
an arbitrary number of external legs. 

To calculate the g ^ gg splitting amplitudes, we have used the unitarity-based method. 
The method is useful in general loop calculations in gauge theories. Its advantage over 
Feynman-diagram calculations is of course most obvious in those calculations which simply 
cannot be done by conventional techniques, such as that of infinite series of loop ampli- 



tudes [|10|, 11 , 54 1 . The present calculation furnishes another example where the method has 



a clear advantage over conventional diagrammatic techniques. While the splitting ampli- 
tudes could be computed by conventional diagrammatic techniques, such a computation at 
two loops would probably require the use of light-cone gauge, because this gauge has simple 
collinear factorization properties, even in the presence of infrared singularities |5^. (An 
analysis in a covariant gauge is likely to be very difficult; for example, a generalization of 
the analysis in ref. [56| would require explicit knowledge of higher-point two- loop integrals.) 
As is well-documented in the literature [^, use of light-cone gauge is fraught with sub- 
tleties and technical complications. Indeed, following the standard methods for dealing with 
light-cone gauge Feynman diagrams, along with any of the popular prescriptions needed to 
avoid ill-defined integrals, would lead to an answer with a surviving dependence on the pre- 
scription parameters, which cannot describe the collinear behavior of a gauge- invariant and 
prescription-parameter-independent amplitude. We will discuss these issues in more detail 
in section ^. The unitarity-based method avoids these complications, and ensures that the 
calculation can be done in a straightforward way using ordinary dimensional regularization. 
The insights furnished by the unitarity method also allow the systematic cancellation of 
ill-defined integrals appearing in more traditional light-cone gauge calculations [59|. 

To evaluate the loop momentum integrals we used integration-by-parts and Lorentz 
invariance [20| identities implemented via the Laporta algorithm |23, Sl| to solve the system 
of equations. With this technique the integrals can be reduced to a set of 13 master 
integrals. The master integrals are computed by constructing a set of differential equations 
in z, along the lines of refs. [18, |20| ]. 

We will organize the calculation in a color-stripped formalism, in which the color 
factors are separated from the kinematic content of amplitudes. At loop level, this leads to 
a hierarchy of terms, from terms leading in the number of colors, down through subleading 
contributions. The leading-color terms, which would dominate in the Nc ^ oo limit, 
correspond to planar diagrams. These leading-color terms contain only a single color trace, 
with additional explicit powers of Nc- The two-loop splitting amplitude enters only into 
the collinear limit of single-trace terms (whether leading or subleading in color). The 
collinear limits of terms with multiple traces depend only on the tree-level and one-loop 
collinear splitting amplitudes. The full color-dressed splitting amplitude also factorizes into 
a color factor multiplied by a factor which is a function solely of the external momenta and 
helicities. The color factor for the g ^ gg splitting amplitude, in the language of Feynman 
diagrams, is always just a Lie algebra structure constant, f"''"^. At two loops, the pure-glue 
result is purely leading-color; there are no subleading-color corrections. Adding quarks in 
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the loops does generate subleading-color terms. However, non-planar diagrams do not give 
any contribution to the g ^ gg collinear behavior at two loops, because their color factors 
all vanish. (This feature will not hold for g ^ gg at three loops, and already fails to hold 
at two loops for splittings with external quarks, q ^ qg and g — > qq.) 

As noted above, Catani gave a formula which predicted the infrared singularities 
of renormalized two-loop amplitudes, which appear as poles through fourth order in e. 
One can take the collinear limit of the n-point formula, and compare it with the (n — 1)- 
point formula, to obtain a prediction for the infrared singularities of the two-loop splitting 
amplitude. The Catani formula is expressed in terms of color-space operators. One term 

" (2) 

arising at order 1/e (denoted by ) contains a rather non-trivial color structure, and 
was known previously only for the four-point case. We have constructed a simple ansatz 
for its generalization to n-point amplitudes. A comparison of the Catani formula with our 
result requires an untangling of color correlations. We have performed this analysis, and 
find complete agreement. One can interpret this agreement as a proof of Catani's formula 

" (2) 

(including our ansatz for ) for the case of two-loop n-gluon amplitudes. The proof is 
inductive in the number of legs n, and requires certain reasonable assumptions about the 
types of functions that can appear in singular terms. (The functions should not vanish 
in all collinear limits. Such vanishing is unlikely to happen for functions which are equal 
to the tree amplitude times pure logarithms or polylogarithms, for example.) This proof 
complements the verification based on resummation given in ref. |62]. 

Given the consistency of our results with Catani's formula, we can subtract the singular 
behavior in e, to obtain a set of relations which control the finite remainder terms in the 
formula. We have checked that these relations are satisfied, up to overall normalization, 
by the finite remainder for the two- loop H — > ggg amplitude |63] for the case of identical- 
helicity gluons. 

In the next section, we review the structure of gauge-theory amplitudes in their 
collinear limits. In section |^, we consider a hypothetical Feynman-diagram calculation, 
both in a covariant gauge and in light-cone gauge. We discuss the difficulties that would 
be encountered in these calculations, and how a unitarity-based method can circumvent 
them. In section ^ we review the unitarity-based sewing method, and present a detailed 
algorithm. (Appendix ^ describes a simple relabeling algorithm, used in the sewing pro- 
cess to bring integrands into a canonical form.) In section |5|, we discuss the application of 
the sewing algorithm to the calculation of the two-loop splitting amplitude. In section |6|, 
we describe the calculation of the required two-loop integrals, and the reduction of tensor 
integrals. We present our results in section 0. These contain the complete set of helicity- 
decomposed g ^ gg splitting amplitudes in QCD with fermions, as well as in = 1 
and AA = 4 supersymmetric gauge theories. We present separately the collinear behavior 
of the finite terms in a two-loop amplitude. In section we present the comparison of 
the singular parts of the color-trivial terms to those predicted by Catani's formula, as well 
as the collinear behavior of the finite remainder terms. (The precise definition of 'color- 
trivial' is given in that section. The full color dependence of the singular terms, including 

" (2) 

the ansatz for if„ , and their collinear behavior, are discussed in appendix ^.) We discuss 
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the dressing of the splitting amplitudes with color factors, and the collinear behavior of 
the color-non-trivial parts of amplitudes, in section ^ We give some concluding remarks 
in the final section. 



2. Splitting Amplitudes 



Gauge-theory amplitudes are singular when external momenta become soft, or when a 
number of momenta become collinear. In these limits, the amplitudes factorize in a univer- 
sal way. In collinear limits, the factorization is governed by splitting amplitudes depending 
only on the legs becoming collinear, and not on the remaining legs of the hard process. The 
surviving hard amplitude, in turn, depends only on the merged leg, whose momentum is 
the sum of the collinear momenta. In this paper, we will consider the splitting amplitude 
governing the behavior of amplitudes as two momenta become collinear. (The tree-level 
behavior as three or four momenta become collinear has been derived by Campbell and 
Glover |4C], Catani and Grazzini |41|, and by Del Duca et al. |6^. The one-loop behavior 
as three momenta become collinear has recently been considered by Catani, De Florian 
and Rodrigo p5[|.) 

We find it most natural to discuss collinear factorization in the context of a trace-based 
color decomposition of the n-gluon amplitudes ||2|, ^, ^, |5| . Although the external gluons are 
in the adjoint representation, this color decomposition is given in terms of traces of matrices 
T"" in the fundamental representation of SU(A'c), which we normalize by Tr(T"T'') = S""^. 
We begin by discussing the behavior of terms leading in the number of colors, Nc, where 
we scale the number of fermions, Nf, with Nc- The full color behavior is a straightforward 
extension, which we defer to Section |^. 

The loop expansion of the n-gluon amplitude is 



L=0 



(47r 



|2-e 



j(L)ai...a„ 



, Ai; . . . ; A;„, An) , (2.1) 



where Oj is the color index of the i-th external gluon, and the factor of [2e~'^'^/(47r)^^'^]-^, 
with 7 = — V'(l) = 0.5772 . . ., corresponds to the normalization convention of ref. |3^]. The 
trace-based color decomposition of the L-loop amplitude is 



Tr(a(l)...a(n))4^)(a(l),...,cT(n)) + 



2.2) 



where A^j"^ are L-loop color-ordered (sub-)amplitudes; where Tr(l . . . n) = Tr(T"^ . . . T""); 
and where a runs over the non-cyclic permutations Sn/Zn of {1, 2, . . . , n}. The latter 
correspond to the set of inequivalent traces. The permutation a acts both on the gluon 
momenta ki and helicity labels Aj, implicit on the right-hand side of eq. (|9.l|). At tree 
level (L = 0), expression (2^) is exact, and has no subleading-color corrections. At loop- 
level {L > 0), there are subleading-color terms containing products of two or more traces 
(see eqs. ( p.4[ ) and (^^)). Also, the amplitudes must be evaluated in D = 4 — 2e dimensions 
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Figure 1: The collinear factorization of a tree-level amplitude. The thick line represents a slightly 
off-shell gluon. 



to regulate the virtual singularities. We extract in a prefactor the leading behavior. 
The normalization of An, and hence of the splitting amplitudes we present below, differs 
from that in refs. |T^. Note that ^4^^''''^ will in general contain terms proportional to 
powers of Nf/N^ and l/N'^. 

The color-ordered amplitudes in eq. (p^) have a universal behavior as legs become 



collinear [S6|. At tree level, the amplitudes behave as 

. . , a^", . . .) * E SplitL°l(^; &^^) • • , • • •) , (2-3) 

A=± 



in the limit where the momenta ka zkp and ^ {1 — z)kp with kp = ka + kb- Here 
is a tree-level splitting amplitude. Legs a and b carry helicities and 
Xb, while the merged leg P carries helicity A. In the sum A runs over the two helicities 
of the intermediate state. The factorization of an n-point tree amplitude into a splitting 
amplitude and an (n — l)-point amplitude (2.3) is depicted schematically in figure |l|. 



The pure-glue tree-level splitting amplitudes are 

SplitL°^(z;a-,6-) = 0, (2.4) 

Split^°^(z;a+,6+) = ^ L = , (2.5) 

v/^(r^(a6)' 

SplitL°^(z;a+,r) = - ^— (2.6) 

^yz{l - z) [ah] 

SplitL°)(.;a-,6+) = - . . (2.7) 

Vz(l - z) [ah] 

These splitting amplitudes are expressed in terms of spinor inner products [|^, (ij) = 
and [ij] = where [i^) are massless Weyl spinors of momentum ki, labeled 

by the sign of the helicity. The spinor products are antisymmetric, with norm | {ij) | = 
I [^j] I — \f^j^ where Sij = 2ki ■ kj. A key advantage of the spinor formalism is that it 
makes the square-root behavior of the splitting amplitudes manifest. The remaining tree- 
level splitting amplitudes, with A = — , may be obtained from the above ones by parity, 
which states (for general loop order L) that 

Split^^} , J z;a-^\b~^'^) = -S-plit^^hz;a^\ b^^) . (2.8) 
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+ 



b 



a 



a+b 



Figure 2: Two types of terms appear in the collinear factorization of a one-loop amplitude. The 
thick lines represent slightly off-shell gluons. 



In quoting explicit results in this paper, we use parity to assume that the intermediate 
state P always has positive helicity, A = -|-. As we shall discuss further in section |9|, 
Bose symmetry implies that the color-stripped splitting amplitude is antisymmetric under 
exchange of its two arguments (including z ^ 1 — z), 



Split - z; h^\a^-) = - Split a^^h^'^) 



(2. 



This relation allows us to obtain the results for P'^ a^b from those for P~^ 



-b+. 



At one loop the structure is similar. In this case the the collinear limits are [10, 56, 66 1 



a\\b 



( SplitL°l(z; , 6"") ■■,P\---) + 

+ SplitL'l(z;a^",6^'')41(...,P\...)) , 



(2.10) 



where Split is a one- loop splitting amplitude. Figure |2| displays eq. ( 2.10| ) schematically. 

Except for Split_(z; a~, b~), which vanishes at tree level (along with Split^(z; a"^, 6^)), 
the ratio of the g ^ gg one-loop splitting amplitudes to the tree-level ones is well-defined. 
The ratio depends trivially on the Lorentz invariant Safe, and non-trivially on z, but it does 
not involve spinor products. For the helicity configuration whose tree splitting amplitude 
vanishes we have, 



SplitL^^(z;a~,6-) = +CT^/4l 



{ab) 



abf (1 

2 



-•Safe 



-26)(2- 

l-e6R 



26) (3 



2e) 



where 



Cr 



e'^ r(l + e)r2(l-e) 



(2.11) 



(2.12) 



2 r(l-2e) 

Notice that in the e-expansion of (fJ-'^ /{—Sab)Y , lii(— s^fe) appears. For time-like kinematics, 
this expression has an imaginary part according to the prescription, 



ln(— s 



ab) 



In Sab - 



Sab > 0. 



(2.13) 



-7- 



The remaining splitting amplitudes, for A = +, are conveniently written in terms of 
their ratios to the tree-level ones, 



SplitL^^(z;a^",6^'') = r^^'^^''^'{z,Sab) x Split (z; a^-, ^^i"). 



(2.14) 



where 



,{!)++, QCD 



^(1)-+,QCD 
S 



and 



Z,S 



{z,s) 



{z,s) 



cr 



cr 



r 



s 



,(1),A^=4 



2 



2z{l - z) 



(1 - 2e)(2 - 2e)(3 - 2e) 



'5 



(2.15) 
(2.16) 



vre 



+ ^2e 



2m- 1 



m=l 



z(l-z)(-s) 



sin(7re 
+ 21nz ln(l - z)-C2 



Li2m-1 



+ 0(e) 



(2.17) 



(l),A^=4 
5 



means that it is 



(Here, is the Riemann zeta function.) The label = 4 on r 
the appropriate function for = 4 super- Yang-Mills theory, for both -|— |- and — h cases 
(see also section [0| ). At one loop this expression also happens to serve as the splitting 
amplitude for pure M = \ super- Yang-Mills theory, for both non-vanishing cases. 



(i)++,.^=i 



z,s) 



^(i)-+,^=i 
s ' 



z,s) 



„(l).A/'=4, 



z,s) 



(2.18) 



This relation will be violated at two loops. The parameter 5^ selects the particular variant 
of dimensional regularization (see eq. ([7.1|)). For = 1 the scheme is the 't Hooft-Veltman 
(HV) scheme, while for = it is the four-dimensional helicity |7^ (FDH) scheme. 
We always quote results for supersymmetric theories in the FDH scheme, which is related 



to, but distinct from, Siegel's dimensional reduction scheme (As noted above, the 

normalization of the splitting amplitudes differs from that in ref . |jl^ ; to recover the earlier 
normalization for Split^^\ replace cr by cr, defined in eq. (2.12) of that reference. To 
recover the earlier normalization of r^^\ replace cr by unity.) 

These one-loop splitting amplitudes were first obtained from the coUinear limits of 
five-point amplitudes [^]. Subsequently they were obtained to all orders in e, as required 
for NNLO calculations |5^, 42, 43 1. We will express divergent parts of the two- loop splitting 
amplitudes in terms of one- loop quantities. Hence we have retained all the terms in the e 
expansion in these expressions. 

At two loops, the subject of this paper, the amplitudes behave as [l66|, 



all?) 



^ rSplitL°i(z; a^" , 6^") 4'1,(. . . , P\ . . 

\ L V 

+ SplitL'i(z;a"%6''')4-i(---,^'' 



A=± 



+ 



+ SplitL'i(^;a^",6^^)4°2,(...,P\...)) , (2.19) 
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(2) 

where Split^j^ is the two-loop sphtting amphtude. One of the goals of this paper is to 
calculate this two-loop splitting amplitude in QCD. To do so we will use the unitarity 
sewing method [|lO|, 11, |l^], as applied to splitting amplitudes [Q. In a previous paper 
we presented the result of this calculation for the special case of = 4 super- Yang- Mills 
theory. For this theory the unique two-loop splitting amplitude has the remarkable property 
of being an iteration of the one-loop result ( ^.17 ), which led to a conjecture that a similar 



iterative property holds for the planar contributions to amplitudes |53|. This conjecture was 



shown to be correct for the four-point amplitude using a previously-derived [72| expression 
for the two-loop integrand. In section ^ we shall present the explicit values of the two-loop 
splitting amplitudes for g ^ gg in QCD, as well as in = 4 and J\f = 1 supersymmetric 
gauge theories. 

As in the one-loop case, we write the two-loop splitting amplitudes in terms of their ra- 
tios rs to the corresponding tree- level splitting amplitudes (when the latter do not vanish). 

(2) 

Taking A = + by parity, we define via. 



Split^^\z;a^-,b^^) = rP^^^'{z,Sab) x SplitL°^(z; 6^'') 



(2.20) 



We give the explicit expansions of the functions r^p'^'''^'' 



as Laurent expansions in e, for 
the various theories in question, in eqs. (|7T^ ), (|7!28|) , (|7!29|) and (^^), and the expression 
for SplitL^^(z;a-,6-) in eq. (ffl^) . 

Although we shall not discuss higher-loop splitting amplitudes here, we remark that 
the obvious L-loop generalization of the collinear behavior. 



a\\b 



E E Split 

1=0 \=± 



u(0 



{z;a 



(2.21) 



can indeed be proven correct ||66[. These splitting amplitudes govern the collinear behavior 
of the entire amplitude, including all multiple color trace terms. We will discuss this more 
fully in section ^. With our normalizations the leading-color contributions to Split are 
of order N^, but in general the splitting amplitudes have contributions of higher order in 
as well as quark-loop contributions of order (Nj/NcY, with p < L, oi higher in 

1/Nl 



3. Difficulties with Feynman Diagram Approach 

Before turning to our calculation of the two-loop g ^ gg splitting amplitudes using 
the unitarity-based sewing method, it is instructive to consider how one would proceed 
using a standard Feynman diagram approach. Heuristically, one might try to 'factorize' an 
n-point amplitude on the collinear pole, i.e. to construct the L-loop splitting amplitudes 
by summing up all L-loop Feynman diagrams with three external legs, one of which is 
off-shell. While this might seem to be the most straightforward approach, a number of 
complications arise in practice, as we shall see in this section. In section |5| we shall sidestep 
these complications using the unitarity-based sewing procedure [10, 11, [l^, |l^], outlined 
in section This method has previously been applied to splitting amplitudes at one 
loop m, ||. 
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3.1 Tree-Level Splitting Amplitudes 

At tree level, factorization works without any subtleties. 
^ That is, we can compute the splitting amplitude directly from 

the Feynman diagram three-point vertex depicted in figure |3|. 
The only diagrams for an n-point amplitude with a pole in Sab 
^ are those containing this vertex as one factor. Such calcula- 

tions have appeared elsewhere, for example in ref. [|67| ]. 
Figure 3. The three-point As two color-adjacent momenta and A;^ become collinear, 

vertex diagram for obtaining p , . ■ , , i-, i ,i i ■ 

,. . we factonze an n-pomt tree amphtude on the Snh Icmematic 

the tree-level splitting am- 

J rr^i .1 • 1 1- pole in terms of a three-point vertex and an in — l)-point am- 

phtude. ihe thick line rep- ^ ^ \ j f 

resents an off-shell gluon. plitude, 

4°)(l,2,...,a,6,...,n) ^ 6,(a)..(b)y^-^ E ^^^^'^^^^ ) >< 



where the kinematics is the same as in eq. ( |2.3D . In setting up the calculation the merged 
leg should be left slightly off-shell. Otherwise the putative splitting amplitude would be 
ill-defined, since the 1/sab pole diverges. This basic structure is independent of the particle 
type, although we have written eq. ( p| ) for the case of an intermediate gluon. Here £^{P) 
is the polarization vector for the gluon P with helicity A. For the case of g ^ gg, the 
ordinary Feynman gauge three-point vertex is 

V'^'^Pika, h) = -J= (r/^^ {ka - hY + rf'P {2kb + kar - r]^' i'^ka + k,r) . (3.2) 



After inserting an explicit representation of the helicity states [y, we obtain from eq. 
precisely the collinear behavior (^]^), together with the explicit values of the splitting 
amplitudes ( |2.4| )-( p77D . In this limit the helicity algebra in the numerator of the vertex 
causes it to vanish like Sab, partially canceling the pole in Sab- Physically, this cancellation 
is due to an angular momentum mismatch between A and A" -|- A^. Overall, we are left 
with the 1/y/sab behavior evident in the splitting amplitudes. Other gauge choices, such 
as light-cone gauge or the non- linear Gervais-Neveu gauge ||7^, give the same final result 
as Feynman gauge. 

3.2 Difficulties at Loop Level 

The simplicity of a conventional diagrammatic calculation at tree level might lead one to 
believe that a similar approach should work at loop level. Such a calculation would involve 
Feynman diagrams of the type depicted in figure ^ This expectation, however, turns out 
to be incorrect. At loop level, the situation is not quite this simple. 

In covariant gauges, such as Feynman gauge, one immediately runs into trouble |^6|| , 
because there are contributions to collinear behavior from Feynman diagrams of the form 
shown in figure 0, with no single-particle pole in Sab- The required pole emerges only after 
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a+b 




a+b 




a+b 



Figure 4: Three non-vanishing Feynman diagrams contributing to the pure-glue splitting ampli- 
tude. The thick gluon line is slightly off shell. 



carrying out the loop integration. Moreover, the e-expansion of some of the integrals has 
discontinuous behavior as Sab — > 0, because of the interchange of this limit with the limit 
e — 0. The appearance of more complicated loop integrals is reflected in the structure of 
the known results for the splitting amplitudes ( 2.15| )-( [2.17] ), which contain polylogarithms 
in the variable z. Such functions simply cannot occur in the integrals encountered in the 
triangle or bubble graphs in figure which in a covariant gauge produce only logarithms, 
in the variable Sab- Clearly, the polylogarithms must come from elsewhere. They do in 
fact arise from non- factorizing diagrams of the type shown in figure || [56|. This peculiar 
absence of factorization is tied to the presence of soft and collinear virtual divergences in the 
theory, reflected as poles in e in the loop amplitudes. At one loop, the non-factorizing pieces 
may be reconstructed using knowledge of all the possible integral functions that appear in 
amplitudes, along with the universal structure of the infrared divergences. Indeed, such 



a reconstruction was used to prove universal factorization of the one- loop amplitudes |56| 
and to compute their explicit values This reconstruction, however, does not generalize 
straightforwardly to higher loops. 

As the next logical choice of a method for calculating 
a splitting amplitude directly from Feynman diagrams, one 




Figure 

gauges, 



might turn to light cone-gauge. This gauge is known |74, 55 1 
to have simple factorization properties. A key advantage of 
light-cone gauge is that only physical states propagate. This 
can help clarify various formal properties. Moreover, light- 
cone gauge vertex integrals do contain polylogarithms of the 
type appearing in the one-loop splitting amplitude. Light- 
cone gauge has a long history and has been used to address 
a wide variety of problems. For example, the first proofs of 
factorization in QCD between the hard and soft parts of a 
process were performed in this gauge ||5^. Another important 
example is the next-to-leading order (NLO) calculation of the Altarelli-Parisi evolution 
kernel in j;-space [74, |7^, Light-cone gauge has also been used for more formal 
purposes, such as the proof of finiteness of maximally supersymmetric (M = 4) Yang-Mills 
gauge theories 



5: In covariant 
diagrams without 
single-particle factoriza- 
tion can contribute to the 
splitting amplitude. 



In describing the collinear limit of amplitudes, light-cone gauge is useful because, as 
a physical gauge, it can prevent the non-factorizing graphs of figure |5| from contribut- 
ing. In a general covariant gauge, such graphs would inevitably mix under residual gauge 
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transformations with the graphs in figure ^. 

In hght-cone gauge, the one-loop sphtting amphtude is given by the sum of the three 
Feynman diagrams shown in figure ^ In this gauge there are no ghost contributions. 
Furthermore, all cactus diagrams, as well as bubbles attached to the massless external 
lines a and b, vanish in dimensional regularization, leaving only the three diagrams shown. 
The light-cone gauge Feynman vertices are the same as those in Feynman gauge. The 
propagator is now 

p'^ + ie 

where the light-cone projector is. 



^M^ = -J%7 (3-3) 



df,u = Vfiu — • (3.4) 

p ■ n 

Here p is the particle momentum, tj is the Minkowski metric, and n is a null vector (n? = 0) 
defining the light-cone direction. 

In carrying out a light-cone gauge calculation, one quickly runs into well-known techni- 
cal difficulties [^] arising from regions of loop integration where the light-cone denominator 
p-n = p"*" vanishes. The light-cone denominators introduce a new set of singularities in the 
Feynman integrals, some of which are not regulated by dimensional regularization. In order 
for a generic light-cone gauge diagram to be well-defined, a prescription for dealing with 
these singularities is needed. For example, the principal-value (PV) prescription replaces 

^ limV ^—^ + ^—.], (3.5) 



p ■ n <5-»o 2 \p ■ n + iS p ■ n — id 

where (5 is a regulator parameter. Another choice, better founded in field theory, is the 
Mandelstam-Leibbrandt prescription 57]. A comparison of these two prescriptions may 



be found in ref. ||78|]. 

The introduction of an additional prescription in a splitting amplitude calculation is 
problematic for a number of reasons: 

• The additional prescription complicates the calculation, and requires the computation 
of more difficult Feynman integrals. 

• At higher loops the validity of the prescription is less clear ||5 

• After expanding in small 5, the results contain factors of In 5, where 6 is the pre- 
scription parameter. In general these cancel only after combining virtual and real 
emission contributions. 

Further to the last point, a calculation of the splitting amplitudes that retains depen- 
dence on a prescription parameter cannot match the splitting amplitudes extracted from 
the collinear limits of scattering amplitudes. An on-shell n-point scattering amplitude is 
gauge invariant, and depends only on the external momenta and on the dimensional regu- 
lator parameter e, not on any light-cone prescription parameter 5. The same is true of its 
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(a) (b) 

Figure 6: One-loop triangle diagrams containing light-cone denominators. The arrow indicates 
the line containing a light-cone denominator. Integral (a) contains an unregulated singularity, but 
does not appear in the unitarity-based sewing method. 



collinear limits. Thus loop splitting amplitudes defined via the collinear limits of scattering 
amplitudes cannot depend on the parameter S. 

If we had been computing the Altarelli-Parisi kernel, which receives contributions 
from both virtual and real emission contributions, then the dependence on 6 would cancel 



between them |75|, |76|, |77|, |79[. However, since we are interested in computing the collinear 
behavior of the virtual contributions on their own, in general the 5 dependence would not 
cancel. This difficulty makes it unclear how the desired splitting amplitudes can emerge 
from a light-cone gauge calculation. 

It is helpful to consider a few examples in order to illustrate the structure of inte- 
grals with light-cone denominators, and to help explain how we will later sidestep these 
difficulties, using the unitarity-based sewing method. 

3.3 One-Loop Light-Cone Integral Examples 

Consider the scalar integrals shown in figure ^. In the figure, an arrow on an internal line 
indicates the insertion of the light-cone factor l/{pi ■ n), with pi the momentum carried by 
the marked line. Such integrals appear in a light-cone gauge calculation of the one-loop 
splitting amplitude, using the Feynman diagrams of figure ^. We shall see that there is a 
big difference between the two types of integrals in figure ^; one is properly regulated by 
dimensional regularization without any additional prescriptions, and one is not |79(| . 
The explicit expressions for the two integrals are, 

jia)(^ S) = -i f I fS 6) 

^ ' ^ J TT^/'^ p'^{p — kiy{p — ki — k2)'^ p ■ n ^ ^ ^ 

where s = {ki + ^2)^ = 2ki ■ k2 and, for consistency with the collinear kinematics, z is 
defined by 

ki ■ n = z {ki + /c2) ■ n , k2 ■ n = {1 — z) {ki + k2) ■ n . (3.8) 

For j(°)(z,s) the light -cone denominator involves the loop momentum between the two 
massless legs, while for J^^\z,s) the light-cone denominator involves a loop momentum 
adjacent to the merged off-shell external leg. 
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In order to evaluate these integrals it is helpful to make use of their properties as n is 
rescaled [^], which imply that 

(-s)^+^ [ki +k2) -n 
■^^'^ ^) = , M+e(t^u, (^) ' (3-9) 

(-S)^+^ (fcl + fc2) • 

where /^"\ J*-^^ are functions to be determined. Without loss of generality we may set 

{ki + k2)-n =-s. (3.10) 



At the end of the evaluation we can use eq. ( |3.9|) to replace one factor of l/(— s) with 
l/[(A;i + k2) ■ n]. (The splitting amplitudes are independent of n, so the light-cone vector 
n must cancel from final expressions.) 

One might be tempted to switch to light-cone coordinates, as is commonly done when 
performing light-cone gauge calculations. In these coordinates we take = p ■ n and 
p" = p ■ n* where n* is dual to n, i.e. 

n^ = (no,n), n* = (no,-n). (3.11) 

In light-cone coordinates, the integral J^^\z,s), for example, is given by 

^"'('.^) = J ,2 _ ■ (3.12) 

which is ill-defined because of the unregulated longitudinal integral / dp^ /p~^ . However, 
with this choice of coordinates we have switched to a version of dimensional regularization 
where only the transverse coordinates are regulated. This does not correspond to covariant 
dimensional regularization: dp'^dp~d^~'^p d^p. This difficulty is relatively minor, and 
may be dealt with by reverting to covariant dimensional regularization. Alternatively, we 
could introduce the principal-value prescription, to justify intermediate steps, but then 
remove it prior to performing the final Feynman-parameter integration (see for example 
eq. (34) of ref. ||7^). In either case, care is required because of the ill-defined nature of 
expressions. 

Even covariant dimensional regularization, however, does not suffice to properly reg- 
ulate all light-cone integrals. Let us first compare the behavior of J*^"^ and J^^^ in their 



momentum-space forms (3.6) and (^.Tj). In momentum space, singularities arise whenever 
two or three denominator factors vanish. The light-cone denominator vanishes when p 
becomes proportional to n, p ^ C,n. In this region p^ also vanishes, but the other two 
denominators do not vanish. The singularity that arises here thus looks very much like the 
collinear singularity that arises when p becomes collinear with fci, p — > C,ki. This singular- 
ity is regulated by covariant dimensional regularization, and so we may expect the same to 
be true for the new singularity that arises in the presence of the light-cone denominator. 
This will indeed turn out to be the case. In contrast, for J^"'\ after shifting the momentum 
p' = p — ki we can see that in the soft region — > we have four vanishing denominators, 
which is indicative of the difficulties that will be encountered in evaluating this integral. 
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We can see the difficulty with J^"^ more explicitly using its Feynman-parametrized 
form. We begin by Schwinger parametrizing, 



lldti 

1=1 



vr 



P 

2-r 



{h +t2 + ts)/ - 2p ■ {t2h + tsh + t3k2) 



— t^p ■ n + + t/^ki ■ n 
To integrate out the loop momentum we perform the shift 



, , t2ki + hki + t3k2 , t4n 
P = P + + ^ , 



(3.13) 



(3.14) 



where T = ti + t2 + 1^. (Note that the Schwinger parameter associated with the light 
cone denominator is absent from T.) Wick rotating and then integrating out the shifted 
loop momentum p' gives 



J('^)(z,, 



i=l 



{t2 + tz)tzs U{t2ki + tzki + t^k2) ■ n 



T 



+ t^s + t^ki ■ n 



-'^ 1=1 



htss 5*4(^3(1 - z) - hz) 



+ 



T 



(3.15) 



where we have used the replacements (|3.8| ) and ( 3.10 ) to obtain the last line. The integral 
over the Schwinger parameter associated with the light-cone denominator is now trivial, 
and yields 



J(«)(^,, 



n dU T- 



D/2+1 



i=l 



sihz - tsil - z)) 



exp 



(3.16) 



As usual we may convert the Schwinger parameters to Feynman parameters by defining 
Ui = ti/T and integrating out the overall scale T, yielding a compact Feynman parameter 
representation, 



jW(z,s) = -r(l + e)(-s) 



-2-e 



1 3 3 

dai 5^1 — aj 



-l-e 



i=l 



aiz — 03(1 — z) 



(3.17) 



The reader will observe that for time-like kinematics, the integrand blows up inside the 
region of integration. For example, if z = 1/2 there is a singularity at ai = 03. This 
singularity could be regulated by analytically continuing in z; in the space-like region, 
where z > 1, it is in fact absent. This is not, however, the only singularity in the integrand. 
It is also singular in the corner of the integration region where 02 —> 1, as can be made 
manifest by changing variables 02 = 1 — f , ai = vu, 03 = v{l — u), with jacobian v, 



j(-){z,s)=Til + e)i-sy 



2-e 



dvdu v' 



-2-2e 



Kl-n)]- 



-l-e 



I — U — Z 



(3.18) 
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The singularity as ?; ^ is stronger than the v^^^^ that would correspond to a logarithmic 
divergence, and would give rise to a pole in e. Formally, analytic continuation in e will 
regulate this divergence (in fact the v integral will not give rise to a pole at all), but this 
requires a large analytic continuation, and effectively happens through the subtraction of 
an infinite constant. (There is of course no associated bare coupling here into which such 
an infinite constant could be absorbed.) Note that this pathology is independent of z, and 
so cannot be cured by analytically continuing in the latter variable. 

The principal- value prescription ( |3.5D is a widely-used method to deal with this prob- 
lem. For this integral it would give rise to a 'naked' 1/5 singularity. Absent such a 
prescription, which as discussed above we must avoid for other reasons, care would be re- 
quired in a complicated calculation to ensure that all integrals are continued in a consistent 
manner, and that these continuations do not violate any symmetries. It would be much 
simpler if we do not have to confront this issue at all. As we shall see, our approach to the 
calculation indeed allows us to avoid integrals like J^"^ altogether. 

In contrast, the integral J^^^ in figure |6|(b) is properly regulated by dimensional reg- 
ularization. To see this explicitly, follow similar steps as in the computation of J^"^ to 
obtain, 

„i 3 



jW(z,.) = r(l + e)(-.)-2- / f\da,6(l-TaM-^ 



a2Z + as 



(3.19) 



In this case, the integrand diverges no worse than a~'^~'' near any boundary, and so di- 
mensional regularization renders the integral finite for small negative e without infinite 



subtractions. Indeed, this integral is well-defined and its value is |43|| , 

.) = 2 I(l±lHlzi) (-.)--.i|-L _ i ,„ . , i . . Li.,1 - .)} + 

+ 0(e), (3.20) 

so there is no need for an additional prescription here. We will also encounter the integral 
— z). In the space-like case, with 1 — z < 0, the latter's integrand will be singular 
inside the region of integration. As explained above, this can be regulated by analytic 
continuation in z from z < 1. 

Although it is not necessary, it is still possible to use an additional prescription for 
dealing with the light-cone denominator singularity in J^^\ Had we, for example, used the 
principal-value prescription we would have obtained instead 

s) = r(i + / n«».<i - E «.) . 

i — 1 j 

(3.21) 

following the discussion in e.g. ref. For 5^0 this integral evaluates to 





6 




6 






- 2 In - In 






s 




s 



In z + 

+ ln^z + Li2(l-z) + 0(e). (3.22) 
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Compared with the result in dimensional regularization ( 3.20 ), the 1/e singularity arising 
from the light-cone denominator has been traded for a In (5 singularity. The In 6 would 
then appear in the result for the splitting amplitude. As mentioned earlier, such a result 
cannot match the splitting amplitude describing the collinear limits of gauge-invariant, 
dimensionally-regulated one-loop amplitudes, because the latter depend only on e, not 6. 
We therefore wish to avoid additional prescriptions, which are in any event unnecessary 
for J^''-' . They would be required only for integrals like J^"^ . 

What distinguishes the integrals (a) and (b) in figure ^7 Could it be that the appear- 
ance of ill-defined integrals like (a) is purely an artifact of the gauge choice and would not 
appear in a more physical construction of the splitting amplitudes? We will see that the 
answer to the latter question is yes. 

A general difficulty with diagrammatic approaches 
is that gauge-invariant results are chopped up and sep- 
arated into gauge-dependent pieces by the decompo- 
sition into diagrams. For example, for a fully on-shell 
scattering amplitude calculation at tree level, each 
light-cone gauge Feynman diagram contains light-cone 
denominators. Yet by gauge invariance their sum must 
be free of such denominators. Unitarity then implies 
that similar cancellations should happen at loop level. 
For splitting amplitudes, one leg is off-shell. Hence 
the cancellation of light-cone denominators is not com- 
plete; nevertheless, there will be a partial cancellation. 
In order to understand which light-cone denominators 




I 

I 



Figure 7: The two-particle cut 
of a one-loop splitting amplitude. 
The cut is represented by the 
dashed line. On the left-hand side 
all legs, including the cut ones, 
are on-shell. On the right-hand 
side the merged leg represented 
by a thick line is slightly off-shell. 



may appear and which should not, we use unitarity. 
On dimensional grounds, the result for the splitting amplitude must have the form 



fie,z) (- 



(3.23) 



for some f{€,z), where the epsilonic power of (— s) follows from the integration measure 
^4-2ep_ In the time-like region, s > 0, we have (— s)"*^ = I — eln(— s) -|- . . . = 1 — e(lns — 
ivr) -|- . . ., so any desired order in the Laurent expansion in e of /(e, z) can be computed by 
extracting the absorptive part of the function to one higher order in e. If we can compute 
the absorptive part, via unitarity cuts, to all orders in e, we can completely determine 
f{e,z). The two-particle cut of the one-loop splitting amplitude is depicted in figure |^. 



Using the Cutkosky rules |8C] we obtain. 



Absp Split^^\z;a,b) 



(3.24) 



5+{ll)5+{ll)Af{a,h,ei,-e^) X d^„(£i)4^(£2) x V^^'t {1^,-1^) 



(2vr) 



D-2 



where ii = p, ^2 = P — ka — k^, the vertex V°'^'^ is defined in eq. (|3.2[) , d^^, is the physical 
state projector given in eq. (p.4|), and S~^{£'^) = Q{£^)S{£'^). This is the only non-trivial cut 
of the one-loop splitting amplitude and therefore yields the complete absorptive part. On 
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(a) (b) 

Figure 8: Sample three-point integrals at two loops containing light-cone denominators. Integral 
(a) is ill-defined without additional prescriptions, but integral (b) is rendered finite by dimensional 
regularization alone. 



the left-hand side of the cut, ^4^^^ is a gauge-invariant amplitude — all legs including the cut 
ones are on-shell. Any light-cone denominator appearing to the left of the cut must there- 
fore be spurious: gauge invariance dictates that light-cone denominators cannot appear 
after combining all Feynman diagrams, because such denominators would not appear in a 
covariant gauge. Because the full one-loop splitting amplitude can be reconstructed from 
its absorptive part, the light-cone denominator appearing in the integral in figure |6|(a) must 
also be spurious (absent from the sum over all diagrams) in the full calculation, not just 
the absorptive part. In contrast, light-cone denominators on the cut lines are introduced 
by the Cutkosky rules via the sum over polarizations across the cut line, 

E (P) = -V,u + ^""^ ^J'"^'^ = -V ■ (3.25) 

a p n 

This physical state projector is the on-shell version of the one appearing in the light-cone 
gauge propagator ( |3.3D . The light-cone denominators that appear here can survive. This 
is the type of light-cone denominator that appears in the integral of figure |6|(b). (Indeed, 
as noted earlier, some light-cone denominators must survive in order to get an answer of 
the sufficient polylogarithmic complexity.) 

The unitarity argument indicates that only light-cone denominators associated with a 
cut line need survive. (In the multi-loop case, denominators associated with lines to the 
right of all cuts can also survive.) Thus it should be possible to perform calculations where 
dangerous light-cone denominators of the sort depicted in figure ^a) do not appear. The 
unitarity-based sewing method, which we present in section ^, has exactly this property. 

3.4 Two-Loop Light-Cone Integrals 

At two loops we encounter a similar situation: some integrals are properly regulated solely 
by covariant dimensional regularization, and others are not. 

As concrete examples consider the two-loop integrals in figure ^, given by 

I (3 26) 

^ ' ' J {27r)D {27r)D p\p - kif{p + qf{q + ki+ k2f{q + k^f{p - ki) ■ n ' ^ ^ 

iWf^ s)=^^ I (3 

^'^ J {2tt)D {27r)D p2(p-kiy{p + q)^{q + ki+k2y{q + kiyp-n- ^ 
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Figure 9: The three-particle cut of a two-loop splitting amplitude. On the left-hand side of the 
cut, the amplitude is gauge invariant since all legs including the cut ones are on-shell. 



In our calculation of the two-loop splitting amplitude from the unitarity sewing method, we 
encounter only integrals similar to L^^\ Following similar steps as at one loop, we obtain 
the Feynman-parameterized form, 

L(")(z,s) = r(l + 2e)(-s)-2-2^ 

X ['f]daa(l-ya,)^, ^ (3.28) 

Jo V ^ ^ z[ae{a2 + as) + a3{a2 + a4)\ - a2ae 

where 

A = (a2 + a4)(a3 + 05) + 06(^2 + 03 + ^4 + 05)- (3.29) 

This integral has an insufficiently regulated divergence in the region 05 — > 1. To see this, 
make the change of variables 05 = 1 — v, 0^=2,3,4,6 = vhj, for which the jacobian is v^: 

L^''\z,s) = r(l + 2e)(-s)-2"2^ 

dv n db, 6(1 -Vfe,).^^-^- ^ .. .(3-30) 

i=2;3 4,6 j ^[^6(62 + 63) +63(^2 + 64)] -^2^6 

where 

Afe = vb3{b2 + 64) + vbe{b2 + 63 + h) + (1 - v){b2 + 64 + be). (3.31) 

Like J(") of the previous subsection, the integral has a power-law divergence for f ~ 
which would require a large analytic continuation, or equivalently the subtraction of an 
infinite constant. This divergence is again independent of z, and hence cannot be cured by 
analytic continuation in that variable. 

On the other hand, for integral (b) in figure]^ we have the Feynman parametrized form 

L^'Hz,s) = -r(l + 2e)(-s)-2-2^ 

. rflda.5(l-Y^a,) ^"^"^"^"f ^. (3.32) 



0206 -I- z[0406 -I- 05(02 + 04 -I- Oe)] 

In this case the integral is well-defined; at all boundaries, the integrand goes like v~^~"^'', 
m = 1, . . . , 4. For example, as 03 = 1 — — > 1, the presence of 03 in the numerator lessens 
the singularity to v~^~'^. 

But need we concern ourselves with the possible appearance of ill-defined integrals 
like L^'^)? As at one loop, light-cone denominators at two loops can be separated into 




Figure 11: Additional three-point integrals at two loops containing light-cone denominators. Both 
integrals are well-defined using dimensional regularization alone. 

two categories, depending on whether they appear in unitarity cuts or not. For example, 
consider the three-particle cut of a two-loop splitting amplitude shown in figure |^ The 
three-particle cuts of the two integrals in figure ^ are shown in figure |l^. Since all legs of 
the five-point amplitude on the left-hand side of the cut, including the cut ones, are fully 
on-shell, then following the same logic as in the one-loop case, the light-cone denominator 
appearing in the integral in figure |8|( a) is a light-cone gauge artifact which can be eliminated. 
(This type of argument cannot be used on the right-hand side of the cut, because the merged 
leg is off-shell.) The light-cone denominator appearing in figure |8|(b), on the other hand, 
is allowed because it corresponds to a physical-state projector on a cut line. 

The use of the unitarity-based sewing method, which we describe in detail in the fol- 
lowing section, allows us to avoid the use of any prescription for light-cone denominators 
in the calculations in following sections. But one can also imagine applying the insights 
above to a more standard diagrammatic calculation in light-cone gauge. One could proceed 
as follows: introduce one of the standard prescriptions for dealing with light-cone denom- 
inator singularities. Then, attempt to combine diagrams algebraically to remove those 
light-cone denominators which by unitarity cannot appear in the desired quantity. Once 
all singularities leading to ill-defined integrals have canceled, one can remove the additional 
prescription (for example, by taking 5 — > in the PV prescription), and only then perform 
the loop integrals. 

The unitarity arguments above can be applied not only to ill-defined integrals, but 
also to rule out certain well-defined integrals. As an example, consider the two integrals in 
figure 0. Both turn out to be well-defined (for integral (a), this is clear from momentum- 
space power-counting). If we examine their cuts, shown in figure however, we see that 
integral (a) has a light-cone denominator to the left of the cut, and hence cannot appear. 
Integral (b) has a light-cone denominator on the cut, and hence is not ruled out by the 
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(a) (b) 

Figure 12: The three-particle cuts of the two-loop three-point integrals of figure [ll|. 
unitarity argument. Indeed, it is the master integral Wedge(2;, s) of figure^, whose expHcit 



expansion in e is given in eq. (6.49) 



In summary, the unitarity cuts point to a method that sidesteps the prescription issues 
associated with light-cone denominators, because only a restricted set of integrals appear. 
In the next section we explain in some detail how to construct loop-momentum integrands 
using the unitarity method. In our calculation of two-loop splitting amplitudes in sections ^ 
and ^, such loop integrals are reduced to a linear combination of master integrals. There are 
several equivalent bases that are convenient for different aspects of the calculation; one of 
the equivalent forms contains only integrals well-defined using dimensional regularization 
alone, and not requiring any additional prescriptions. All singular boundaries lead to 
logarithmic singularities in the e — > limit, just as for the integral J^*^ discussed in the 
previous subsection, or the integral L^^^ above. 

4. Review of the Unitarity-Based Sewing Method 
4.1 Overview 

The unitarity of the scattering matrix in a quantum field theory is the statement that 
probability is conserved. It is an essential property of any sensible and consistent theory. 
It relates the non-forward part T of the scattering matrix S to its square, —i{T—T^) = T^T, 
where T is defined via 5 = 1-1- iT. In Feynman diagrams, unitarity is expressed by the 
Cutkosky rules |80, which express the 'imaginary' or absorptive part of a diagram^ in 
terms of phase-space integrals over products of lower-loop diagrams. The product is given 
by 'cutting', replacing specified sets of propagators by delta functions in the propagator 
momentum. Loop amplitudes are computed, of course, by summing over appropriate 
collections of Feynman diagrams. Their absorptive parts are given by sums of products 
of lower-loop diagrams. Collecting all diagrams on each side of a cut into amplitudes, we 
see that the absorptive parts of loop amplitudes are just sums over products of lower- loop 
amplitudes. 

This observation is particularly powerful in gauge theories (and in gravity as well). In 
gauge theories, there are extensive cancellations between different diagrams in the com- 
putation of scattering amplitudes for on-shell states. These simplifications can be made 
manifest at early stages of a tree-level calculation using the spinor helicity representa- 
tion for gauge-boson polarization vectors. (These simplifications may be understood using 



^By 'imaginary' we mean the discontinuities across branch cuts. 
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twistor space Q.) The final answers in massless theories are particularly simple, some- 
times simpler than the expression for a single Feynman diagram out of the hundreds or 
thousands that contribute, and have a natural expression in terms of spinor products. We 
may then express the cut of an on-shell one-loop amplitude, given by a product of on-shell 
tree amplitudes (or a sum of such products), in simple form as well. This simplicity carries 
through order-by-order in perturbation theory. The sewing technique aims to exploit this 
simplicity, by turning the process around, and building loop amplitudes out of their cuts, 
in turn given by lower-loop amplitudes. 

The full amplitude can in principle be reconstructed using dispersion relations. For 
general gauge theories in four dimensions, the dispersive reconstruction of an amplitude 
suffers from an additive ambiguity related to divergent ultraviolet behavior. One can 
add a rational function, free of cuts, to the amplitude. This problem has traditionally 
hampered the use of dispersion relations to obtain complete amplitudes. It is solved in 
massless theories^ through the use of dimensional regularization, which effectively tames 
the ultraviolet behavior of the bare integrand p^ ] and thereby removes the need for explicit 
subtractions. This represents a third role for the dimensional regulator beyond its usual 
roles as a regulator for ultraviolet and infrared divergences. The sewing technique we review 
is equivalent to the use of dispersion relations in dimensional regularization, although for 
practical purposes it is preferable to make use of ordinary Feynman-integral techniques for 
performing the necessary integrations rather than doing explicit dispersion integrals. (At 
one loop, for example, knowledge of the complete decomposition of n-point integrals in 
dimensional regularization in terms of a basis of known integrals |^3| reduces the problem 
to an algebraic one.) Sewing back together cut amplitudes, with the cut lines on shell but 
treated exactly in D = 4 — 2e dimensions, will reproduce the full gauge-theory answer. 
(The different ways of continuing the amplitude to D dimensions correspond to the use of 
different variants of dimensional regularization, such as CDR HV |6^, or FDH ||8|, [7C|].) 
One can then expand in e to obtain the answer through 0{e^), including the rational terms. 
A more pedestrian way to understand how the rational terms are included properly is to 
observe that in dimensional regularization, these terms are not purely rational, but rather 
are rational functions of the momentum invariants, multiplied by (— s)^*^ at L loops, for 
some invariant s {e.g. as in eq. ( 3.23| )). At ©(e), this factor contains a logarithm, and 
hence an imaginary part for s > 0. Only when the e ^ limit is taken at the end of the 
calculation, does the term become purely a rational function, free of discontinuities. If one 
calculates the cuts to all orders in e, the full amplitude can also be reconstructed to all 
orders in e. 

Yet another way to understand the uniqueness of cut reconstruction is as follows: In 
reconstructing a loop amplitude from a single cut, there is always an ambiguity with respect 
to adding terms in the numerator of the loop-momentum integrand which vanish on the 
cut in question, e.g. those terms proportional to l\ and I2 where £1 and £2 are momenta 
crossing a two-particle cut. However, after merging the information from all cuts, the 



^In massive theories, there is an additional source of ambiguities (from masses inside bubbles on ex- 
ternal legs). When there is only one mass in a calculation, this problem can be resolved through simple 



adjustments [Q. 
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(b) (c) 

Figure 13: Color-ordered Feynman diagrams for tlie one-loop four-point amplitude in Trcfr' field 
theory: (a) the box diagram (b) the s-channcl bubble diagram (c) the s-channel triangle diagrams. 
The two t-channel triangle diagrams and the i-channel bubble diagram are not shown explicitly. 
Bubbles on external legs vanish in dimensional rcgularization and are not shown here either. 

only remaining ambiguities are those terms which are free of cuts in every channel. For 
the all-massless theories we consider, such integrals always take the form of external leg 
corrections; that is, they are scale-free integrals which vanish identically in dimensional 
regularization. 

4.2 Sewing at One Loop 

Consider first the sewing method at one loop. Before explaining it in generality, it will be 
useful to examine the procedure in a simple example. In each example here and in later 
subsections, we will make contact with standard methods by starting with an amplitude 
expressed in terms of conventional Feynman diagrams. In the examples, we will work 
in a massless Tr cj)^ field theory (with (p transforming under the adjoint of SU(A'^c))) but 
as discussed above the method applies to general theories, and indeed is relatively more 
powerful precisely in field theories with many redundant variables in their covariant form, 
such as gauge theories and gravity. 

Let us start by considering the four-point one-loop amplitude. The full amplitude has 
Bose symmetry, which just as in the gauge-theory case we can exhibit most concisely by 
rewriting it as a sum over color permutations of a more basic quantity. The leading-color 
contributions (leading in a 1/Nc expansion) involve only planar diagrams. We will focus on 
these contributions in the examples. In particular, we focus on the coefficient of the color 
trace Tr(r"i T^^T'^^T"'*) given by color-ordered diagrams ||6^, |l^ with the 1234 ordering 
of legs. 

There are seven color-ordered diagrams contributing to this one-loop partial amplitude. 
The ones with cuts in the s channel are depicted in figure The ordered amplitude has 
cuts in only two channels, s and t. If we examine the s channel, we see that only four of the 
diagrams contribute to the cut: the box, two triangles, and one of the bubbles. Similarly, 
four diagrams contribute to the cut in the t channel. The s-channel cut may be obtained 
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Figure 14: Diagrams depicting the s-channel cut of the one-loop four-point amphtude in Ti cjr' 
field theory. 



by replacing the propagators cut in figure 14 via 

1 



(4.1) 



-|- le 

This replacement converts the loop integral to one over the phase space of the two cut legs, 
which are placed on shell. We can also see that the sum of terms factors, so that on each 



side of the cut we obtain a tree amplitude as the sum of diagrams, as shown in figure |15 



In each channel, the cut is thus given by a phase space integral of the product of two tree 
amplitudes, 



A«(l,2,3,4) 



(27r)^-2 



-6+{f^)6+{£l)d^{£i + £2 + ki + k2) X 
x^(0)(i^2,£2,^iM^°^Mi,-^2,3,4), 



where 



A(°)(1,2,3,4) 



\Si2 SiiJ 



(4.2) 



and we have suppressed powers of the three-scalar coupling. 

The sewing procedure reverses this process. We start, for example, in the s channel. 
Multiply the tree amplitude on the left-hand side of figure ^ by that on the right-hand 
side, 

(1, 2, £2, ^i)^(°) (-^1, -^2, 3, 4) = -(— + —) (i (4.3) 

where ii = —£.2 — — ^2- Put in the two propagators crossing the cut, and integrate over 
the loop momentum £ = £1, to yield. 



d 



D, 



1 1 \ 1 / 1 

+ 



1 



(27r)^ Vsi2 sun^KsM SiJ{£ + ki + k2Y 



d^£ 



+ 



sl^£^[£ + ki + k2f si2£^{£-kiY{£ + ki + k2Y 
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(4.4) 
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Similarly, from the t-channel cut, we obtain. 



/ 



1 



1 



£2(i + k2 + k3y 



+ 



1 



+ 



+ 



1 



1 



+ 



1 



(4.5) 



S23{i + k2)H^{e + k2 + k3y {l + k2?P{l-kif{l + k2 + k3f_ 

The first three terms have no cut in the s channel, but the last term does: it is given by 
the residue of the poles as ^ — A;i and £ + ^2 simultaneously go on shell, 

1 



(4.6) 



Alternatively, we can shift 



2{£ - ki) ■ ki2{(i - ki) ■ ki 

■ I + ki, upon which the last term becomes identical to the 
last term in eq. (4.4). 

We cannot simply add the contributions from the s and t channels, because this would 
correspond to the sum of eight diagrams, double-counting the box diagram ^(a), given 
by the last terms in eqs. (^^) and ( |4.5| ). Accordingly, we must find a function which has 
the correct cuts in all channels. We can do this either before or after integration, although 
in general it is easier to do it before loop integration. One way is simply to sum both 
contributions, and then remove the overlap: terms in one cut channel which also have a 
'cut' — in the sense of having the propagators which give rise to a cut — in the other 
channel. One can alternatively think of this as 'merging' the two expressions, taking a 
term if present in either cut or in both, but taking it only once in the latter case. 

The net effect is to drop one of the two equivalent terms; we obtain the sum of the 
remaining terms for the ordered one-loop amplitude, 

^«(1,2,3,4) = 



d'^l 
(27r)^ 



+ 



'12 



i^{t + ki + k2Y Si2p{l-kiY{l + ki + k2)' 



+ 




Figure 15: The s-channel cut of the one- loop four-point amplitude in Trc/)^ field theory, seen as a 
product of tree amplitudes. 
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si2(,i + h)H^{£ + h + k2f ii + ki)H^e-hy{e + ki + k2)'^ 
111 1 



+ 



+ 



i^{e + k2 + k^Y S23 - kxf{t + k2 + ksY 
1 1 



(4.7) 



+ 



{i + k2)H^{e + k2 + k3Yy 



exactly as would have emerged from a Feynman-diagram computation. Of course, in a (p^ 
field theory, there are no cancellations between different diagrams, so the sewing method 
is also equivalent in complexity to the usual approach. In gauge theories, the sewn on-shcll 
tree amplitudes are much simpler objects than the one-loop diagrams, and so the sewing 
approach helps minimize the complexity of intermediate steps. 

In the above examples, the procedures for sewing and removing any overlaps or double- 
counting are completely mechanical. Note that none of them make any reference (to use 
an old-fashioned language) to double dispersion relations. Indeed, only in an abstract 
sense are dispersion relations used at all, since we do not perform the dispersion integrals 
explicitly, but rather implicitly via construction of appropriate Feynman integrals. We will 
next explain how to formalize these procedures, and then give an algorithm which can be 
used to implement them in practice. 



To formalize the sewing procedure, introduce the basic promotion operator J A. It 
will be applied to products of amplitudes (or to terms from a product). It represents the 
combined operations of summing over hclicity states (and over different particle states if ap- 
propriate), re-expressing spinor products in terms of the cut momentum, multiplication by 
the cut-crossing propagators, and completion of dot products in denominators to standard 
propagator denominators. It does not introduce the phase-space integral over the product 
of amplitudes. The result of the promotion operation is an integrand which depends on the 
external momenta and on the loop momentum. Note that the sum over intermediate states 
must in general be carried out in Z) = 4 — 2e dimensions, and that it is a sum only over 
physical states. This implicitly introduces a physical projection operator. In calculations 
of full amplitudes, the resulting operators leave little trace, but in calculations of splitting 
amplitudes such as the one we carry out in the present paper, these operators will give rise 
to light-cone-like denominators in integrals. 

We will also need to introduce the cut-projection operator Vs, which yields the part 
of its argument that has a cut in the s channel, where s denotes an arbitrary invariant 
of m consecutive external momenta. At one loop, it extracts the joint pole term in two 
propagator denominators jrp, where (^1-1-^2)^ = s. It corresponds to requiring that a pair 



(any pair) of propagators yielding a cut in the s channel be present in the diagram. 

In a computation with n external massless momenta, there are in general n{n — 3)/2 
independent invariants in D dimensions. Denote by C the ordered set of these invariants. 



(We will denote the number of elements in C by nc-) We will sew the channels in the 
specified order, with the notation Sj G C denoting the j-th invariant in C. The optimal 



C = {S12,S23,S123,---}- 



(4.8) 
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ordering (from the viewpoint of computational efficiency) depends on the process and the 
particle content of the theory, but of course the final answer is independent of this ordering. 
Let Kj = + ■ ■ ■ be the momentum whose square is the given invariant sj . The first 
momentum (within the cyclic order of external momenta) we will denote aj, and the last 
momentum by bj. The momentum before aj will be labeled aj — 1, and the one after bj, 
bj + 1. 

The analytic behavior in different invariants is independent even if the invariants are 
related by Gram determinant conditions arising from the restriction to four dimensions. 
Thus even if we take all external momenta to be in four dimensions, we must still take the 
full set of £)-dimensional invariants. 

The complete integrand of the one-loop amplitude ^^^^(1, . . . ,n) is then given by the 
sum over all channels. 



X ^(°) {£ + Kj,bj + 1,..., aj - 1, -£) . (4.9) 



That is, we sum over all channels, each time removing all terms already found in previous 
channels. The amplitudes in this equation must in general have the sewn legs {i or i + Kj) 
in D dimensions. Whether the external legs are taken to be in four dimensions or in D 
dimensions depends on the variant of dimensional regularization employed. In practice, it 
is best to use a four-dimensional scheme, and convert later if necessary. 

Indeed, there arc several practical aspects not addressed by the formal expression 
above. These include questions of diagram labeling, classification, and the use of a basis 
for organizing numerators of terms in the integrand. As the formal expression hints, none 
of these tools are intrinsically required by the unitarity-based sewing method. The cut 
projection can be performed by extracting residues of poles; and the promotion involves 
simple algebraic manipulations. Furthermore, non-manifestly vanishing expressions still 
vanish, and do not affect the final answer. In a practical calculation, however, wc would 
like the cut projection to be simple, ideally just amounting to the identification of the 
formal coefficient of a pole. We would like to avoid the appearance of complicated expres- 
sions which actually vanish. Furthermore, the integrand produced by the sewing method 
will ultimately be fed to an integration machinery which does require the corresponding 
diagrams to be labeled and classified by topology. We therefore might as well incorporate 
these aspects into an algorithm. 

To do so, start with all color-ordered graphs with a maximal number of propagators 
containing the loop momenta. We will call these 'parent' diagrams. All other topologies can 
be obtained by canceling propagators, that is by multiplying the numerator by an inverse 
propagator. These we will call 'daughter' diagrams; below, we will also include the parent 
in its set of all daughters. At one loop, if all external particles are massless, there is in fact 
only one parent diagram. (Recall that we arc restricting attention to processes where all 
internal masses vanish.) If we have two or more massive external particles (for example, 
W bosons), we will have different parent diagrams corresponding to the different ways of 
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attaching the massive legs to the loop. At two loops, there are n{n + l)/2 planar parent 
diagrams when all external particles are massless. (Integrals that reduce to products of 
one- loop integrals, such as bowtie integrals, have intrinsically two- loop integrals as parent 
diagrams.) Examples are shown in figure 

Each cut will in general start with a different labeling of any given parent diagram, 
because the loop momentum i may denote a different propagator. The algorithm will make 
use of a simple relabeling operation, reviewed in appendix to bring these into canonical 
form. The labeling of propagators must incorporate a notation for the parent diagram, 
because the algebraic relations of dot products to inverse propagators differ from diagram 
to diagram. 

The external gluon legs we may choose to treat using formal polarization vectors £i{ki), 
or using the spinor-helicity method. The external fermion legs we may choose to treat using 
formal spinor wavefunctions u{ki), u{ki), or again using a helicity basis. Either choice 
(or a mixture) may be employed with the algorithm we will present below. The basis 
one should use for expressing numerator polynomials depends slightly on the external leg 
treatment, because there can be different numbers of independent invariants from which 
the polynomials are built. The algebraic processing of expressions will also be somewhat 
different. In all cases, the basis at one loop will contain all inverse propagators containing 
the loop momentum. These are sufficient to express all dot products of the loop momentum 
with external momenta. Note that Levi-Civita tensors involving the loop momentum can be 
converted to Gram determinants (and thence to dot products) by multiplying by another 
Levi-Civita tensor involving only external momenta. (The latter object is just another 
constant as concerns the manipulations we perform.) We will also need the square of 
(— 2e)-dimensional components of the loop momentum, (^_2e)^. We must add dot products 
of the loop momentum with formal external polarization vectors (if any). These will not 
give rise to expressions that can cancel propagator denominators. (We could have taken 
these to be in D dimensions, after all, in which case they would clearly be independent.) In 
all algebraic manipulations, one should be sure to use momentum conservation, eliminating 
one external momentum, and re-expressing invariants in terms of an independent set, in 
order to avoid the appearance of zero in obscure forms. (For formal expressions £j ■ ki, one 
should pick a momentum other than kj to eliminate, so as to impose the on-shell conditions 
too.) 




(a) (b) (c) (d) 

Figure 16: Examples of parent and daughter diagrams: (a) a parent diagram at one loop; (b) 
a daughter diagram at one loop, which is a daughter of the diagram in (a); (c) and (d) parent 
diagrams at two loops. 
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For example, in computing the an n-point gluon amplitude, we can pick the standard 
labeling to have the loop momentum between legs n and 1. If we treat all external legs in 
the spinor-helicity basis, then the basis set will simply be, 

{f, {£ - hf, {i-h- k2)\ ...,{i-h K-lf, {i-2ef}; (4.10) 

if we choose to treat legs 1, . . . ,j using formal polarization vectors, we should add 

{ei ■£,..., ej-£} (4.11) 

to this set. 

Amplitudes with fermions in a loop arise from sewing amplitudes with pairs of external 
fermions on either side of a cut. Sewing will include closing a fermion loop. (The usual 
minus sign must be included explicitly.) This yields a spinor trace, which can be expanded 
in terms of dot products and Levi-Civita tensors. The latter can be converted to dot 
products as described above. Internal fermions thus do not require any new basis elements. 

Amplitudes with external fermions will contain in different terms, factors of a 'spinor 
string' consisting of an external spinor wavefunction (either formal u, u or in the helicity 
basis), a product of gamma matrices dotted into various D- or four- vectors (momenta, 
polarization vectors, the light-cone vector, or other spinor strings), and an ending spinor 
wavefunction. Roughly speaking, we need to perform sufficient manipulations on these to 
ensure that no difference of two such objects contains a factor of an inverse propagator. 
A basis for spinor strings involving spinor wavefunctions can be obtained by commuting 
loop momenta to the left; commuting /_2e) if present, to the next position; and commuting 
formal polarization vectors (if present) to an ordered sequence following them. The spinor 
string can end with either a spinor carrying an external momentum or the light-cone vector, 
or another external fermion wavefunction. Alternatively, one can convert the spinor string 
to a trace by multiplying by appropriate spinorial factors involving only external momenta 
and spins, and then expanding the trace into dot products and Levi-Civita tensors as above. 
(The spinorial factors are just spinor products in a helicity basis.) In this case, no spinor 
strings are needed in the basis. 



The projection Vsj of the promoted integrand in eq. (O) back onto the same channel 
which was sewn, ensures that no (spurious) terms lacking a cut in the sewn channel are 
generated. This projection is not really needed in the purely formal expression, but ensures 
that algebraic manipulations when working in the basis required for a practical algorithm 
do not create unwanted terms. 

The basis will be used to identify terms that have cuts in different channels. A term 
with uncanceled propagators corresponding to the cut channel will have a cut in that 
channel. 

For a practical algorithm, one may proceed as follows: 

1. Form the ordered set of all independent channels in D dimensions. (That is, in 
determining the independence of different invariants, one should use only momentum 
conservation, and not integer-dimension-specific Gram determinant relations.) 
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Figure 17: Sewing together two tree-level amplitudes to produce a contribution to the one-loop 
integrand. 



2. Associate a labeling of internal lines to each distinct Feynman-integral parent which 
has a daughter appearing in the integrand. (Because of cancellations in gauge theo- 
ries, not every topology that appears in the set of the usual Feynman diagrams for a 
given process will necessarily appear in the sewn integrand.) 

3. For each parent integral, form a basis set for expanding numerators, consisting of the 
inverse propagators; the square of the (— 2e)-dimensional components of the loop mo- 
mentum, (£_2e)^; and dot products of the loop momentum with any formal external 
polarization vectors. For external fermions with formal wavef unctions, spinor strings 
should be added as described above. 

4. Initialize the integrand's value to zero. 

We will form the integrand by iterating over channels. For a color-ordered amplitude, 
each channel corresponds to a consecutive set of external momenta, ka^, ■ ■ ■ , kf,. , the cut 
invariant being sj = {ka^ + . . . -\- kb^)"^. 
For each channel, 

5. Sew. Form the product of the two on-shell tree amplitudes A^^\li,aj, . . . , 6j, ^2) and 

, n, 1, . . . , aj — 1, —^1), where £1 -|- £2 + Kj = 0, summing over the 
products of amplitudes for the different particle types and helicities that can circulate 
in the loop. The sum must in general be performed in Z) = 4 — 2e dimensions. Use 



polarization- vector identities (eq. ( |3.25D , etc.), 



t ■ n 

a 

u{£)u{£) u{-i)u{-i) (4.12) 

to express everything in terms of the cut momentum i. Here, n is the light-cone 
reference vector. (If all external legs are on shell, different may be used for dif- 
ferent cut legs if desired, and indeed the light-cone denominators can be removed 
algebraically.) The dimensionality of ?? — the value of rjji, — depends on the variant 
of dimensional regularization (see eq. ([7.1|)). Multiply by —1 when a fermion loop 
is created by sewing. Put in the propagators crossing the cut, for each cut leg. 



-30- 



Complete dot products in the denominator to form propagator denominators adja- 
cent to the cut, 2i ■ ki —>■ ±(£ it fcj)^, the + corresponding to ki on the left side of 
the cut, the — sign to those on the right side of the cut. Rewrite the expression in 
terms of the basis set, and expand sums so that each term can be classified as the 
daughter of a single parent diagram. As mentioned above, the sewing procedure for 
a given cut is ambiguous with respect to numerator terms which vanish on the cut in 
question, i.e. those terms proportional to if and After merging the different cuts 
(in step 9 below), the only remaining ambiguities will be those which have no cuts in 
any channel. Such integrals have the topology of external leg corrections, and vanish 
for dimensionally-regulated all-massless theories. The sewing operation is depicted 
schematically in figure 17. In eq. ( [4. 91) , it corresponds to the promotion operator J. 



6. Put into canonical form. Use momentum conservation to reduce the number of cut- 
crossing momenta (now loop momenta) appearing in spinor traces, and then expand 
spinor traces and all dot products in terms of the basis. (As explained above, this 
expansion may not always be necessary.) When using the spinor-helicity method for 
external polarization vectors, or explicit helicity states, we will obtain spinor strings of 
the form {j^ ' ' ' ' ' ' is ) • explained above, one can complete these to a trace, 
then convert the trace to dot products and Levi-Civita tensors. The latter should be 
converted to Gram determinants (and thence to dot products) by multiplying (and 
dividing) by a Levi-Civita tensor involving only external momenta. (Spinor strings 
involving only external momenta need not be manipulated, obviously.) 

7. Relabel. For each term in the sewn expression for a given channel, relabel the mo- 
menta to the standard labeling for its parent integral. Where required, insert factors 
of squared momenta in the numerator and denominator to match a 'parent' diagram. 
In some cases, it will be possible to obtain different parents by inserting different 
factors; it doesn't matter which one is picked. 

8. Clean. Remove all terms which have no cut in the current channel. (Such terms 
might have been introduced by earlier algebraic manipulations.) That is, using a 
canonical basis as described above, remove any terms which do not contain both cut 
propagators. This step corresponds to the operator Vg. in eq. (^). 

9. Merge. Remove all terms in the current-channel sewn expression that already appear 
in the net integrand. That is, remove any term which has cuts in a previously- 
processed channel. Using a canonical basis as described above, it suffices to pick 
out and remove terms that have a pair of propagators corresponding to a cut in a 
previously-processed channel. This step corresponds to the operator Hti ~ T^si 
in eq. (O). 



10. Accumulate. Add the remaining terms in the current-channel sewn expression to the 
integrand. This step corresponds to the sum in eq. (^ 

11. Continue with the next channel at step 
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In special cases (for example, massless supersymmetric theories at one loop), it may 
be possible to compute the cuts using four dimensional helicity states, and to make use of 
spinor-helicity simplifications for the cut-crossing momenta (It is always possible 

to use such simplifications for the external momenta at an early stage of the calculation.) In 
this case, one must re-express spinor products involving ii and £2 in terms of dot products 
of these momenta with other vectors, {aii) [£i b] —>■ {a~\ i\ \h~). (This is always possible 
because the cut-crossing momenta will appear with opposite phase weight in the two tree 
amplitudes on either side of the cut.) It may also happen that some channels are redundant; 
this will happen when all integrals appearing in the answer have cuts in multiple channels. 
In the computation of splitting amplitudes, one side of the cut will have an on-shell tree 
amplitude, while the other side contains a tree-level splitting amplitude. The sewing pro- 
cedure will introduce physical-projector denominators. These resemble light-cone gauge 
denominators. In the computation of on-shell amplitudes, these physical-projector denom- 
inators will disappear algebraically. (This is a consequence of gauge invariance.) In the 
computation of splitting amplitudes, in contrast, these denominators will survive into the 
integration, and in fact play a crucial role in obtaining the right sort of integral. However, 
as discussed in section ^, at one loop they will only arise in lines corresponding to cut 
momenta. As we have seen in the example of the 1 — > 2 splitting amplitude, integrals with 
such projectors have no singularities beyond those regulated by dimensional regularization. 

In general, the integrand has cuts in many different channels. Indeed, the resulting loop 
amplitudes are expressed in terms of polylogarithmic functions that have discontinuities 
in several different invariants. Accordingly, any given term in the integrand may have 
combinations of propagators leading to cuts in different channels, and thus might emerge 
from sewing in any of those channels. The different ways of sewing must all yield the same 
answer. The simplest way of seeing this is to consider a gedanken calculation of a one-loop 
amplitude from Feynman diagrams. We can extract the cut of the sum of all diagrams in 
any given channel. Putting the cut legs on shell, the sums of diagrams on either side will 
yield a product of tree amplitudes. Because of their origin in a unique expression, however, 
reconstructions of the analytic functions from different channels will yield the same answer. 

In terms of the cut projection operator introduced earlier, this amounts to the state- 
ment that 

J^i(s2)^2(s2) = V,^\Az{sx)A^{s{). (4.13) 

where Ai^s) is an abbreviation for the amplitude with two adjacent legs satisfying 
2l\ ■ £.2 = s. Checking the consistency of different cuts in intermediate steps, as expressed 
by this equation, is a good way of verifying the correctness of code implementing the 
algorithm. 

The consistency of different cuts also shows that the order of evaluation of channels, 
i.e. the ordering in C in eq. ( |4.8| ), does not affect the final result. For example, if we were 
to evaluate two channels in the order {si, S2}, we would obtain, 

1^3(51)^4(51) + (1 - Vs,)\Ai{s2)A2{s2). (4.14) 
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(Recall that the projection back onto the sewn channel is not really needed in the formal 
expression, and has been omitted here.) In the other order, we would obtain 

\Ai{s2)A2{s2) + (1 - Vs,)\a^{s^)A^{s^). (4.15) 
The difference of the two evaluations is 

Vs,\A3isi)A4isi) - Vs,\Ai{s2)A2{s2), (4.16) 



which vanishes using eq. ( 4.13 ) 



Many terms (all, in an AA = 4 supersymmetric theory) may have cuts in more than one 
channel. The cut consistency condition can be used to cross-check these terms, and is often 
of great utility in debugging computer code implementing a calculation. Note that when 
using the spinor-helicity techniques, the use of non-trivial (Gram determinant) identities 
may be required to show cut consistency. 

Each sewn and integrated cut is gauge invariant independently when all external legs 
are on shell. Its absorptive part is, after all, equal to the phase-space integral of a gauge- 
invariant non- forward matrix element. The dispersion integral of such a quantity is gauge 
invariant as well. At the integrand level, this is reflected in the disappearance of any 
dependence on the light-cone vector n introduced by the sum over gluon polarizations. 
(One typically needs to make use of momentum conservation to see this explicitly.) The 
same statement is not true if we are computing an object with off-shell legs, as is the case 
for the splitting amplitude. In that case, some light-cone denominators will survive. As 
we will discuss in a later subsection, others will cancel algebraically, and it is possible to 
predict in advance which must disappear. 

The procedure we have described above is not the only way to merge information from 
different cuts to yield a single function with the correct cuts in all channels. One can also 
imagine performing the merging step after integration. That is, one could integrate the 
sewn integrand in each channel separately, and then search for a function whose branch- 
cut discontinuities match those found in each of the separate invariants. This requires 
the use of a nonredundant basis of master integrals. Operationally, one would rewrite the 
result of integration in each channel in terms of this basis, throwing away integrals with 
no discontinuity in the given channel. One could verify cut consistency by checking that 
the coefficients of a given master integral with discontinuities in multiple channels are in 
fact the same in the different computations in these channels. In combining channels, one 
would then take the result from any of the channels (that is, pick one, rather than adding 
together the different contributions). This is in fact the procedure that was presented in 
ref. |]ll[. It may be advantageous for some calculations, though in general it will require 
the computation of superfluous integrals. As with the basic sewing procedure, there are 
adjustments that would be required in a two-loop computation, which we will discuss in 
the next subsection. 

4.3 Sewing at Two Loops 

As for the one-loop case, we start our two-loop discussion with the example of a four- 
point amplitude. We again consider the contribution with a given ordering of the external 
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Figure 18: Color-ordered Feynman diagrams for the two-loop four-point amplitude in Tr (f)^ field 
theory: (a) planar box diagrams (b) triangle-in-box diagrams (c) bubblc-in-box diagrams (d) s- 
channcl ladder triangle diagrams (e) s-channel bubble-in-bubblcs. (f) s-charmcl bubblc-in-trianglc 
diagrams (g) s-channel triangle-pair diagram (h) s-channel lizard-eye bubble (i) s-channel double 
bubble (j) s-channel triangle-bubble diagram (k) s-channel triangle-in-triangles. The corresponding 
f-channel diagrams for (d-k) are not shown explicitly. 



legs, and restrict attention to planar diagrams. In the one-loop example, we saw that we 
must pay attention to potential double-counting in assembling contributions from different 
channels. At two (and higher) loops, we must also confront potential double-counting in 
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Figure 19: The s-channel two-particle cuts of the planar double box. 

contributions to a given channel. To understand how this arises, start once more with the 
Feynman diagrams for this planar ordered amplitude, shown in figure 

Consider in particular the planar double box shown in figure p!8|(a). It has both two- 
and three-particle cuts. The two-particle cuts, shown in figure contain a product of a 
one-loop amplitude and a tree amplitude, with four external legs apiece. There are two 
separate cuts. 

The three-particle cut in the planar double box is a sum of two terms, shown in 
figure 1^. It corresponds to a product of two five-point tree amplitudes. The complication 
here arises from the fact that a given term in the two-loop amplitude may contribute to 
both terms in the two-particle cut; both terms in the three-particle cut; or to both two- and 
three-particle cuts. In such a case, when we reconstruct the original integrand, we must 
count it only once. A simple analogy would be an integrand of the form X^; since cutting 
is analogous to differentiation, we would have 

[^']|cut = 2^lcut^> (4.17) 

and promoting the latter back to an integrand requires a factor of 1/2 just as it would 
for integration. Note that if we denote the cut momenta by ^1,2, 3) then the terms which 
contribute to both two- and three-particle cuts necessarily contain the propagators 

^ (4.18) 



2 



using the labeling of legs in the three-particle cut. 

We will focus here on the leading-color (planar) contributions to the two-loop am- 
plitude, though of course an analogous procedure applies to the subleading contributions 
(which include non-planar integrals). Let us begin the sewing procedure with the two- 
particle cuts in the s channel. We must take the one-loop amplitude (|4.7D , with legs (3,4) 
replaced by legs {£2 , ^1 ) , and £ replaced by p in the integral, and sew it to the tree amplitude 
^(0) (-^1,-^2, 3, 4), obtaining 

. r d^p d^£ 
V (2^(2^ 
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Figure 20: The s-channel three-particle cuts of the planar double box. 




Figure 21: The s-channel 'double' two-particle cuts of the planar double box. 
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si2 {p + kiYp'^{p -I- A;i -I- {p + kiYp'^{p - lY{p + ki + k2Y 

11 1 1 



' ^kxYfp^{p-t-kxY {l^kxYp\p-kxY{p- i-kxf 
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1 



[i + kxY {p + ^2) V(p - ^ - ^1)^ 



1 



1 



+ 



1 



{i + ki + k2YP\su {i-k^y 



where we have relabeled ii £. Restricting attention to those terms with explicit powers 
of 1/5^2 or l/sf2 5 we have 

d^p d^e 



1 



(27r)« (27r)^ 52^(^ + ^:1+^2)2^2 



+ 



1 



S12 P^{p + ki + A;2)2 p^{p — tYip + ki + k2Y 
1 



+ 



+ 



{p + kiYp^ip + ki + k2Y p^{p + ki + k2Y{l - k^Y 



+ 



(4.20) 



These terms are sufficient to illustrate the issues associated with potential double-counting. 

There is a similar contribution from sewing ^'^'^^(1, 2, £2, ^1) to ^1, — ^2, 3, 4). If 

we now take the cut of the expression ( 4.ig| ), however, we discover that we can cut not only 
the i loop, but also the p loop. The corresponding terms would appear not only in this 
sewing, but also in the other contribution (with the loop amplitude on the right-hand side 
of the cut). This would be a source of double-counting. To correct for it, we must subtract 
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Figure 22: The s-channel three-particle cuts of the hzard-eye bubble. 



from the sum of the two contributions those terms that can be cut both ways. These are 



the terms extracted by a 'double cut' (shown in figure 21), again in the sense of requiring 
that the propagators giving rise to a two-particle cut be present in both the i and p loops, 



(- + -)(---)(---;' 



(4.21) 



or promoted back to a two-loop integral. 



{2tt)D (27r)^ s\^{i + kx + fe2)2^V(p + ki + 



1 



+ 



1 



+ 



+ 



+ 



(4.22) 



which is identical to the terms in eq. ( [4.20| ), obtained from the first of the two-particle cut 
contribution. Counting the contribution only once then gives eq. (1^ as the result of 
combining the two two-particle cut contributions. 

Next, we must consider the three-particle cuts. We begin with the product of two 
five-point amplitudes. 



1 



+ 



1 



+ 



+ 



+ 



1 



S^2^3'5£il ■S^i£2*12 ■S£i£2'S^32 Si^l^Su S£j^lSi^2 



(4.23) 



and ^W(-^i,-^2,-^3,3,4). 

Let us focus on the terms containing a factor of l/s^2- These are 



■I 



1 



{27r)D {27r)D sf^ ijelili + £3 + ki + k2f 

1 2 

+ 



(4.24) 



+ 



+ ki + k2YY {h + ki + k2Y{£i + ki + k2Y [{£i + ki + k2YY)' 



The first and last of these terms correspond to the diagrams of figure 11(e). The middle 
term corresponds to the diagram of figure llB(h); but it appears in the product with a 
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Figure 23: The i-channel cut of the planar double box. 



(superfluous) factor of 2. This factor is due precisely to the fact that this term can be 



cut in two different ways, corresponding to the two terms depicted in figure |2^. We must 
remove this double counting, by subtracting those terms which contribute twice; this leaves 
us with 



■I 



d^li d^is 1 



1 



{27r)D {27r)D sj^ + ^3 + A:i + fea)^ 

1 1 

+ 



(4.25) 



+ fcl + ^5 



\212 



+ 



(^3 + A;i + A:2)2(£i + A;i + fc2)2 {{h^kx^k^YY 



Finally, we must combine the two- and three-particle cuts. Again, we can add the two 
contributions, and remove terms which appear in both, for example by removing terms in 
the three-particle cuts which have (any) two-particle cut. Of the terms listed explicitly in 
eq. ( f4.25| ), only the middle term has a two-particle cut. Removing it (thereby performing 
the required merging), and relabeling i\ I, ^ —p — ki — k2, we obtain. 



d^p d^£ 1 



1 



(27r)^(27r)« sj^ P{p + ki + k2y 

1 1 1 



+ 



Si2P^{£+ki+k2f p^{p-lf{i + ki + k2f 
1 1 



+ 



+ 



+ 



{p+kifp^{l+ki+k2Y p^{l-kif{i + ki+k2Y 

1 1 



+ 



(4.26) 



+ 



+ ••• 



again in exact agreement with the terms that would emerge from a Feynman-diagram 
calculation. (The reader may wonder why the first and last terms in eq. (I^D, which 
correspond to the diagrams in figure 18(e), do not contain two-particle cuts. Naively, 
these diagrams do contain two-particle cuts; but a closer inspection shows that such a cut 
would have a bubble on an external line; such bubbles are scale-free and hence vanish in 
dimensional regularization.) 

Note that not all contributions will require these subtractions to match overlaps; for 
example, the t-channel cut of the planar double box shown in figure |2^ has only one 
contribution, a three-particle cut. 

In practice, it appears to be more efficient to start with the two-particle cuts, and then 
merge in additional terms from the three-particle cuts. In addition to the cut projection 
operator Vs, which yields the part of its argument that has any cut (two- or three-particle) 
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in the specified channel s, we will also make use of four additional projection operators. 

(2) (3) 

These are Vs , which extracts terms containing a two-particle cut; Vs , those containing 

(2) 

sl three-particle cut; V u , those containing a 'double cut', that is whose two-particle cut 



can be cut again (for example, figure pID ; and "P^x"*) those containing a contribution to both 
terms in the three-particle cut (corresponding to terms which have all the propagators cut 
in figure 21 and the center propagator in addition). Note that these projection operators 

(2) 

do not necessarily remove terms with other cuts; for example, Vs may yield an expression 
containing three-particle cuts in addition to the two-particle cut. In cases where there is 
only a single two-particle cut, or where there is only a single contribution to the three- 
particle cuts, V^'^^ and V^^^ are understood to vanish. 

Using these projection operators, we can define a complete promotion operator for 

JO 

a given channel, 

[ A( )({A}, a„ . . . , bj)A(-^{bj + 1, . . . , a,- - 1, {-£i}) = 



lp{2) 

2 -II 



A'^'\£,aj,...,bj, 



X ^(i)(^ + Kj,bj + l,..., aj - 1, -£) + 



(4.27) 



+ A«(^,, 



A^^^£ + Kj,bj + l,...,aj-l,-£) 



+ 



1 - -rS 

2 



A^'>\£i,£2,a,,...,bj,-£i-£2 
xA^^\£i+£2 + Kj,bj + l,. 



,aj - l,-£2,-£i) 



where A^^^ denotes the integrand for the one-loop amplitude. Here we consider only the 
planar case, so that all cut-crossing momenta are color- adjacent; but the construction 
generalizes in a straightforward way to non-planar amplitudes. As in the example discussed 

(2) (S) 

earlier, the role of the factors of 1/2 in front of and V).J is to remove double-counting 
that occurs in cutting, ensuring that each term contributes to the integrand with the correct 
coefficient. 

Using this complete promotion operator, the complete integrand for the two-loop am- 
plitude A(-^\l,... , n) is given by an expression very similar in form to that for the one-loop 
amplitude 



/(2) 



%\^A(-\{£,},a{s,),...,b{sj))x 
X A^-\b{sj) + l,...,a{sj) -!,{-£,}), 



(4.28) 



but with the complete promoted cut being used where the basic one was in the one-loop 
case. 

A practical algorithm is similar in structure to that at one loop, but requires additional 
merging over two- and three-particle cuts. The first three steps (construct ordered set of all 
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(a) (b) 

Figure 24: The two contributions wherein one sews a one-loop amphtude to a tree-level amplitude 
to produce a contribution to the two-loop integrand. 



channels, determine labelings, initialize integrand to zero) are the same. As in the one-loop 
case, from a formal point of view, the relabelings below are not required by the method, but 
since we will ultimately be feeding the resulting integrands to an integration machinery, 
we might as well incorporate the standardization at an early stage in the calculation. 

In the two-loop case, 'irreducible numerators' appear. These are dot products of loop 
momenta and external momenta that cannot be written as linear combinations of inverse 
propagators. Such terms must be added to the basis of expressions described just before 
the one-loop algorithm. With this modification, the algorithm again continues after the 
first four set-up steps by iterating over all channels. For each channel {ka^, . . . ,ki,.), 

5. Sew first two-particle cuts. Form the product of the integrand for the on-shell one-loop 
amplitude ^^^-'(^i, , . . . , 6j, £2) and the tree amplitude A^^\—i2, bj + 1, . . . ,n,l, . . . , 
aj — (where ii + £2 + Kj = 0), summing over the different particle types 

and helicities that can circulate in the loop. The sum must in general be performed 
in = 4 — 2e dimensions. Use the polarization- vector identities ( 4.12| ) to express 



everything in terms of the cut momenta ii,i2 and the light-cone vector n. (As in the 
one-loop case, this will introduce physical projectors into the integrand.) Multiply by 
—1 for each fermion loop is created by sewing. Put in the propagators crossing the 
cut, ^ for each cut leg. Note that in general four-dimensional Fierz identities may 
not be used. Fermion traces over internal fermion lines should be performed first, 
re-expressing the result in dot products, with summation over gluon polarizations 
performed afterwards. Rewrite terms using the basis expressions, and expand sums so 
that each term can be classified as the daughter of a single parent diagram. Complete 
dot products in the denominator to form propagator denominators adjacent to the 
cut, 2£ ■ ki — > ±{i lb ki)'^, the -|- corresponding to ki on the left side of the cut, 
the — sign to those on the right side of the cut. The sewing operation is depicted 
schematically in figure |2^(a), and corresponds to the promotion of the second term 



in brackets in eq. (|4.27D 



6. Put into canonical form. Use momentum conservation to reduce the number of cut- 
crossing momenta (now loop momenta) appearing in spinor traces, and then expand 
spinor traces and all dot products in terms of the basis. (As explained above, this 
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expansion may not always be necessary.) When using the spinor-helicity method for 
external polarization vectors, or explicit helicity states, we will obtain spinor strings 
of the form (jf • As explained earlier, one can complete these to a 

trace, then convert the trace to dot products and Levi-Civita tensors. The latter can 
be converted to Gram determinants (and thence to dot products) by multiplying (and 
dividing) by a Levi-Civita tensor involving only external momenta. (Spinor strings 
involving only external momenta need not be manipulated, obviously.) 

7. Relabel. For each term obtained in step ^ relabel the momenta to the standard 
labeling for its parent integral. Where required, insert factors of squared momenta 
in the numerator and denominator to match a 'parent' diagram. In some cases, it 
will be possible to obtain different parents by inserting different factors; it doesn't 
matter which one is picked. 

8. Sew second two-particle cuts. Form the product of the on-shell tree amplitude 
A^^\ii,aj, . . . , 6j, £2) and the integrand for the on-shell one-loop amplitude ^(^)(— ^2, 
6j + 1, . . . , n, 1, . . . , Oj — 1, —(-i) (where £1 -|- £2 + Kj = 0), again summing over the 
different particle types and helicities that can circulate in the loop. As in step |5|, 
put in fermi minus signs and the cut propagators, complete dot products to propa- 
gators, and rewrite terms using basis expressions. The sewing operation is depicted 
schematically in figure 24(b), and corresponds to the promotion of the first term in 
brackets in eq. ( |4.27 ). 



9. Put into canonical form. Use momentum conservation, the cut condition, and ex- 
pansion of spinor traces, along the lines of step ^. 

10. Relabel. For each term obtained in step ^ relabel the momenta to the standard 
labeling for its parent integral, inserting factors of squared momenta in the numerator 
and denominator where required. 

11. Merge two-particle cuts, removing double- counting. To the result obtained in step |^, 
add those terms obtained in step |l^ not present in the result from step |7[ (Equiv- 
alently, add the results of steps ^ and |l^ and then subtract terms present in both 
expressions. These latter terms are those with two pairs of propagators present, i.e. 
those cut through in the 'double' two-particle cut.) This step corresponds to the 
operator (1 — ^v'^J)) in eq. ( 4.27| ), and yields the sewn two-particle cut in the given 



channel. 

12. Sew three-particle cuts. Form the product of on-shell tree amplitudes ^^''-'(^25 ^i, Oj; 
. . . , 6j,£3) and A(°)(-^3, 6j + l, . . . ,n, 1, . . . , aj -1, -h, -£2) (where £i+e2+£3 + Kj = 
0) , summing over the different particle types and helicities for the ii . As in step |5[ put 
in fermi minus signs and the cut propagators, complete dot products to propagators, 
and rewrite terms using basis expressions. This operation is depicted schematically in 
figure 25, and corresponds to the promotion operation in the last term in eq. ( 4.27] ). 
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Figure 25: Sewing together two tree-level amplitudes to produce the three-particle cut contribution 
to the two-loop integrand. 



13. Put into canonical form. Use momentum conservation, the cut condition, and ex- 
pansion of spinor traces, along the lines of step E. 



14. Relabel. For each term obtained in step 13, relabel the momenta to the standard 



labeling for its parent integral, inserting factors of squared momenta in the numerator 
and denominator where required. 

15. Remove double- counting. Multiply by one-half those terms (if any) which contribute 
to both terms in the three-particle cuts. These are the terms that contain all five 
propagators l/ij, l/i^, 1/1% l/(£i + £2?, and 1/(^2 + h? ■ As before, in gen- 
eral, identification of these terms may require expansions of numerators, and use of 
momentum-conservation identities. This step corresponds to the operator (1 — ^V^^^) 
in eq. ( ^.27] ), and yields the sewn three-particle cut in the given channel. 

16. Merge two- and three-particle cuts. To the two-particle cut obtained in step |ll], add 
those terms in the three-particle cut not present in the two-particle cut. (Equivalently, 

add the two expressions, then subtract those terms in the three-particle cut also 

(2) 

present in the two-particle cut, corresponding to the action of the operator Vl in 



eq. ( |4.27] ).) This yields the complete sewn expression (4.27) in the current channel 



17. Clean. Remove all terms which have no cut in the current channel. (Such terms might 
have been introduced by earlier algebraic manipulations.) Using a canonical basis as 
described above, remove any terms which do not contain both cut propagators. This 
step corresponds to the operator Vsj in eq. (4.28). 



18. Merge. Remove all terms in the current-channel sewn expression that already appear 
in the net integrand. That is, remove any term which has cuts in a previously- 
processed channel. Using a canonical basis as described above, it suffices to pick out 
and remove terms that have a pair or triplet of propagators corresponding to a cut in 
a previously-processed channel. This step corresponds to the operator Jlti ( ^ 



in eq. (4.2. 



19. Accumulate. Add the remaining terms in the current-channel sewn expression to the 
integrand. This step corresponds to the sum in eq. ( [4.28| ). 
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20. Continue with the next channel at step 



As in the one-loop case, one could alternatively do the merging after integration rather 
than before. One again needs a nonredundant basis of master integrals; but here, one 
needs to adjust the coefficients of some integrals to account for the double-counting issues 
discussed in the earlier example, and handled in steps and |l^ of the two-loop algorithm 
above. The master integrals must be chosen so that each is associated with a definite 
overall correction for double-counting; and one must keep track of the original set of cut 
lines in each term, alongside the integral result. In performing integral reductions, one 
must eliminate integrals in which cut propagators are cancelled. In the computation of a 
given cut in an amplitude, one first needs to merge the different two-particle cuts, as there 
is now more than one contribution. Those master integrals with two contributions to the 
three-particle cuts in the given channel would have their coefficients decreased by a factor 
of two. The three-particle cuts must then be merged with the two-particle cuts to obtain 
the full set of terms for the given channel. One can check cut consistency between different 
channels just as at one loop. Master integrals with discontinuities in two or more channels 
must appear with the same coefficient in each channel. The merging of contributions from 
different channels also proceeds in the same manner as at one loop: for those master 
integrals with discontinuities in multiple channels, one would take the result from any of 
the channels (that is, pick one, rather than adding together the different contributions). 

4.4 Merging with Legs Off Shell 

When computing splitting amplitudes, we have an off-shell leg in the problem, and not 
all light-cone denominators cancel from the final integrand as they do for fully on-shell 
scattering amplitudes. For two-loop splitting amplitudes, one side of the cut will have 
an on-shell tree or one-loop amplitude integrand, while the other side contains a tree- 
level splitting amplitude or one-loop splitting amplitude integrand. The three-particle cuts 
require the use of 1 — > 3 tree- level splitting amplitudes, which should be calculated in 
light-cone gauge. This introduces additional light-cone gauge denominators into certain 



integrals. As at one loop, sewing gluons across the cut using eq. ( 4.12 ) will introduce 
similar denominators. These denominator factors will survive into the integration, and in 
fact play a crucial role in obtaining the right 'complexity' of integrals. However, as we 
shall discuss in section ^ they will only arise in lines corresponding to cut momenta or 
connected directly to the off-shell vertex. 

The presence of non-canceling light-cone denominators does not alter the merging 
procedure described above: these denominators simply go along for the ride. Following the 
discussion in section ^ we may expect all surviving light-cone denominators to produce 
integrals properly regulated by covariant dimensional regularization, and this is indeed 
the case in the calculation we have performed. In contrast, a light-cone gauge Feynman 
diagram approach would contain ill-defined diagrams requiring an additional prescription 
such as the PV or ML prescriptions. The ill-defined contributions may cancel in the sum 
over all diagrams before integration, if a great deal of care were taken to align momenta 
properly across different diagrams, as guided by the unitarity cuts. 
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(a) (b) (c) 

Figure 26: The three-particle cut (a) contains a Ught-cone denominator on the leg indicated by 
an arrow, but the two-particle cut (c) will not contain one on that leg. Cut consistency demands 
that the light-cone denominator indicated in (a) must cancel from all terms where the two-particle 
cut (c) does not vanish. The light-cone denominator arising in the three-particle cut indicated in 
(b) is compatible with the ones in the two-particle cut (c) and therefore does not need to cancel. 



As discussed in section 4.2, if all legs are on shell then ah light-cone denominators 
cancel from the integrands of each cut. Even if some legs are off shell, some light-cone 
denominators (introduced by the physical-state projector a cut) can be canceled prior to 
integration by combining information from different cuts. The possibility of these cancella- 
tions is dictated by cut consistency. If a given term in an integrand has a cut in more than 
one channel then it must appear with the same coefficient in each such cut. A light-cone 
denominator absent in a term in any one of the cuts must also cancel in all other cuts. 
This cancellation may, however, not be manifest. To cancel the denominators, in general, 
momentum conservation rearrangements are required. We note that in performing the 
calculation it is generally helpful to explicitly cancel as many light cone-denominators as 
possible, to reduce the number and complexity of integral types that need to be evaluated. 
(This goes beyond the automatic absence of dangerous light-cone denominators discussed 
at the end of section ^.) 

As an example, consider the contributions to the two-loop splitting amplitude depicted 



in figure 26, with both two- and three-particle cuts. In the cuts a light-cone denominator 
will appear on each sewn gluon line. The three-particle cut contains contributions with the 
light-cone denominators indicated in figure p6|(a) and (b). However, if these contributions 
also have the two-particle cut indicated in figure |26|(c), then the light-cone denominator 
indicated in (a) must cancel algebraically in the three-particle cut after suitably combining 
terms using momentum conservation. The reason is cut consistency: in terms which have 
both two- and three-particle cuts, we must obtain identical results. In the two-particle 
cut (c) everything to the left of the cut is gauge invariant, and cannot contain light cone- 



denominators. (If we consider, instead of figure 26, a contribution for which either of the 
two cut propagators in (c) are absent, then we cannot determine from that two-particle 
cut whether the denominator (a) should cancel or not. However, by considering the other 
contribution to the three-particle cut we find that the absence of the top cut propagator 
in (c) allows the denominator in (a) to be present; but in the absence of the bottom cut 
propagator alone, the denominator in (a) has to cancel.) More examples of arguments of 
this type can be found in the discussion of the types of integral topologies encountered in 
figure ^ in section |^. 
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(a) (b) (c) 

Figure 27: The three contributmg cuts. 

5. Generation of Splitting Amplitude Integrands 

Using the unitarity-based method described in the previous section we have generated 
the integrand for the two-loop sphtting ampHtude. Here, instead of sewing together two on- 
shell amphtudes, we sew an on-sheU amphtude to a sphtting amphtude. The tree sphtting 
amphtudes are given by a current with one leg off shell Q, or equivalently, an amplitude 
with one leg off shell divided by the squared momentum of that leg. In general, the current 
should be evaluated in light-cone gauge. It can be computed either recursively ||7|, or 
via color-ordered Feynman rules. In a computation of the one-loop splitting amplitude, 
only the 1 — > 2 tree splitting amplitude enters; in the computation of the two-loop one, 
both the 1 — > 2 and 1^3 splitting amplitudes enter. The 1^2 tree-level splitting 
amplitude is basically just a vertex, and so has the same expression in both covariant and 
light-cone gauges; but the expression for the 1 — > 3 amplitude is not the same in different 
gauges [P, m. 

The 1 — > 3 splitting amplitude governs the universal behavior of tree-level gauge-theory 
amplitudes in triply-collinear limits, where three color-adjacent momenta ki^2,3 become 
collinear, and all invariants ii23, si2) S23 are comparably small. These limits have been 
described previously in refs. |^, 

There are three basic types of contributions we need to consider: 

(a) two-particle cuts with an on-shell one-loop four-point amplitude integrand sewn onto 



the tree- level 1 — > 2 splitting amplitude, depicted in figure |27|(a); 



(b) two-particle cuts with an on-shell tree-level four-point amplitude sewn onto the one- 
loop 1^2 splitting amplitude integrand, depicted in figure p7|(b); 

(c) three-particle cuts with an on-shell tree-level five-point amplitude sewn onto the tree- 
level 1 — > 3 splitting amplitude, depicted in figure |27|(c). 



In our calculation, there is a cut in only one channel (the invariant S12 in figure 27), 
so cross-channel projection and consistency issues do not arise. However, the three con- 
tributing cuts do need to be constructed and merged as described in the previous section, 
as there are terms that are common to different contributions. As in the one-loop case, 
the light-cone denominators inserted by the physical projection operators on the cut are 
crucial to getting the correct answer. This is also true for the light-cone denominator on 
a line connected to the off-shell vertex, contained in the 1 — > 3 splitting amplitude. Its 
absence would result in different and inconsistent integrands emerging from the two- and 
three-particle cuts. 
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Figure 28: Examples of non-planar contributions to the two-loop g ^ gg splitting amplitude, all 
of which have vanishing color factors. 

In principle, non-planar topologies could enter into the g ^ gg splitting amplitude, for 
example the crossed triangle graphs shown in figure |2^. However, it turns out that all such 
non-planar graphs, for both pure-glue and fermion-loop contributions, have vanishing color 
factors: one simply dresses the diagrams with their color factors and performs the color 
algebra to demonstrate this. Accordingly, we only need to sew amplitudes into planar 
configurations for our calculation. As we shall show in section ^, the only color factor 
that arises in the pure-glue contributions is the leading-color one, namely C\ = N^, for 
the simple reason that there are no other color Casimirs at this order. The fermion- 
loop contributions can be divided into leading-color {CaNj = N^Nj) and subleading-color 
{{Ca — 2CF)Nf = Nf/Nc) terms. Both types arise from planar diagrams. Each planar 
diagram contributing to the subleading-color terms can be drawn with the virtual gluon on 
the inside of the fermion loop, whereas each for the leading-color terms can be drawn with 
the virtual gluon on the outside. (Some diagrams with bubble insertions can be drawn both 
ways, and are proportional to Cp = {N^ — l)/{2Nc).) We discuss the full color dressing of 
the splitting amplitudes in section ^ 

For splitting amplitudes with external quarks, g ^ qq and q ^ qg, the non-planar 
color factors no longer vanish. Non-planar two-loop three-point integrals with light-cone 
denominators will be required. However, these integrals should be amenable to the same 
methods used in the present computation. 

In carrying out the calculation, we group terms in each sewn cut into different inte- 
gral topologies, according to their propagators and light-cone denominators. To reduce 
the number of independent topologies, we perform a partial-fraction decomposition when 
certain multiple light-cone-denominators appear. For example, if pi +p2 = ki + k2 for two 
internal momenta pi, p2, we substitute 

1 1 
pi ■ n p2 ■ n {ki + k2) ■ n 

Although the physical-state projectors for a two-particle cut generate two different light- 
cone denominators, the use of eq. ( |5.1| ) allows us to consider only one at a time. For 
triplets of light-cone denominators associated with a three-particle cut, only two have to be 
considered at a time. We also use a symmetry of the integrals under ki <^ k2 {z ^ 1 — z) 
to restrict the number of integral topologies to those described in figure ^ in the next 
section. 



+ 

pi ■ n p2 ■ n 



(5.1) 
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After relabeling momenta circulating in the sewn cut (see appendix to match the 
labeling used by the integration routine for a given topology (see the next section), we 
obtain an integrand of the form 



(5.2) 



ignoring factors that come out of the integral such as £i ■ ej. Here p and q are the loop 
momenta; njP? shorthand for the set of both Feynman propagators and light-cone 
denominators for the topology. The polarization vectors for particles i = 1,2, P are denoted 
by £i, and are taken to satisfy the light-cone-gauge condition £i ■ n = 0, as well as the 
transverse condition £i ■ ki = 0. (Note that although kp is slightly off shell, any terms 
arising from ep ■ kp ^ will not be sufficiently singular in the fcp — > limit to contribute 
to the splitting amplitude, and hence we may as well set ep ■ kp to zero.) 

In the dot products of polarization vectors with loop momenta, £i ■ p and • q, the 
D-dimensional loop momenta are effectively projected into four dimensions, because the 
external physical polarizations are four dimensional. Following the discussion of ref. we 
write the loop momenta in these dot products as linear combinations of four independent 
momenta: ki,k2,n, and the dual vector 



Then we have, for example, 

£i ■ p = c^£i ■ ki + c^Ei ■ k2 + c^Ei ■ n + c^Ei ■ V , 



(5.3) 



(5.4) 



where 



4 



4 



2ki 


k2ki ■ n 




1 


2ki 


k2k2 -n 




1 


2k2 


nki ■ n 




1 



2ki ■ k2ki ■ nk2 ■ n 



— k2 ■ np ■ ki + ki ■ np ■ k2 + ki ■ k2P ■ n 
k2 ■ np ■ ki — ki ■ np ■ k2 + ki ■ k2P ■ n 
k2 • np ■ ki + ki • np ■ k2 — ki ■ k2P ■ n 



1 



2ki ■ k2ki ■ nk2 ■ n 



p ■ V . (5.5) 



After this substitution we use momentum conservation to express each integrand in terms 
of a basis consisting of: 

1. inverse propagators; 

2. light-cone denominators and numerators; 

3. irreducible numerators (dot products of loop momenta with ki and k2 which cannot 
be written in terms of inverse propagators); and 



4. dot products of loop momenta with the dual vector ( |5.3D , coming from and c?. 
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Following the merging procedure described in the previous section, at the integrand level, 
we obtained a single expression with correct cuts in all channels. Using the formulae in 
section 4 of ref. we then replaced dot products of loop momenta with the dual vector, 
p ■ V and q ■ V, by other elements of the basis. (These formulae are valid at the level of 
integrals, not integrands, so they should only be applied after checking cut consistency.) 
Typical propagator momenta, light-cone dot products, and irreducible numerators that 
appear in the calculation are given in eqs. ( |6.4| ) and ( |6.6| ). This produces an expression 
ready to be integrated. We describe the integration method in the next section. 



6. Integrals 
6.1 Introduction 

All the integrals encountered in our computation of the two-loop splitting amplitudes are 
3-point integrals with one external massive leg, and two massless legs. The massless legs 
carry momenta ki and k2. The massive leg carries momentum kp = ki + k2 with kp = 
{ki + k2)^ = 2ki ■ k2 = s. We consider the time-like case, s > 0. For the space-like case, 
we take the splitting to be {—ki) {—kp) + k2, with fcp = s < 0. Here {—kp) and k2 
carry longitudinal momentum fractions x and 1 — x respectively, with x = 1/z. The space- 
like case may be obtained from the time-like case using analyticity. As there is no other 
dimensionful parameter in the problem, the dependence on s is determined by dimensional 
analysis to be oc (— s)"^"^. To reach the physical range < x < 1, it is also necessary to 



continue the momentum fraction z to values larger than 1, as we shall discuss in section [7.4. 

In addition to standard propagator factors of the form in the denominator, and 
tensors in the numerator, there can also be denominator factors from light-cone gauge (or 
physical state projection, of the form l/{p ■ n) where is the light-cone gauge vector, 
n? = 0. The vector n is also used to define the collinear momentum fraction z, according 
to the kinematic relations ( |3.8| ). We also rescale n so that it obeys eq. ( |3.1[1| ), kp ■ n = —s. 

Two-loop 3-point integrals containing light-cone denominator factors have not been 
encountered previously, and contain non-trivial dependence on z. For the case of a gluon 
splitting to two gluons, we shall see that the non-planar topologies all have vanishing color 
factor, and hence the corresponding integrals do not need to be calculated. The integrals 
from the planar topologies can be reduced to a set of 13 master integrals, using identities 



based on integration by parts (IBP) |£0| and Lorentz invariance |20]. The master integrals 



can be evaluated as Laurent expansions in e with the aid of differential equations in z, 
along the lines of refs. |18, pC| ]. 

The Laurent expansions of each master integral, through the order in e required to 
obtain the splitting amplitudes to order e*^, can be expressed in terms of logarithms, plus 
the polylogarithmic functions defined [86| by, 

Un{x) = Y.1^= -^^n-l{t), (6.1) 

-•—1 Jo 

/ -Ml-t). (6.2) 
Jo ^ 



i=l 

rx 

U2{x) 
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Here we need Li„(x), for n = 2,3,4, and the argument x can be z, 1 — z or —(1 — z)/z. 
(For n < 4, identities relate some of these polylogarithms to each other.) For the order 
terms in the sphtting amphtudes, the corresponding set of functions (with n extending up 
to 5) is not sufficient; instead harmonic polylogarithms |87, pO| are required. 



In section 6.2 we describe the procedure for reducing the two-loop three-point integrals 



with light-cone denominators to master integrals, and give a list of master integrals required 



for the g ^ gg splitting amplitude. In section 6.3 we illustrate how to derive differential 
equations for the master integrals, and present the differential equations. We also give the 
Laurent expansions of the master integrals. 

6.2 Reduction Procedure 

A typical integral encountered is the two-loop nested double triangle integral with two 
light-cone denominators inserted, 

L{ui,U2,us,...,ug)= I -^-^ [I , (6.3) 

1=1 ^ ' ' 

where D = 4 — 2e, 

Pi = q, P2 = q + ki+k2, P3=P, Pi=p-ki-k2, 

P5 = p-ki, P6=p + q, 

P7 = Q-n, pl = 2q- ki, pj = {p - ki - ■ n, (6.4) 

and the are integers. The momentum routings for this integral are depicted in fig- 
ure p9|(a). The off-shell external momentum ki + k2 flows in from the right of the diagram, 
and splits into on-shell momenta ki and k2 flowing out to the left. The uncircled numbers 
adorning the internal lines label the internal momenta pi, whose squares are the 'ordinary' 
denominators pf,i = l,2,...,6, appearing in eq. (|6.3|) . The circled numbers correspond to 
the light-cone denominators, in this case pj and pg, which are linear in n. Associated with 
each circled number is an arrow on an internal line. The arrow is a reminder of the direc- 
tion of the internal momentum which is Lorentz-contracted with n to form the light-cone 
denominator. One of the in eq. ( |6.3| ) is not shown in figure |2^(a): pg = 2q ■ ki. In this 
topology, it only appears in the numerator, i.e. fg < 0. It will appear as an irreducible nu- 
merator in the tensor integrals arising from the numerator algebra generated in evaluating 
the cuts. Its presence is also required to close the IBP equations. 

We also encounter the 'pentabox' integral topology (but with one ordinary prop- 
agator cancelled), again with two light-cone denominator insertions, 

f d^p d^q -y-T 1 
LPB{ui,U2, ^^3, . . . , ^9) = / n ' (6.5) 



i- 



where 



pi = q, P2 = q + ki + k2, P3=p, P4 = q + ki, 

P5=p-ki, P6=P + q, 

P7 = {q + ki + k2) ■ n, pI= p - n, p^ = 2p ■ k2. (6.6) 
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(a) (b) 

Figure 29: (a) The planar double triangle integral defined in eq. ( |6.3[ ). (b) The planar double 
triangle integral LPB{vi) defined in eq. ( |6.5| ). In each case fci and k2 are the two outgoing massless 
external lines. The internal lines, or 'ordinary' propagators, labelled with an uncircled integer i, 
i G {1, 2, . . . , 6}, carry momenta p^. The special light-cone denominator factors pi ■ n which can be 
present are marked with an arrow on an internal line, and the corresponding integer label is circled. 
The arrow serves to remind one of the direction of the loop-momentum used in their definition. 
Propagators 8 (in case (a)) and 9 (in case (b)) do not appear in denominators in these integrals. 
The expressions for pg and pg are given explicitly. 



The momentum routings for this integral are depicted in figure 29(b). In this case pj and 
Pq are the circled, light-cone denominators, and pg only appears in the numerator, < 0. 

Besides the two types of integrals depicted in figure there are several more types, 
which differ from either L or LPB only in the location of the light-cone denominators. All 
the light-cone configurations are shown in figure 30. The leftmost two cases, L and LPB, 
are equivalent to figure ^9|(a) and figure p9|(b). The momentum routings and the labels for 
the external legs and for propagators 1 to 6 are exactly the same as in figure ^9], so we have 
suppressed them in figure^. The remaining cases define the integrals -D(t'i), F{ui), G{i'i), 
H{i>i), J{vi) and M{yi), by analogy to eqs. ( |6^ ) and (^). As in the L and LPB cases, 
one extra propagator, linear in the loop momentum but not a light-cone denominator, is 
required for tensor integral reduction and for the IBP equations to close; the expression 
for this propagator is shown explicitly in the figure. 



The symbol 'x' on a line in figure 30 indicates that that ordinary propagator never 
appears in the denominator, when the light-cone denominators indicated by the arrows are 
present. As discussed in section 4.4, the unitarity cuts guarantee these facts. For example, 
suppose the propagator marked with a 'x' in the J topology were present. Then such 
an integral has a 2-particle cut to the right of the circled propagator 8. (If propagator 
2, which also must be present in the 2-particle cut, did not appear in the denominator, 
then the integral would become a massless external leg integral, which vanishes trivially in 
dimensional regularization.) But we know that no light-cone projectors should appear to 
the left of such a cut, so either the circled propagator 8 or the x -marked propagator must 
be absent in J. If the circled propagator 8 is absent from the denominator, we consider 
the integral to belong to the F topology instead. 



Note that we have not explicitly shown integrals that are related to the ones in figure 30 
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L D G H 




LPB F J M 

Figure 30: List of required light-cone denominator insertions for the g ^ gg splitting amplitude. 
The leftmost two cases are equivalent to figure ^(a) and figure ^9|(b). The labels for the external 
legs and for propagators 1 to 6 are suppressed. The expression for the extra propagator required for 
tensor reductions is the same for L, D, G and H topologies: p\ = 2q-ki. The extra propagator for 
LPB, F, J and M topologies is p| = 2p-k2- For the D, G, H, J and M topologies, the ' x ' signifies 
that the marked ordinary propagator is not required in the denominator, only in the numerator. 



by the external leg permutation fci <-> /c2i which also takes z <-> 1 



z. 



We always use this 



permutation to map such integrals into those shown in figure |30 



There are three independent external four- vectors in the light-cone integrals, ki, k2 and 
n. Due to momentum conservation, the same number of vectors appear in the two-loop 
four-point integrals needed for 2^2 scattering, ki + k2 ^ + k^. In fact, n may be 
thought of as a fictitious momentum vector /cs, if desired. Then the derivation of IBP and 
Lorentz invariance identities follows straightforwardly from previous work on four-point 
integrals |6^, 15, 20, 19]. The IBP equations for the L topology are derived by considering 



vr 



DI2 ^DI2 Qa^^ 



1 



(6.7) 



where the pf are given in eq. ( |6.4| ), a'^ G {p^^fl^^ are the two independent loop momenta, 
and can be any of the five available vectors, b^^ G {p^,q^,k^,k2,n^}- Thus for each 
value of {I'i}, there are 2 x 5 = 10 IBP equations, which are linear in the L{vi). They can 
be written as 



= {D -2v^-Vi-u^- uq- uq) + 1/4(5 - 3^)4+ - z^53"5+ + 

+ -3-)6+ + ugs9+ , 

= - + i^3(l" - 6^)3+ - u^is -2 -3 + 6^)4+ + 1^5(1- + 3" 

-z.6(l~-3-)6+-z/9 7-9+, 
= -z/3 + 1/5 + z^35"3+ + U4{s + 5^)4+ - 2/53^5+ - 

- 1^6(3" - 5" + 8^)6+ + uq zs 9+ , 



6^ + 8^)5" 



(6.8) 

;+- 

(6.9) 

;6.io) 
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= - 1/5 - i^3{s - 4" + 5^)3+ - 2^45^4+ + 1/54^5+ + 

+ i^eil- -2- + 4- -5- + 8^)6+ + 1/9(1- z)s 9+ , (6.11) 

= 2/3(5 - 9-)3+ - 1/49^4+ + z/5((l - z)s - 9^)5+ + i^eis - 7" - 9^)6+ , (6.12) 
= j/i - z/6 + 1/1(3- - 6-)l+ - 1/2(5 - 1- - 4- + 6-)2+ + - 3^)6+ + 

+ z/7(s - 9^)7+ - 1/3(3- - 5^)8+ , (6.13) 

= D - 2i/i - 1/2 - z/6 - z/7 - i/8 + z/2(s - l-)2+ - 1/6(1- - 3-)6+ , (6.14) 

= -z/i8-l+ - z/2(s + 8-)2+ - 1/6(3" - 5- + 8-)6+ + z/7 zs 7+ , (6.15) 
= i/i - z/2 + z/i(s - 2- + 8-)l+ + U2{1- + 8-)2+ + z/6(l- - 2- + 4- - 5- + 8-)6+ + 

+ 2/7(1 - z)s7+ - Z/8S8+ , (6.16) 

= -z/i7-l+ + z/2(s-7")2+ + z/6(s-7- - 9-)6+ + 2/3 8+ . (6.17) 

Here are operators taking z/j — > z/j it 1; for instance, in eq. ( |6.8| ) the expression —2/53-5+ 
is shorthand for the term —v^L{vi, z/2, z/3 — 1, 2/4, 1/5 + 1, ^6, 2/7, fs, 1/9). 

Three Lorentz invariance identities for each value of {z/j} can also be derived, by 
requiring 

where 





— k^k2 


^2 ^1 ' 




= k'^n" 


- n^'kl, 


^31/ 


= ki^n" 


- n^k^. 



(6.19) 



The three identities are 



= (1/2 + 1/4 - - i^s) - 1^2(5 + 1" + 2 8-)2+ - 1/4(5 - 3- + 2 5-)4+ + 2/53-5+ + 
+ z/7(zs + z 1- - z 2- + 8-)7+ - z/9(zs - (1 - z)3- - z 4' + 5^)9+ , (6.20) 

= z(i/5 - z/7 + z/g - 1/9) + z/2(7- + 8-)2+ + 1/4(5 - 3- + 5- - 9-)4+ - 

-z/5z3"5+ + z/9zs9+, (6.21) 

= z/2 + 2/4 - (1/5 + 2/8)z - (2/7 + 2/9)(l - Z) - 1/2(5 + 1- - 7- + 8-)2+ - 

- 2/4(5- + 9-)4+ + z/5((l - z)s + 2 4- - 9-)5+ - 

- us{zs + z 1" - z 2- - 7-)8+ + 2/9(1- z)s 9+ . (6.22) 

There are analogous sets of equation for the other integral topologies: LPB, D, G, H, F, 
J and M. 

The next step is to solve the linear system of IBP and Lorentz equations for each 
topology. We use a Gauss elimination algorithm first introduced by Laporta |23]. We have 
used a customized version of this algorithm written in MAPLE |^9| and FORM pO| . 
After performing the reductions, we obtain 13 master integrals for the g ^ gg problem: 

2 2 

L(l, 1, 1, 1, 0, 0, 0, 0, 0) = Spec(5) = (^3^ ' 
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L(l, 
L(l, 
L(0, 
L(l, 
L(l, 
L(0, 
L(l, 
L(l, 
L(l, 
L(l, 
LPB{0, 
F(0, 



0,0,1, 
1,0,0, 

1,1,1, 
0,0,1, 
1,0,0, 
1,1,0, 
1,1,0, 
1,1,1, 
1,1,1, 
0,0,1, 

1,1,1, 
1,1,1, 



0,1,0,0,0 
1,1,0,0,0 
0,1,1,0,0 
1,1,1,0,0 
1,1,1,0,0 
1,1,1,0,0 
1,1,1,0,0 
1,1,1,0,0 
1,2,1,0,0 
1,1,0,0,1 
1,2,1,1,0 
1,1,1,1,0 



LPB{0, 1, 1, 0, 0, 1, 0, 0, 0) = Sset(s) , 
Btri(s) , 
Fish(s) , 

-LPB{0, 1, 1, 1, 0, 1, 0, 1, 0) = Wedge(z, s) , 

LBtri(2:, s) , 

WedgeF(z, s) , 

Zig(z,s) , 

Ptri(z, s) , 

Ptri2(z,s) , 

-LPS(0,l,l,0,l,l,0,l,0)U^i_, = -LPBWedge(z, s) , 
LPBDtri2(z,s) , 

FDtri(z,s). (6.23) 



(Btie(s) is an integral with two independent one-loop triangles, which is trivially related 
to Spec(s).) 

The master integrals are depicted in figure ^ The first three of the master integrals 
have no light-cone denominators present, so they do not depend on z and were encountered 



long ago in the computation of the quark form factor |91]. The next master integral, Fish(s), 
has a light-cone denominator, but nevertheless does not depend on z. The remaining master 
integrals depend on z. The next six come from the L topology, although they can also 
come from other topologies. In the L topology labeling the have light-cone denominator 
7, but not 9, present. The next master integral after that, LPBWedge(z, s) can come from 
either the L or LPB topology, but in the L labeling it has light-cone denominator 9, but 
not 7, present. Finally, the last two master integrals come from the LPB and F topologies 
and have two light-cone denominators present, 7 and 8. 

In intermediate steps, two other master integrals appear, J(0, 1, 1, 1, 1, 1, 1, 1, 0) = 
JDtri(2;, s) and M(0, 1, 1, 1, 1, 1, 1, 1, 0) = MDtri(2;, s). However, there is a partial fraction 
equation relating FDtri, JDtri and MDtri. Consider the identity 



1 



1 1 



+ 



-1 



{p + q) ■ n q - n p ■ n q ■ n {p + q) ■ n p ■ n 



It implies, after inspecting figure |30|, that 



JDtri(z, s) = FDtri(z, s) + MDtri(l - z, s). 



(6.24) 



(6.25) 



After using this relation to eliminate JDtri, we find that the coefficient of MDtri cancels 
out of all g ^ gg splitting amplitude expressions. 

6.3 Differential Equations for Master Integrals 

Since s is the only dimensionful scale entering the master integrals, their dependence on s 
is fixed simply by dimensional analysis. The z-dependence of the master integrals can be 
determined by differentiating their Schwinger-parametrized form with respect to z. This 
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Spec(s) X^3G3 — Sset(s) 




Btri(s) 




Fish(.) XX?- 



Wedge(z, s) 




LBtri(z, s) 




WedgeF(2;, s) 




Zig(z, s) 




Ptri(2;, s) 




Ptri2(2;, s) 




LPBWedge(z, s) 




LPBDtri2(z, s) 




FDtri(z, s) 




Figure 31: The planar master integrals required for computing two- loop leading-color splitting 
amplitudes. An arrow indicates a light-cone propagator, a dot indicates a doubled ordinary propa- 
gator, and a 'x' indicates an omitted ordinary propagator. 

produces an integral where some of the Vi and the dimension D have been shifted by integers 
(D is always shifted by an even integer). The effect of shifting D can also be converted 
into a shift of z/j. Integrals with shifted z/j are again reducible to master integrals. Thus 
we first derive 'dimension-shifting' relations, and then differential equations in z, for the 
master integrals. 

Consider a case where all the are non-negative. Schwinger parametrization is 
achieved by inserting (for i/j > 0) 

1 (-1)-' 



(6.26) 



into the loop-momentum integrals, Wick rotating, and then performing the integrals over 
p and q, to arrive at the generic result 



1 ^(^'.(-irf^) X [^(*.) 



-D/2 



exp 



{-s)Qx{z,ti) 



(6.27) 



Here X stands for one of the integral topologies, riu is the number of positive I'i, and 
Qx{z,ti) is a cubic polynomial in the ti. Also, 



'^p^q ~l~ TpTpq -\- TgTi 



(6.28) 



where Tp, Tq, Tpq are the sums of Schwinger parameters along the lines carrying loop 
momenta p, q, p + q, respectively. For the L, D, G and H integrals, 



?p — ^3 + *4 + *5 ) Tq — ti + t2., 



'-pq 



t6- 



(6.29) 



-54- 



For the LPB, F, J and M integrals, 



Tp = t3 + t5, Tq = ti+t2 + U, Tpg = t6. (6.30) 

From eq. ( |6.27| ) we see that shifting D ^ D — 2 is equivalent to inserting a factor of A 
into the Schwinger-parametrized result. Breaking up A into monomials in the Schwinger 
parameters ti, these integrals can be rewritten using shifted indices I'i. For example, 

Ptri(z, s)\d^d-2 = L{1, 1, 1, 1, 1, 1, 1,0, 0)|d-.d-2 
= L(2, 1,1, 1,1, 2, 1,0,0) +L(1,2, 1,1,1, 2, 1,0,0) +L(1, 1,2, 1,1, 2, 1,0,0) + 
+ L(1, 1,1, 2,1, 2,1, 0,0) +L(1, 1,1, 1,2,2,1, 0,0) +L(2, 1,2, 1,1, 1,1, 0,0) + 
+ L(1, 2,2, 1,1, 1,1, 0,0) +L(2, 1,1, 2,1, 1,1, 0,0) + L(l, 2, 1, 2, 1, 1, 1, 0, 0) + 
+ L(2,1,1, 1,2,1,1, 0,0) +L(1, 2,1, 1,2,1,1, 0,0). (6.31) 

The latter combination of integrals can be reduced to a linear combination of Ptri(z,s), 
Ptri2(^;,s), and the eight master integrals preceding them in eq. (|6.23|) or figure We 
have worked out such dimension-shifting relations for all of the master integrals, but we 
refrain from presenting them all here. We shall give two examples in eqs. (|6l^ ) and (|6ll ) 
below. 

Once the z-dependence of Qx{z, ti) in eq. ( |6.27] ) is known, the derivative of X(i/j) with 
respect to z may be computed from the parametric representation (6.27). The polynomial 
Qx{z,ti) is always linear in z, 

Qxiz,ti) = Q''xiU) + zQ'xiti), (6.32) 

where and are independent of z. For the integrals encountered in the calculation, 
the Qx{z,ti) have a definite sign for < z < 1. Hence all the master integrals are real 
in the time-like region, apart from the overall factor of (— s)~^^. For the first ten master 



integrals in eq. ( 6.23 ), Q\{ti) = t^t^tj, unless any of v^, fg, V7 vanishes. In these cases, 
Q\{^i) — ^"^^ the master integral is independent of z. Thus Fish(s) is independent of z, 
despite containing a light-cone denominator. For the non-trivial cases where v^uqUj > 0, 
using eq. ( |6.27 ), the derivatives of the master integrals are simply obtained by inserting 



(—s)Q^/A into the Schwinger-parametrized forms. (There is an extra minus sign from 
Wick rotation.) For example, 

— Ptri(z, s) = — L(l, 1, 1, 1, 1, 1, 1, 0, 0) = (-s)L(l, 1, 1, 1, 2, 2, 2, 0, 0)\d^d+2 ■ (6.33) 

The shift D ^ D + 2 comes from the factor of 1/A. The effect of this shift is found by 
inverting the D ^ D — 2 shift computed by inserting one factor of A. 

For the integral LPBWedge(z, s), we have Q\{ti) = ti{t3 + tQ)tj. For LPBDtri2(2;, s) 
and FDtri(z, s), we find Q\{ti) = [t^{U + h) + U{t3 + te)]*? - + ^e) + i5(i2 + t6)]i8- 

Reducing the right-hand side of equations like ( |6.33D to linear combinations of master 
integrals, we obtain the following set of differential equations: 

|-Wedge(z, .) = 4p^Wedge(z, .) + ~ '^f^ ~^'^^\~ Sset{s) , (6.34) 
oz z[l — z) e^s^ z[l — z) 
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— LBtrifz, s) 
oz 

d 

— WedgeFfz, s) 
oz 



^Zig(.,.) 



dz 



Ptri(2;, s) 



^Ptrif 



dz 



X l-3e 

-e— -Latiiiz, s) ; 

z{l-z) ^ ' ^ s z{l - z) 

2(l-20-,:(2-3c) 



Btri(s) , 



(6.35) 



+ 



2z{l - z 
3 1 - 2e 



-WedgeF(z,s) + 



LBtri(2;, s) + - 



2(1 - 

3 (1 - 2e)(l - 3e) 



Zig(^, s) + 



2z{l-z) ' ' 4 es z{l-z) 
(1 - 2e)(l - 3e)(2 - 3e), 



Btri(s) 



2e2s2 z{l-z 
3 e 



-Sset(s) , 



(6.36) 



2s z(l - z) ""^"'"^ ' 4 _ ^ 

3 (1-2.)(1 3.)(2 3.) 
2 e^s'^ 2(1 — z) 



(6.37) 



-1 



3-2z 



3 + 6e-2z^ , , 8e(l-2e)^ (6). x 
Ptri(z, s) + — ^ -PtTi^2 {z, s) + 



+ — Zig(z, s) - f WedgeF(2:, s) + 

s s"^ (1 — Z) 



+ 



2:^(1 — z) 
2e(l-z)(9-4z) 



Wedge(z, s) + ^^^—^^Fish{s) + 
2(1 — 2) 



+ (l-2^)(l-3^KS'-lfo + 5.')Btri(,) _ 



e 2^(1 — 2) 
(1 - 2e)(l - 3e)(2 - 3e) (9 + 262) 



2^ 



12 e 

Sset(s) Btie(s) 

s z 



, (6.38) 



-1 



3-22 



2 (1 - 2 



-Ptri(2,s) + 



^3(l-e)-2(3-2.)2 + 22^p^^.(a)^^^^) 



2(1 - 2) 
Zig(2, s) 



4s (1-2) 



WedgeF(2, s) 



4 (1 - 2) 

ll(l -2e) _ 1 ^^^g^^^^ _ 

4es 2(1 — 2) 2s 2 

3 ^. w X 5(1 - 2e)(l - 3e)^ , 

Fish(s) - „ „ ;\ -^Btri(s) + 



2 s 2(1 -2)^^""'"' 8e2s2 2(l-2 
2(l-3e)(2-3e) (2 - e - 2e2) 



s^ 2(1 — 2) 
3-4e2 

Btie(s) 



-Sset(s) + 



LPBWedge(2, s) = 



2(1 - 2e) 2(1 - 2 

2 + e , , (1 - 3e)(2 - 3e)^ 
— -LPBWedge 2, s - ^ , -^Sset s , 

2(1 — 2) €S^ 2(1 — 2) 



(6.39) 
(6.40) 
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oz 



1 - 2z 



LPBDtril^Vz, s 



z{l-z) 2 

1 - 2e /LPBWedge(z, s) LPBWedge(l - z, s] 



(1-36) A" 



+ 
+ 



(1 - 3e)s z{l - z) 
(1 - 2e)(2 - 3e) (1 - 2z) 



22 (1 - zf 

Wedge(z, s) — Wedge(l — z, s) ) + 
Sset(s) , 



+ 



d_ 
dz 



FDtri(z,s) 



g2 g3 2;2(1 _ 

^ ~ ^ ^-FDtri(z, s) - 2 rZig(^, s) + 



(6.41) 



+ 



z(l-z) 
6(1 - 2e) 



s z(l — z) 
z LPBWedge(l - z,s)- LBtri(z, s) 
1 - 3e, 



2es 



-Btri(s) 



(6.42) 



In eqs. (6.38), ( 6.3S| ) and ( |6.41D , we have performed a change of basis: we exchange the (4 — 
2e)-dimensional integral Ptri2(2;, s) = Ptri2(2, s)|D=4-2e for its (6 — 2e)-dimensional version, 
Ptri^^^(z,s) = Ptri2(z,s)| D=&-2e, and similarly for LPBDtri2(2;, s). All of the integrals in 
this new basis are well-defined using ordinary dimensional regularization, following the 
same analysis as in section 3^. The relations required to change the basis are examples of 
the dimension-shifting relations discussed previously, 



Ptri5^^(z,s) 



3(l-z) + e(9-8z) 
2(l-2e)(H-2e)z(l-z 
(1 - 3e)(2 - 3e) 



-Btie(s) 



9 - 472 + 38z^ + e(72 - 373z + 270z^) + 



k3(l + e)(l + 2e)2s3^2(i_^) 

+ e^(171 - 964z + 652^^) + 2e^(63 - 379z + 2462^)1 Sset(s) 



1 -3e 



16e2(l + e)(l + 2e)2s2^2(i_^) 



54 - 92z + Uz^ + e(7 - 8z)(54 - 31z) + 



1459z + 578^2) + 8e^(72 - 130z + 53^^) Btri(s) 



+ 



3 - 42 + e(9 - lOz) 
4(l-2e)(l + 2e)sz(l-z) 
1 - z 



4(l-2e)(H-e)(H-2e)sz2 L 
1 



8e(l + e)(l + 2e)sz2(l-z) 

+ 6^(288 - 403z + 126z^)l LBtri(2, s) + 



Fish(s) + 

9 - 82 + 6e(9 - 7z) + €^(63 - 46z) Wedge(z, s) 
2(27 - 37z + 9^2) + 3e(90 - 125z + 36^^) + 



+ 



8(l-2e)(l + e)(l + 2e)s(l-z) 
1 

8(l-2e)(H-e)(H-2e)z(l-z) 

+ e^(120 - 145z + 42z^)l Zig(z, s 



2z - 3e(9 - Sz) - €^(45 - 34z) WedgeF(z, s) 
2(9 - llz + 32^) + e(102 - 125z + 362^) + 
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s(3-2z + e(9-8z)) 
8e(l - 2e)z(l - z) 



Ptri(2:, s) 



s2 (1 + e) (3 - 2z) 
8e(l-2e)(l + 2e)(l-z) 



Ptri2(2,s) , 



(6.43) 



and 



LPBDtr4^^(z,s) 



2e(l 



;i-2e)(2-36) 



4e3(l + 2e)2s3 22(i_^)2 

5z(l-z)-2e(l-12z(l- 
2z) . . . 1 



4e2fl 



2(l + 2e)(l 
1 - 2e 



3e) s 



rWedge(z, s) 



-z + 2e(l 



8z{l-z)) Sset(s) 



2z) 



o2 



3e) s z 



LPBWedge(z,s) 



2(l + 2e)(l 
1 -2e 



;■ Wedge (1 



1 + e 



e(l-3e)s(l 
LPBDtri2(z,s) . 



- 3e) s (1 - z): 
LPBWedge(l 



z,s - 



(6.44) 



6 (l-3e)(l + 2e)(l + 3e)' 

We solve the differential equations ( |6.34D -( |6.42 ) as Laurent expansions in e, beginning 
with the simplest cases, such as eq. ( ^.34D , which only depends on Wedge(z,s) and the 
previously known integral Sset(s). In a given differential equation (or a coupled pair of 
equations, in the case of eqs. ( 6.3^ ) and ( 6.39D ), we start with the most singular terms in 
the Laurent expansions, and proceed until we reach the order required for the expansion of 
the splitting amplitudes through C(e^). This order corresponds to transcendental weight 
4, where each ln(x) or vr has weight 1, and Li„(a:) and (n (Cn = Lin(l)) have weight n. 
We insert into the differential equation an ansatz which is a linear combination of such 
functions, where x can he z, 1 — z, or —(1 — z)/z. We adjust the coefficients in the 
linear combination until the differential equation is satisfied. The constants of integration 
were determined by different techniques. In some cases the integrals could be performed 
analytically, at least for one value of z. It is also possible to evaluate the Schwinger 
parameter integrals numerically by Monte Carlo integration. To fix the constant, it is 
often sufficient to check that for a certain limiting value of z, z = or z = 1, an integral 
remains finite. 

We first give the Laurent expansions for the previously known master integrals, which 
constitute the inhomogeneous terms for the first set of non-trivial differential equations: 

^_^^_i_,r(i + 6)r2(-e)i2 



Btie(s) 



r(l - 2e) 



i-s) 



-2e 



Sset(s) = (— s 



,l-2e 



1 , C2 , 14C3 , 21 ■ 
r(-l + 2e)r3(l-e) 



+ Oie) 



r(3 - 3e) 



-2e 



(1 - 2e)(l - 3e)(2 - 3e) 



1 

2^ 



Btri(s) = -(-s) 



_2,r(6)r(26)r2(i-6)r2(i 



16 
~ "3 

26) 
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(6.45) 



+ O(e^) , (6.46) 



r(2 



i-s) 



-2e 



(l-2e)(l-3e) 



- 2e)r(2 
1 C2 



2e^ 



+ 



3e) 
13 



41 , 



+ Oie^). 



(6.47) 
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For the remaining integrals we just give the Laurent expansions, omitting an overall factor 
of exp(— 27e). The z-independent integral Fish(s) can be performed analytically, and the 
result is 

i-s^'l (2 , 7C3 55 



Fish(s) = -- 



2e2 + 27 + T^^ 



(6.48) 



The remaining, ^-dependent Laurent expansions are 



Wedge (z, s) 



LBtri(z, s) 



Zig{z, s) 



-2e 



s{l-z) 



2e3 2e 



16 



-Cs 



ln(z) + 2e Li2 



+ 4e2 Lia^-i-^) +8e^ Li4(- 



1 - z 

z 

1 - z 

z 



(6.49) 



(-s)-2^ f 1 Inz 1 



s(l-2e)|2e3 2e2 ^ 2e 



Li2(l-2;)+ln^2; + C2 



+ 



+ Li3(z) + ^Li3(l-z) + 

1 13 16 

+ - ln2 2 ln(l -z)--\n^z- -(2 Inz- —(3 + 



+ e 



-Li4(2:) - ^Li4(l - z) - Li4^-^— - 

11 16 23 

+ - InS - - In^ zln(l - z) + yCslnz - — C4 



-^(Li2(^)+ln2ln(l-2)-21n2 2) + 



(6.50) 



WedgeF(2;, s) = — 



1 

+ - 
e 



L!!{^(Li2(i-.)-c2) + 

-3Li3(^) + ^(Li3(l - z) - Cs) - 2(Li2(l - z) - C2) Inz 



^ln2 2ln(l -z) 



llLi4(z)- -Li4(l-^)-3Li4(- 



l - z 



+ 



+ 4 1nz(^Li3(l-z)-C3 + 2Li3(z)j + 

+ Li2(z) ^(Li2(z)-C2) -21n2z + 31nzln(l-z) 

-^In^ z + lln^ z ln(l - z) + | In^ z In^ (1 - z) - 
8 6 2 

- ^C2 (in' ^ + lnzln(l - z)) - ^(4^ 
i-s) 



(6.51) 



-2e 



2; 



1 In z 

+ 



1 



+ 



4e4 ' 2e3 2e2 
7 



— (ln^z+-C2) + 



1,3 7^ , 
-ln''2;+ -C2ln2; 



-7tC3 



+ 
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+ 6 Lu{z) - 6 Li4(l - z) - 6 U4 (--^) + 6 Li3(l -z)lnz + 



+ 3Li2(z 
3 



-Li2(z) +lnzln(l - 2;) + C2 
/13 



^ln^z + ln^zln{l-z) + 



+ - In^ z ln2(l - z) - C2 ( Y ln^z-3\nz ln(l - z)) 



4 487 



(6.52) 



Ptri(z, .) = -^^1 ^ - + ^ (Li2(l - .) + 2 In^ . + ^(2) + 



+ 



8 Lis (z) — 3 In 2; Li2 (z) 
4 



25 



- - In'* 2; + In^ 2;ln(l - - 5(2 Inz Ca 

3 6 

- 14Li4(z) - 8Li4(l - z) - 8Li4(-^— ^) - Li3(l -z)lnz + 
5 



Li2(2) 
1 



U2{z) - 51n2;ln(l - z) + 31n^ 2; + 6C2 



2 

2, , , . 5 



+ 



Ftvif \z,s) 



1 1 



+ - InS + - In^ zln(l - -z) - - In^ zln2(l - z) + 
3 3 2 

/ \ 28 5 1 

+ C2 (in^ 2; + 6 In 2 ln(l - z)j + y Cs In 2 + -(4 L 

4Li4(2) - 3Li4(l - z) - AU^l^-^-^^ + 

+ 4(^Li3(l - 2) - Cs) ln2 + Li2(2)(^Li2(2) + 21n2ln(l - 2)) 
z In^ z - I In z ln(l - z) - \n^{l - z) + 2(2) - 4(4 

QLuiz) - 2 In 2 (Lis (2) - Ca) + Li3(l - 2) ln2 + 
+ ^Li2(z)(Li2(2) + 21nzln(l - z) - 2(2) + 



(6.53) 



In^ 



1 - z 



+ 



1 19 
+ - In^ z ln2(l - z) - C2 In z ln(l - 2) - — C4 



(6.54) 



LPBWedge(z, s) = 



i-s) 



-2e 



1 ln(l - 2) 1 



2s(l-2)(l-2e) e3 



+ -(Li2(z) + -ln2(l-2;)-C2) + 



1, ,0.. . 1, . 35 



+ Li3(2) + Li3(l - ^) + § ln2 ln2(l - z) - - ln\l - z) - y Cs + 



+ e 



Li4(-^Y^) -C2(Li2(z)+lnzln(l-z)- ^In^z) + 
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+ l(lnz - ln(l - z))' + |c3ln(l - z) - 



(6.55) 



LPBDtr4^^(z,s) 



-2e 



1 



2s2^(i _ z){l - 3e) I e 



In z Lis (z) + ln(l - Li2(l - + 
+ ^ In 2 ln(l - z) (\n z + ln(l - z)j 



+ 3 



(In z + ln(l -z))( Lis (z) + Li3(l - z) - Ci 



U2{z) Li2(l -z)-- (ln'^zU2{z) + ln2(l - z) Li2(l - z) 
1 / 9 

- In z ln(l - z) (in^ z + In^ (1 - z) - - In z ln(l - z) + I2C5 



(6.56) 



FDtri(z, s) 



(-S) 




z^ 


1 




+ - 


3 


e 





2e 



1 



In z 1 



g3 g2 



C2 



^-2^ + ^ 3Li2(l-2) + 21n^z-^ + 



3Li3(z) +6Li3(l -z) + 31nzLi2(z) 
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-ln-^z + - In^ zln(l - z) - 6(2 Inz- —(3 



+ 



+ 6 Lu{z) + 3 Li4(l - z) - 9 U4 
-31nz(^Li3(z)+Li3(l-z)) + 



1 - z 



+ -Li2(z) Li2(z) -ln^z + 21nzln(l - z)) + — In"^ z 



7 



24 



2 In^ z ln(l - z) + ^ In^ 2 ln2(l - z) + ^lnz'^ + ^(s In z - 



34 

y 



41 



(6.57) 



The differential equations ( |6.34D -( p.42 ) are valid to all orders in e. Whereas eqs. ( |6.49 )- 
(6.57) only give the solutions through transcendental weight 4, we have obtained the so- 
lutions through weight 5, up to integration constants. Beyond weight 4, the ordinary 
poly logarithms Li„ are insufficient to describe the solution space. Instead, one can use the 
harmonic polylogarithms (HPLs) |]87| , 20 1, denoted by H(m^; z). For transcendental weight 
w the vector is a string of w entries, which (for our application) can only take on the 
values or +1. The HPLs obey the differential equations. 



d_ 

dz 



H(m^; z) = f{a; z)R{rhw-i; z) , 



(6.58) 



where a = ruyj is the leftmost component of m^, m^_i is obtained from rhyj by omitting 
that component, and 



/(O;^) 



1 



/(+i;^) 



1 

1 - z 



(6.59) 
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These 2"" HPLs suffice because the only true singularities in z on the right-hand side of the 
differential equations are those given in eq. ( |6.59 ). Actually, the differential equations ( 6.3^ ) 



and ( 3.39! ) for Ptri(z, s) and Ptri2^^(2:, s) contain factors of 1/(3 — 2z^ on the right-hand 



side. However, these factors are artifacts of the change of basis ( |6.43| ) from Ptri2(z, s) to 
Ptri5^^(z,s). 

The 32 HPLs at w = 5 can be written in terms of ordinary logarithms and polylog- 
arithms, plus three more functions. Although a few of the z-dependent master integrals 
(Wedge, LPBWedge, and LPBDtri^°0 do not require the three additional functions at the 
u) = 5 level of their expansions, the generic master integral requires all of them. 

7. Splitting Amplitude Results 

With the integrands obtained as described in section ^, and using the integrals obtained in 
the previous section, we can express the results for the splitting amplitudes in a variety of 
gauge theories. We will present results for QCD as well for = 4 and M = \ super- Yang- 
Mills theory. For QCD, both the integrands and the combination of master integrals are 
too lengthy to present here. We will present only the results expanded in e through O(e^). 
The corresponding expressions for the AA = 4 maximally supersymmetric Yang-Mills 



(MSYM) theory are relatively simple, and we shall present them in detail in section 7.1 
We also give the results expanded through ©(e*^), which we presented previously |5^. The 
MSYM results also serve as useful building blocks for the splitting amplitudes in AA = 1 
super- Yang-Mills theory. In section |7.2| we give the M = \ results in terms of master 



integrals, and expanded in e. The M = \ results in turn form useful building blocks for 



representing the QCD results in section 7.3 



In the numerator algebra for the loop-momentum polynomials, the dimensionality of 
the metric r/^,y appears, 

< = I), = 4 - lebn. (7.1) 

As mentioned in section setting = \ defines the HV scheme, which has 2 physical 
states for external gluons, but 2 — 2e physical states for internal gluons. Setting = 
defines the FDH scheme, which has 2 physical states for both internal and external gluons, 



matching the number of fermionic states, and which preserves supersymmetry [^, 70|. We 
therefore present the splitting amplitude results for supersymmetric theories in the FDH 
scheme. 

The QCD results are given in both the HV and FDH schemes. The CDR scheme is 
often used in computations of unpolarized cross sections, or amplitude interferences. It has 
2— 2e physical states for external as well as internal gluons. The HV-scheme results could be 
converted to the CDR scheme by including splitting amplitudes for gluons carrying epsilonic 
helicities. In practice, the CDR results for the scattering amplitude interferences required 
for unpolarized cross sections agree, through two loops, with the sum over helicities of the 



HV results, after the infrared singularities are removed using the Catani formula] 15, |33| . 
The tree and one- loop amplitudes in the Catani formula are different in the two schemes, 
but the finite remainders agree. Thus the collinear limits of the finite remainders in the 
CDR scheme should be obtainable from the HV-scheme results in section %. 



-62- 



The results we will present in this section are for bare ('unrenormalized') splitting 
amplitudes. We will discuss their renormalization in section ^, where we also present a 
discussion of the collinear behavior of remainders after subtraction of infrared divergences. 

7.1 M = A Super- Yang-Mills Theory Results 

In super- Yang- Mills theory, scattering amplitudes are heavily constrained by supersymme- 
try Ward identities (SWT) |]92| , and this has implications for the splitting amplitudes. 
First of all, amplitudes with only or 1 negative-helicity gluon vanish at any loop order L 
in any super symmetric theory []92|, 



^(L),susY(,±^2+,,...,n+) = 0. (7.2) 

Next consider the color-ordered amplitudes with two negative-helicity gluons, 

^(L),susY(i+^2+,...,z-,...,r,...,n+), (7.3) 

known as maximally helicity-violating (MHV) amplitudes. In MSYM, the M = A SWI 
imply that eq. (|7.3| ) is completely independent of the cyclic position of i and j, up to a 
trivial overall spinor-product factor of (ij)^ ^3|. In = 1 super- Yang- Mills theory, this 
relation holds at tree level, but is violated at one loop. 

What are the implications for supersymmetric splitting amplitudes? First consider the 
vanishing amplitude with 1~ in eq. ( [7.21 ). Let two cyclicly-adjacent legs a and h become 
collinear, so that the amplitude factorizes on the MHV amplitudes, 

X;Split«'^^^^(z;a+,6+) x (l-,2+, . . . , p-, . . . , n+) . (7.4) 

1=0 

Since the amplitudes A^^-i^ ™ ^'^^ non- vanishing, the helicity-flip splitting amplitudes for 
P~ — > a"*" 6"*" must vanish to all orders in a supersymmetric theory. Using parity, we have 

SplitL^^'^^^^(z;a-,6~) =0, (7.5) 



for all L. This result includes the tree-level vanishing (2.4), since tree-level n-gluon ampli- 
tudes are effectively {M = 4) supersymmetric. 

Due to parity, we need to present results below only for the cases where P has positive 
helicity, — > a^"^b^'', which we shall denote by 'AaAf,'. Two of these four cases, — h and 
H — , are related to each other by exchanging legs a and b, which also exchanges z <-> 1—z. As 
just noted, the case vanishes in any supersymmetric theory. The two independent non- 
vanishing supersymmetric splitting amplitudes are for P+ a^b^ (++) a-nd ~^ a^b^ 
(— -|-). Since these two tree-level splitting amplitudes are nonzero, we define loop ratios 
rlf^ as in eqs. ( |2.14| ) and (2.20), or for L loops. 



' s 



Split5J(a^% 6^') = rf^^^'^'iz, Sab) x SplitL°l(a^% 6^") . (7.6) 

In = 4 supersymmetry, the +-|- and — h cases are related to each other by the fact 
that the expression ( [7.3[ ) divided by (ij)'^ is independent of i and j. Thus the collinear 



-63- 



limit P + 1)^ is essentially the same as that of 1^2^, up to overall spinor- 

product factors which are the same as at tree level. In other words, a universal rs function 
describes both cases: 



„(L)++,Ar=4, 



Z,S 



^{L)-+,Af=A 



(z, s) = r 



(L),Ar=4 
S 



{z,s) . 



(7.7) 



Using the Bose symmetry relation ( |2.9[ ) for P^ — > a'^b'^, the universal function r^^^"*^ ^(z, s) 
must be symmetric, 

rf^''^=\z,s)=rf^'^=\l-z,s). (7.8) 

This symmetry is not manifest in the one- loop expression ( |2.17| ), but it is easily demon- 
strated using polylogarithm identities. 

After carrying out the unitarity-based sewing procedure for M = 4 super- Yang- Mills 
theory, the integral required for the two-loop g ^ gg splitting amplitude can be written 
as the sum of four terms, corresponding to the L, D, LPB and F integral topologies, plus 
the same four terms with z — > 1 — z: 



„(2),A/'=4. 



{z,s) = i^=\z)+ir\i-z) 



where 

and 
h -- 

Id -- 

Ilpb - 

If -- 



li=\z) = Il + Id + Ilpb + If, 



s2 f d^p d^q 
7r^/2 7r^/2 



.. + (i-.)(pf-p?)+pi + *<i--'>+''^> 



s" I d^p d^q sz + pi TT- 1 

n;;2 



^=2P^ 



s\ f d^p d^q 2 A 1 

2 f d^p d^q ^ 2^ A 1 



(7.9) 
(7.10) 

i=l 

(7.12) 
(7.13) 
(7.14) 



The labeling of the propagators p'^ is given in figure tol and, in the case of the D and F 



topologies, figure |30|. Note that the z <-> 1 — z symmetry ( [7.81) is manifest in eq. (|7.9|) . The 
result (7.9) is given in the FDH scheme, 6r = 0. 



We can see explicitly from these expressions how the appearance of the light-cone 
propagators is restricted by unitarity. In eq. (7.11), whenever pg appears in the denom- 



inator, P2 appears in the numerator. The appearance of P2 eliminates any two-particle 
cut (through lines 1 and 2) or any three-particle cut (through lines 2, 3 and 6) entirely to 
the right of propagator 9, in the terms containing that light-cone denominator. Exactly 
the same argument (after flipping the L integral over by letting ki ^ k2) explains why p^ 
cannot appear in the denominator of the D integrand in eq. ( 7.12| ) along with p^. Finally, 
in eqs. (7.13) and ( [7.14] ), p\ never appears in the denominator. If it had, there would have 
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been a two particle cut (through lines 1 and 2) entirely to the right of propagator 8, which 
is again forbidden by unitarity. 

The result of reducing eq. ( [7.101 ) to master integrals is 



1 



16 3 - 2z 



4sz{l- z)(s Ptri(z, s) + 2(1 - 2e) Ptri^^^ {z, s) 



--z{3-z)[s Zig(z, s) + 3 WedgeF(z, s) 
- 2^—^ (7 - z) LBtri(z, s) - 



4(1 - z){5- 3z) Wedge(z, s) - 4z Fish(s) 
(1 - 2e)(l - 3e)(2 - 3e)(5 - Uz) 



g3 g2 



(l-2e)(l-3e)(16-9z) 



s 



■ Sset(s) — 
Btri(s) -2s(3-4z)Btie(s) 



+ 



+ 2s 



1 - 2e r 



4s^z 



(1 - 2)LPBWedge(z, s) - zLPBWedge(l - z, s) 



-szFDtri(z, s) + (1 - 3e)(l - z)LPBDtrif ^(z, s) 
3 



(7.15) 



Inserting the Laurent expansions ( |6.45| )-( |6.57| ) into eq. ( [7.15 ), and adding the terms 
with z ^ 1 — z, we obtain, 



,(2),A^=4 



{z,s) 



' S 

U/i^X^M 1 1 



4 V -s 



^2e4 
1 

+ - 

e 



7 — <^,-^ Hz)+Hl-z) +^ ln'{z)+ln^{l-z)+C2 + 



1 



2(Li3(z)+Li3(l 



^[ln\z)+\n\l-z) 



+ (ln(z) ln(l -z)-2 C2 j (ln(z) + ln(l - z)) - ^Cs 
- 2(Li3(z) + Li3(l - z) - yCa) (ln(z) + ln(l - z)) + 
+ i (in^(z) + ln4(l - z)) - (ln{z) ln(l - z) - 2C2) (in'(^) + ln2(l - ^)) 



2C2ln(z) ln(l-z) + -C4 



(7.16) 



This result can be cast into the suggestive form |53], 



4'^'"^^^.; z, s) = i (r«'^-^e; z, .)) ' + /(e)r«'^-^2e; z, .) + 0(e) 



(7.17) 



where 



fie) = (^'(l - e) - V'(l))/6 = -(C2 + Cse + C4e' + • • •) 



(7.18) 
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with ip{x) = (d/dx) Inr(x), V'(l) = ~7- The iterative relation ( 7.17| ), along with a similar 
type of relation for the two-loop, leading-color, four-gluon amplitude in = 4 super- Yang- 
Mills theory, has led to an ansatz based on collinear limits for the two-loop, leading-color, 
n-gluon amplitude in = 4 super- Yang-Mills theory ||5^ . At least for the case of maximal 
helicity violation, we expect the following relation to hold. 



M(2)(e) = 1 (M(i)(e))' + f{e) M(^\2e) - + 0{e) 

AL) / .(0) 



(7.19) 



where M^^^(e) 

For n = 4, we found that eq. ( 7.19| ) was violated at order |53|. To be more spe- 
cific, [Mi^^(e)]^ contains at 0{e^) only two-types of Lis polylogarithms, Li^[—s/u) and 
Li^^—t/u), whereas (e) contains in addition the independent function Li5(— s/t). We 



have now examined the z-dependence of the order terms in r^g'^"^~^{z, s). Again we 
find functions not present in the square of the corresponding one-loop quantity. At order 
[^s ^"'^"'^(-2) -5)]^ contains Li5(— (1 — z)/z), but no other Lis functions. On the other 
hand, r^g^' ~^(z,s) contains all three Li5(2;) functions {x = z, 1 — z, —(1 — z)/z) as well 
as all three additional non-Li„ functions required (in the basis we used). The violation of 
eq. ( 7.19| ) at 0{e) is consistent with the intuition that conformal symmetry underlies this 
result, and so it should hold only for D — > 4. 



7.2 M =1 Super- Yang-Mills Theory Results 

For pure Af = 1 super- Yang-Mills theory, the ratio r^^'^^'''^'"''^~^{z,s) now depends on the 
helicity configuration. The two independent non-vanishing cases are P"*" — > a^b'^ (++) 
and a^b^ (~+)- The linear combinations of master integrals obtained are again 

relatively simple, and may be expressed as. 



,(2)++,Ar=i 



{z, s) = r 



{2),J\f=4 



iz,s)+I^+'^=\z)+l 



A 



(1 



(7.20) 



where 



ir^^-\z) 



1 - 2e 
-|- LBtri(2:, s 



z(s Zig(z, s) + WedgeF(z, s) - 2Wedge(l - z,s)) + 



(7.21) 



and 



(2)-+,Ar=i 
s 



(z,s) 



j2)++,Ar=i 



'5 



1 2 s z e (^Ptri^^^ {z,s)- Ptrif ^ (1 - z, s) 



1 -2e 



1 - z 



s Zig{z, s) + WedgeF(z, s] 



(1 - z) (^s Zig(l -z,s)+ WedgeF(l 



z,s 
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z 



2 £2^2 1 



LBtri(z, s) + LBtri(l - z, s) - 

2(2-3e)Sset(s) + seBtri(s) 



\-z 

1 - 3e 1 - 2z 



(7.22) 



These expressions, like the N = \ results, are evaluated in the FDH scheme with r\^^ = 
D, = 4, or 5r = 0. 

At one loop, it happens that the splitting amplitude in = 1 super- Yang- Mills theory 
is the same for — -|- as for -|— |-, 



^(l)-+,Ar=l 



(z, s) = r 



(l)++,Ar=l 



z,s). 



(7.23) 



This relation is spoiled at two loops, but we expect the difference to be finite as e ^ 0, due 
to the one-loop relation. A form of eq. ( [7.221) which makes this property more manifest 
can be obtained by switching basis from Zig(z, s) to 



and defining 



I^=^{z, s) = {l- 3e)Zig^''^(z, s) + ^^WedgeF(z, s) + e sPtri^''^(z, s). (7.25) 

1 2c 



Zigf (z,s) =L(1,1,1,0,1,2,1,0,0) 



_D=6-2e 



(7.24) 



(6), 



Then 



j2)-+,Ar=l 



S 



++'^=\z, .) + \ sz(l^=\z, s) - -z,s)]. (7.26) 



(z,s)=r^ 

The function I^^^{z, s) is finite as e — > 0: Ptri2^^ is finite, so it gives a vanishing contribu- 
tion. There are 1/e poles which cancel between Zig2 and WedgeF, yielding 



2 Li3(z)-C3 -Inz Li2(z)-C2 +0(e 



s(\-z)V 

The result of performing the expansions in e of eqs. (7.20) and ( 7.26] ) is. 



(7.27) 



J2) ++,Ar=i 



z, s = r 



(2),Ar=4, 



+ 



3/^2 



2 \ 2e 



1 1 



2e^ 



+ ^ ln(z)+ln(l-z)-l + 



1 

+ - 

e 



- (^ln(z) - ln(l - z)) \ 2 (^ln(z) + ln(l - z) - l) - ^ 



- 12(^Li3(z) + Li3(l - z)j + 2(^ln(z) Li2(z) + ln(l - z) Li2(l - z)j + 
+ ^ (ln2(z) + ln2(l - z) - 41n(z) ln(l - z) + ^(2) (ln(z) + ln(l - z) 

- 2(ln(z) - ln(l - z)) ' + 4(ln(z) + ln(l - z) - l) + ^Ca - C2 | , 



(7.28) 
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S2)-+,M=l 



{z,s) 



j2)++,Ar=i 



{z, s) + 



+ ^|2(Li3(l -z)- C3) - ln(l - z) (Li2(l -z)- C2) - 
2(Li3(z) - C3) - ln(z) {Lh{z) - C2) 



1-z 



(7.29) 



7.3 QCD Results 

We now present the g ^ gg splitting amplitudes in QCD, for a general value of the 
regularization-schcmc parameter Sr. First we introduce some auxiliary functions describing 
the fermion-loop contributions and residual dependence on 5r. 

The functions /l describing the leading-color fermion-loop contributions to the QCD 
splitting amplitudes are similar to the differences between the J\f = I and J\f = 4 splitting 
amplitudes, so we write them as 

2 



friz,s) 



9 



S2)++,JV=1 



+ 



1 ffl 



2\2e 



12 V-s 



1 



zil-z) 



(z, s) - 
1 



(2),Ar=4 



{z,s) 



+ 



6 



+ 



1 

+ - 
e 



17 

18 



(1 - z) (ln{z) + ln(l - + ^) + 1 ^ln{z) + ln(l - ; 
+ z{l - z) (ln(z) - ln(l - ^)) ' - y (in(^) + ln(l -z))-^^ + ^ 
- ^ (ln(z) - ln(l - z))' + - z) Hz) +{^+z) ln(l - z) - 



3 

C2 _ 151 
6 54 



(7.30) 



and 



frizes) 



4'-^'''-\z,s)-r^S 



+ 



+ 



12 



2 \ 2e 



+ 



1 1 

z(l + z) 



^'^'^=\z,s) 
i(ln(z)+ln(l-z))- 

2(Li3(z) - Ca) - ln(z)(Li2(z) - C2) 



17 
18 



+ 



(1 - z)^ 

2(r^N^-^)-2^^) + ri 

- (ln{z) - ln(l 



3 V 



ln(z) 



+ -(ln(z)+ln(l-z; 



The functions /^l describing the subleading-color fermion-loop and /2 describing the 
double fermion-loop contributions are very simple: 



1 



(7.32) 
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fs^{z,s) = 0, 



2 \ 2e 



1 16 
e 3 



(7.33) 
(7.34) 

f^+{z,s) = 0. (7.35) 
There are also some functions describing the residual Sr dependence of the QCD results: 

A++{z,s) = A-+{z,s) + ^z{l-z)(^J^y l_ln(^)_ln(l-^) + ^ , (7.36) 
A-+(^,.) = ^(zi)''!"^ + J(ln(^) + M1 - ^) - ^) - [Hz] - HI -z)y + 

+ ^(ln(.) + ln(l-.))-|-f|. (7.37) 

Then the unrenormalized two-loop QCD splitting amplitude rs factors are given for both 
++ and — h by 



„(2)a,QCD _ (2) a, 



a = ++ or — h. 
For the helicity-flip case, — a~b~ ( ), we similarly define 



(7.38) 



^(ln(z) + ln(l-z))-2(^ 



ln(z) ln(l-z)\ 187 Sr 



+ 



— z 



6^3 



fsL 

A-{z,s) 



18 V -s 



2\ 2e 



12 



2 \ 2e 



1 16 
e 3 



1 2 
- - ln(2;) - ln(l - z) + - 

C o 



(7.39) 
(7.40) 

(7.41) 
(7.42) 



Then the helicity-flip splitting amplitude in QCD is 



Split 



(2),QCD. . 



(z;a ,b ) = \/ z{l - z) 



jab) 
[abfl 



1 



N/ 



+ SrA- 



(7.43) 



Wc remark that the functions fsL and /2 arc so simple because no light-cone projectors 
on internal propagators arc required to compute them. Clearly /2, representing the double 
fermion-loop {Nj'^) contribution, has no internal gluon lines, and only the simple master 
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integral Btie(s) can appear. There is an internal gluon line in the graphs contributing 
to the subleading-color single- fermion loop (Nf^) contribution fsL- However, this gluon 
always appears inside the fermion loop. All such graphs contribute with equal weight to 
fsL, exactly as if the gluon were a photon. If we think of this gluon as a photon, it becomes 
clear that its gauge transformations can be separated from those of the off-shell external 
gluon with momentum kp, and indeed, a covariant gauge could have been used for it. 
Therefore, in this special case we could have used a covariant propagator for the internal 
gluon, instead of a light-cone gauge propagator, which means that only the ^-independent 
master integrals given in eqs. (|6.45|) -( |6.47|) can appear in fsi- We have discussed in detail 
in sections ^ and ^ how unitarity can prevent certain light-cone denominators from ever 
appearing. The fs l terms are examples of how color can occasionally do the same sort of 
thing. 

Equations ( 7.3^ ) and ( 7.43| ) give the dependence of the QCD results on Nf and Nc- The 



results can also be written in terms of general group Casimir constants (after multiplying 
by the extracted factor of A'^^ = C^), using the substitution rules: 



2CATRNf 



^^2{CA-2CF)TRNf, 



{2CATRNff . (7.44) 



To recover pure N = 1 super- Yang-Mills theory, which contains one Majorana fermion in 
the adjoint representation, we set Cp = Ca, Tr = Ca, and Nf = 1/2. Thus Nf/Nc — > 1 
and Nf/N^ —1, so that for 6r = each term in eq. (|7.43| ) vanishes, and all but the first 
term in eq. ( f7.38 ) vanish, as required. 



Motivated by ref. let us examine the degree of transcendentality of the above 
results, where ln(x) or vr has weight 1, and Li„(x) and Cn (Cn = Lin(l)) have weight n. 
Leading transcendentality is 4 units for the 0{e^) terms in the expansion, or more generally 
4 — m units for the 0{e~"^) terms. The result in = 4 super- Yang-Mills theory, eq. ( 7.16| ), 



carries this maximum weight uniformly. The M = 1 super- Yang- Mills results (7.28) and 



( |7.29 ) are equal to the AA = 4 result, plus terms of lower weight. The additional functions 



/£, f§^, and A" entering the QCD results (fL3^ are also of lower weight. Hence the 
maximum-weight terms in the QCD, M = 1 and AA = 4 two-loop splitting amplitudes all 
coincide. It is easy to see that the same result (with 4 — m replaced by 2 — m) is true for 
the one-loop splitting amplitudes to O(e^) (in fact to all orders in e). So both types of 
two-particle final-state contributions to the NNLO Altarelli-Parisi kernel — two-loop times 



tree, and one- loop squared — respect the transcendentality relation of ref. [52|. 
7.4 Continuation to Space-Like Region 

The preceding splitting amplitudes are for time-like kinematics, kp = s > 0. In the 
remainder of this section we briefiy discuss the space-like case, s < 0. In the application 
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of splitting amplitudes to NNLO evolution kernels, the time-like region is relevant for the 
evolution of fragmentation functions. The space-like region is relevant for the evolution 
of parton distribution functions, which play a key role in the prediction of collider cross 
sections. (As mentioned in the introduction, the NNLO kernels for unpolarized space-like 
evolution have been computed very recently |48|.) 

The time-like results can be continued analytically into the space-like region. We take 
the space-like splitting process to be 



(-A^i) 



-Ai 



+ h 



A2 



(7.45) 



with /cp = s < 0. It is obtained from the time-like process by crossing leg 1 into the initial 
state, and leg P into the final state. The physical helicities of legs 1 and P flip under this 
crossing, but we retain the uncrossed labeling. Relative to (— the vectors {—kp) and 
k2 now carry longitudinal momentum fractions of x and 1 — x respectively. Comparing the 
ratio {ki ■ n)/{kp ■ n) between the time-like and space-like cases, we identify z = 1/x. In 
principle, all we need to do to obtain the space-like splitting amplitude results from the 
above time-like formulae is to let s be negative, and substitute z — > 1/x, where < x < 1. 
For example, the tree-level splitting amplitude Split ^''(z; a~^, b^) given in eq. ( |2.5D becomes 
the space-like splitting amplitude for {—ki)^ {—kp)^ ~^ •> which we denote by giving 
Split(°) the argument x. 



SplitL°^(x;a+,6^ 



\/l — X {a b) 



(7.46) 



(Note that in the construction of the space-like Altarelli-Parisi kernel, after squaring the 
x-dependent part of the splitting amplitude there is an additional factor of x; see e.g. 
eq. (6.23) of ref. ||7^.) There are overall phases stemming from the factors -y/l — z and (a b), 
but they can be associated with the external states, and we neglect them in eq. ( [7.46| ). 

At the loop-level, the preceding formute for r^\z, s) have been been written so that 
the logarithms and polylogarithms are manifestly real for < z < 1. To continue them 
to z > 1, one can apply polylogarithm identities so that the only function which is not 
manifestly real is ln(l— z) = ln(l— x)— InxitzTr. For example, consider the one-loop splitting 



amplitudes r^P"^ ^(^, s) in eq. ( 2.17|) , which contain Li2m-i{j^)- These functions develop 



branch cuts for z > 1. In this case, though, we can use the (non- manifest) z <-> 1 



symmetry of r 



(l),^=4. 



z,s 



eq. (^), to let 



1 



z 



in eq. ( p. 171) before substituting 



z — > 1/x. (This amounts to using an infinite sequence of polylogarithm identities.) 
obtain, 



We 



's 



(x,s) 



cr 



1 



-(1 



ne 



sin(7re) 



+ ^2e2™-iLi2™-.i(l 



m=l 



(7.47) 



In eq. ( [^.47 ), the factor {^'^ /{—s)Y no longer contains any in terms, as s is now negative. 
However, in the expansion of [—(1 — x)]~^ we must set ln[— (1 — x)] = ln(l — x) it ivr. 
The presence of imaginary parts in the space-like region is a bit surprising, but it can be 
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traced back to integrand denominators that vanish in the interior of the integration region 
for z > 1 — for example, the denominator of J^^\l — z) in eq. ( |3.19| ), as discussed in 
section 3.3. The appearance of an imaginary part can be verified from the explicit collinear 
behavior of one-loop amplitudes. Prom this verification, one also learns that the sign of 
the in term is ambiguous; it depends on whether the leg color- adjacent to leg 2 is incoming 
or outgoing. Fortunately, this sign is the same for all imaginary terms, and it drops out 
of the interferences required for the computation of evolution kernels. The tt^ factors from 
the product of two iir terms do survive and are unambiguous. (Note that s{l — x) has the 
same negative sign in both space-like and time-like regions, so that the tt^ terms from the 
first term in eq. ( 7.47 ) are the same in both regions.) 

At the two-loop level, eq. ( [7.47 ) and the squaring relation ( [7.17 ) provide a convenient 
way to continue r^g'^"^~^{z, s) to z > 1 (through 0{e'^)). We refrain from giving the explicit 
continuations of the M = 1 and QCD results, though they are straightforward to carry out. 



8. Comparison to Catani's Formula and Finite Remainders 

In this section we compare the pole terms in the two-loop splitting amplitudes to the expec- 
tation based on Catani's formula [^8| for the two-loop infrared divergences of renormalized 
amplitudes, in order to establish their mutual consistency. More importantly, by examining 
the order terms in the splitting amplitudes, we derive relations that the finite remainders 
in the Catani formalism must obey in the collinear limit. 

8.1 Singular Term Comparison 

Catani's general formula includes color-space operators which have a fairly intricate struc- 
ture in the trace-based color decomposition. In this section, to simplify the analysis we 
shall restrict attention to the single-trace coefficients aH^'^ in the n-gluon amplitude, given 
explicitly in eq. (|2.2| ), and to the terms in Catani's formula obtained by making the re- 
placement, 

Nc 

Ti ■ Tj — > — — 1, i,j color-adjacent, 

— > 0, otherwise. (8.1) 

These terms include all the leading-color terms, as well as certain of the subleading-color 
terms (including all of the color dependence in ' at order 1/e). We refer to these terms 
as color trivial. In appendix [A.4[ we perform the full color-space analysis. 

First we need to remove the ultraviolet divergences from the splitting amplitudes pre- 
sented in section 0. The relation between the bare coupling a" (implicitly used above) and 
renormalized coupling as(/u) = (7^(/i)/(47r), through two- loop order, may be written as |^8| 

< S, = as{n) I?' 

where /i is the renormalization scale, and = exp[e(ln 47r-|-^/;(l))]. The first two coefficients 
appearing in the beta function for QCD, or more generally SU(iVc) gauge theory with 
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flavors of massless fundamental representation quarks, are 



ba 



6 



17CI - {IOCa + GCF)TRNf 
6 



(8.3) 



"0 



l)/(2iVc), and Tr = 1/2. In pure M = 1 super- Yang-Mills 
3 



where Ca = Nc, Cp ■ 
theory, the values are 

In = 4 super- Yang- Mills theory, a conformal theory, the values are of course 



0. 



One can define a 'perturbative expansion' of the g ^ gg splitting amplitude as 



Split^(as(^)) = g{p) 



Split 



(°)+^Split«+f^ 



27r 



V 27r 



Sphtg) +0{al{fi)) 



(8.4) 
(8.5) 

(8.6) 



where Split^^ is the L*'^ loop contribution. Equation (8^) is equivalent to the following 
MS renormalization prescriptions at one and two loops. 



Split« = 57^ SphtW-l^SphtW, 



(8.7) 



Split^^^^ 



S-^ Split 



(^)-||.-Spiita) 



3 hi 



£2 Ar2 



4eA^2 



Split 



(0) 



Here Split*-^^ governs the collinear limits of unrenormalized amplitudes, while Split^^ con- 
trols the limits of renormalized amplitudes. All quantities appearing in the remainder of 
this section are renormalized. In particular, A^^'^ refers to renormalized L-loop amplitudes, 
in contrast to the unrenormalized ones used implicitly in previous sections. The renormal- 
ization of the splitting ratios r^^'' follows simply from eqs. ( ^.7| ) and ( |8.8D . We will denote 
the renormalized splitting ratios by 

The full color-space forms of Catani's formula at one and two loops are given in 
eqs. ( A.15| ) and ( A.38| ). The color-trivial terms, defined by the replacements ( |8.lD , are 
obtained by letting the operator I^n'^ iji'^l, in its action on the single-trace term 

^n^''(l, 2, . . . , n) in the amplitude. The factor of is extracted from I^^^ for consistency 
with the normalization of A^^"^; we want both a!^^ and to be independent of N^. as 
A^'c — > DO in the pure-glue case. At one loop, we have 

1 e-^'^(i) 



where 7^ is given in eq. ( A. IT)) 

r(2) 



2r(i 



1 , 2^1 

e2 ^ AT, e 



Similarly, /„ is given by 
+ 



e+^^(i)r(l - 2e) 1 

r(i - e) n'c 



4er(l 



1=1 



2t 



^ + i^K.s)/«(2e) + 

(8.10) 
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where is given in eq. (|A3^ ) and iff ^ is given in eq. ( |A.33| ). 
The singularities of the single-trace coefficients are, 

a{1) ^ f(l) 4(0) , /l(l)fin 

n n n ^ n ' 

4(2) = f (1)4(1) -U f(2)4(0) , 4(2)fin 



Ul) 
^.12) 



Things become a bit simpler if we write the finite remainders in Catani's formula as mul- 
tiples of the corresponding tree amplitudes: 



4^)fi°(l, 2, . . . , n) = fW(.,,) X 40)(1, 2, . . . , n) . 
Equations ( |8.1l| ) and ( ^.12 ) become 



^.13) 

^.14) 
^.15) 



If the tree amplitude vanishes, we cannot perform this step. However, in this case the 
entire analysis is much simpler, essentially equivalent to the one-loop analysis. 

Now consider the collinear limits. We will provide two independent forms for the pole 
terms in e of r^''. The first form is similar to Catani's formula ( 8.10| ) for /^^\ The second 
form retains finite terms, so it can be used to predict the collinear behavior of the Catani 
finite remainders. 

In the first derivation, we note that the term by itself does not have particularly 



enlightening collinear limits, because of contributions coming from 0{e) parts of I,^^'*. 
ever, this obscure behavior is balanced by similar behavior in the In^ A^'^ term in eq. 
suggesting that it is best to combine the first two terms in this equation: 



How- 



+ 



e+^'^(i)r(l - 2e) 1 



r(i-6) 

e-^V'(i) 



4er 1 



E ■ 



(2) 



6) ^ iV2 

1=1 



2f 



+ 



e+^^(i)r(l - 2e) 1 



+ 



r(i-6) 



r -fvRs. 

e 



n H, 



(2) 



.16) 



and 



4e iV2 ^ 

In the last step we added finite pieces, in particular ones proportional to [-Fi^'*]^ 
K^.^F^^\ Two singular terms proportional to (1/e) x 60 -^n^'* cancel against each other. 

Now the combination In^ [e)+Fli'^ is just An^ (e) /An^ , so it has simple collinear limits. 
Indeed, inserting eq. ( |2.14| ) into eq. ( 2.1C| ), we find that 



,(1) * f{i) 



(1) 
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We use the behavior ( 8.17 ) to take the colhnear hmit of eq. ( |8.16| ), and compare the 
result with eq. ( 2.19 ). The terms quadratic and hnear in i^}_i{e) + Pj^}^ belong to An^ x 
Split^''^ and An^ x Split^\ Most of the remaining terms belong to An^ x Split^^ = 
An^ X Split '•''^ xr^g\e). Collecting them, we see that the divergent parts of the two-loop 
splitting amplitude should be given by 



1 



.(1)/ 



1 bom. . e+^'^(i)r(l - 2e) 1 /ft. 



r(l-6) 



N, 



:2i + K^\f^p^2e) + 



+ 



1^ 
4e 



+ O(e0) . 



^.18) 



(2) 

In the second derivation of the pole terms in (e), we shall retain the finite terms. 
Consider first the one-loop case. Equation ( ^l7|) can be rewritten as 



.(1) 



/(I) . 



a\\b 



f(i) _ 
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which allows one to predict the singular terms in r^^ ^ in terms of the collinear behavior of 
/(I). 



.(1) 



e pole 



J(l) 



a\\b 



.20) 



But we can also solve eq. ( p. 191 ) for the collinear behavior Fn\a\\bi using also eq. ( |8.2q ). 
For definiteness, we will assume that a = 1 and 6 = 2. We have then 



.21) 



1112 



where 



^(^^(z,SnP,sp3,s) = r 



■(1) _ ^(1) 
s ' s 



e pole 



^.22) 



Note that the Mandelstam invariants involving the gluons which are color- adjacent to a = 1 
and h = 2 appear, namely gluons n and 3. In evaluating eq. ( ^.22| ), r^^'^le pole is given by 
eq. ( ^.20 ), including all terms at order e''. 

At two loops, the collinear limit as a || 6 of eq. ( ^.15 ) is, using eq. ( |2.19| ), 



a\\b 



■^^^ + WiW + F^'^) + „^ X Split(o) x4°l, 

J a 



f(2) _i_ f(l) /f{l) _i_ T^yi-) \ I T^K^) I 
-'n-l "T -'n-lV-'n-l ^n-l) ^ n-1 



n{2) 



_|_ J. 



X Split(o) x4°l, 



Solving eq. ( 8.23| ) for r^g \ we find 



.23) 



^(2) 



-'n-l -'n-lV-'n-l -^n-lJ 'S V-'n-^l -^^n-l 



(2) ^(l)/fa) 
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Equation (8.24) contains the one-loop finite remainder Fn^^ multiplied by a singular factor 
In\ However, the pole terms in should not depend on any one-loop finite parts. 
Therefore we use eq. ( |8.19| ) to eliminate the In^ Fn^ term. After a little algebra, we obtain 



.(2) 



a\\b 



f(2) 
'n~l 



Now it is clear that the e pole terms in r ' have a universal form. 



(2) 



(1) \2 



e pole 



(1) \ fW 



-'n-l V'5 -'n-1 /'-'n-1 • 



(1) Af(l) 
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We have checked that the pole terms in this expression are equivalent to those in eq. 
for AA = 4 and = 1 super- Yang- Mills theory, and for QCD. We have also verified that they 
agree with the singular parts of the g ^ gg splitting amplitudes from our explicit results, 
eqs. (|7l^ ), (fO^ ), (fOgp , (173^ ) and (ffl^ ). Along with the full color-space discussion 
in appendix A. 4, this shows that Catani's formula for the singular behavior of two- loop 
amplitudes is completely consistent with the collinear limits. 

8.2 Finite Remainders 

The next step is to use the finite parts of the two-loop splitting amplitudes to deduce the 

(2) 

collinear limits of the finite remainders Fn in the color-trivial parts of Catani's formula. 
Again letting a = 1 and 6 = 2, we rearrange eq. ( ^.24|) , with the help of eq. (|8.26| ), to get 
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1||2 
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n-l + 
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.(1) 
s 



t pole 



or 



where 
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J(2)/ 



^^^)(z,s„P,sp3,s) = r 



■(2) _ -(2) 

s ' s 



e pole 
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i.29) 
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Here r^g^\^ poie is given by eq. ( |8.26| ), including all terms at order e''. Equation ( 8.2g| ) 
provides a useful check on finite remainders of two-loop scattering amplitudes, as any two 
external gluons become collinear. 

Now we present the values of ^^^^ and ^^^^ for the various theories we have been 
considering. The values of ^^^^ are: 



^(i),Ar=4 

^(1)++,QCD 
^{1)-+,QCD 



ln(l — z) In 



-SP3 



+ ln(z) In 



-•SnP 



+ ln(z)ln(l-z)-C2 , (8.31) 



(1) -+,M=1 ^ ^{l),Af=4 _ % 



2Nr 



-(ln(z)+ln(l-z)+ln(-s)), 



{l),Ar=4 
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■(ln(z) +ln(l - z) +ln(-s)). 
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To get the proper analytic behavior of such expressions, one should apply the prescrip- 
tion ( p.l3D for logarithms of time-like invariants, after expanding the logarithmic ratios, 
ln((-si)/(-S2)) = In(-si) - ln(-S2). 

Next we the quote the values of For AA = 4 super- Yang-Mills theory, ^^^^ obeys 
an iterative equation in terms of as a consequence of eq. ( [7.17| ), 



t(2),Ar=4 ^ 1 r 
^ 2 



(l),Ar=4 



C2C 



{l),A/'=4 _ ^ 

32 



(8.35) 



For J\f = 1 super- Yang-Mills theory, we first define the auxiliary function 

V{z,snp,sps,s) = -2(Li3(2;) + Li3(l - z)) + ^(lnzU2{z) +ln(l - z)Li2(l - z] 

1 / 3 \ 8^ 

- In z ln(l - z) - -C2 (Inz + ln(l - z) + In(-s)) + —(3 ■ 



24 



lnzln(^^^)(ln(-s„p) +21n(-s)) + 



+ ln(l - z) In 



-SP3 



ln(-sp3)+21n(-s) + 



+ In z ln(l - z) ( ln(-s„p) + ln(-sp3) - In(-s) I + 



+ 2(ln2zln^ """^^ 



j +ln^(l - z)ln 



Then the two functions required in eq. ( 8.2S| ) are given by 



-SP3 
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e'^-^'^=' = ^(2)++'^=! + U 2fLi3(l - z) - C3) - ln(l - z)(Li2(l -z)-C, 



1 - z 



2 Li 



3(z)-C3) -ln(z)(Li2(z)-C2) 



. (8.38) 



Equation ( 8.38| ) follows directly from eq. ( |7.29| ), because eq. ( ^.26 ) for the subtraction term 

_(2) 

'"g |e polo is the same for ++ and — h. 

For QCD, we define one more auxiliary function. 
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In terms of this function, we have 
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We can check some hmiting properties of these results as z ^ and z ^ 1, using simple 
facts about soft limits of amplitudes. For example, in the limit z — > 0, leg 1 becomes soft. 
The helicity of the hard leg should be conserved in the soft limit, and the limit should 
be independent of the helicity of the soft leg. Thus the cases P'^ — > 1"'"2"'" (++) and 



— > 1^2+ ( — h) should behave identically as z ^ 0. The tree splitting amplitudes ( ^I^ ) 
and are the same in this limit (up to an external phase associated with the soft external 
state). Hence the rs factors, and also the Catani-subtracted functions, for ++ and — h 
should behave identically as z — > 0. This property is obvious for = 4 super- Yang-Mills 
theory; the two rg factors are identical for all z due to the = 4 supersymmetry Ward 
identity. For pure Af = 1 super- Yang- Mills theory, one can inspect eq. ( [7.29 ) or eq. ( 8.3^ ), 
recalling that Li„(l) = (n, to see that the quantity in braces indeed vanishes as z ^ 0. 



(2) oi 

For QCD, the factors in section |7.3| are not written in a particularly convenient way 
for checking the limit. However, one can easily compare ^(2)++, QCD -^^ 

(|O0| ) with 

^(2) h.QCD ( 8.41| ). They do approach the same limit as z — > 0, namely the limiting 

behavior of 7*^'-^^. 

For the -|-+ case, the limit z ^ 1 is the same as the limit z ^ 0. For the — -|- case, in 
the limit z ^ 1 the helicity of the hard leg, now leg 1, flips from -|- to — ; hence this splitting 
amplitude should be suppressed. Equation (|2.7D shows that it is suppressed by two powers 
of 1 — z at tree-level. In the = 4 and M = 1 supersymmetric cases, there is no additional 
1/(1 — z) singularity in the rs or ^ factors. (There is an apparent 1/(1 — z) singularity 
in the quantity in braces in eqs. ( 7.29| ) and ( |8.38| ), but again recalling that Li„(l) = Cn, 
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one sees that it cancels.) So, up to logs, the soft behavior is the same at the loop-level as 
at tree-level. In the case of QCD, the limiting behavior of ^^^^ h,QCD -^^ 

( |8ll| ) looks 

quite singular as z — > 1, since powers of 1/(1 — z)^ and 1/(1 — z)^ appear. However, these 
cancel, and the actual behavior is 

^(2) -+,QCD + (8.42) 

^ 12 V Nj 1-z ^ ^ 

Taking into account the (1 — z)^ behavior of the tree-level splitting amplitude, the one 



power of 1/(1 — z) in eq. (|8.42| ) still means that the soft limit z — > 1 for — h is suppressed 
by a power of 1 — z, relative to that of ++, as expected from the helicity flip on the hard 
line. 

For the term corresponding to the helicity-flip splitting amplitude, P"*" 1~2~, we 
should not remove the tree-amplitude factors. Instead we write 



4^)^°^ = Splitg)jz; l-,2-)A«f + 7^0^^ ^(2)-,QCD^(0)^ 



.43) 



where 



^(2)— ,QCD ^ _^ 



'^(i),^=4 ^ ^ ln(l - z) _ 

1 — z z 

(Inz + ln(l — z) + 31n(— sj) H h 



2Nc ' " 18 6 9A^, 



-Nf . (8.44) 



For this helicity configuration, the tree-level splitting amplitude vanishes, so the one-loop 
renormalization is trivial: Split^]j(z; 1~, 2~) = Split^^(z; 1~, 2~). The expectation as 
z — > (or equivalently, z ^ 1) is that ^'•^^ ,QCD g}jQ^u have no power-law 1/z singularity, 
since its prefactor z(l — z) in eq. ( ^.431 ) has only one power of z suppression relative to 
the ++ case. Indeed, in eq. ( p.44| ) ln(l — z)/z is finite as z; ^ 0. 

Equations ( |8.21 ) and ( |8.29| ) govern the collinear behavior of finite remainders Fn^^ for 



the case where only one intermediate helicity A can contribute. If both helicities contribute, 
one needs to sum over the two, according to eqs. (|2.10| ) and ( ^.19 ). Before doing this, it is 



best to multiply back by the tree amplitudes A^}_^{X), because they are different for the 
two terms. 



9. Dressing with Color 

In this section we present the full color structure in the collinear limits. We will do this 
in the context of the trace-based color decomposition discussed in section |2| (and reviewed 
in refs. |67, |l3|). An alternative, color-space language has been used by Catani to pre- 
dict |38| and describe the infrared structure of two-loop amplitudes. In appendix ^ we 
re-express the collinear behavior in the color-space language, in order to demonstrate that 
the divergent parts of our splitting amplitudes are fully compatible with the structure of 
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infrared singularities predicted by Catani |^^. The color dressing here applies equally well 
to renormalized and unrenormalized amplitudes, but in the discussion in appendix all 
quantities are understood to be renormalized. 

At tree level the trace-based color decomposition for an n-gluon amplitude is 

^(0)aia2...a„ ^ ^ Tr(r"-(i) . . . T^-w (cj(l), . . . , cj(n)) , (9.1) 

where the are tree-level color-ordered partial amplitudes and the notation is defined 
below eq. ( |2.2[ ) . On the left-hand side we have suppressed the dependence on helicities and 
momenta but have left the color dependence explicit, since that is what we focus on in this 
section. 

The fully color-dressed tree amplitude has the following factorization property as the 
momenta of legs 1 and 2 become collinear, 

where Split!i^{ is the tree-level splitting amplitude ( |2.3D defined for the color-ordered am- 
plitude, and there is an implicit sum over the helicity A of the intermediate state. Here 
jaxa2ap -g SU(iVc) Structure constant corresponding to our normalization of the gener- 
ators, Tr(r'^r^) = b"-^: 

so that f"'^'^ = ^/2f°''"^, where f"'^'^ is the structure constant with conventional T"" normal- 
izations. 

By inserting eq. (|9.1| ) into both sides of eq. ( |9.2[ ), and using SU(A^c) Fierz identities on 
the right-hand side, we can see that this equation is equivalent to the collinear behavior 
of color-ordered amplitudes given in eq. ( |2.3| ). A given term in the amplitude ( |9.1| ) will 
contribute to the collinear limit only when the collinear legs 1 and 2 are cyclicly adjacent 
in the associated color trace. Alternatively, we can derive the color structure of eq. ( |9.2| ) 
directly from the fact that the only diagrams contributing to the collinear limit are of the 
type shown in figure where a = 1 and 6 = 2. The factor of J"i"2ap precisely the color 
factor of the vertex joining legs 1, 2 and P. 

At one loop the full color decomposition of an n-gluon amplitude is, 

Ln/2J+1 

_4a)aia....a„^ ^ ^ Gr„;,(a(l),a(2),...,c7(n))A„;e(^(l),a(2),...,a(n)), (9.4) 

C=l lj£Sn/Sn:c 

where [x\ is the largest integer less than or equal to x. The leading-color structure factor 

Gr„;i(l,2,...,n) = iVeTV(r'^i ...T'^"), (9.5) 
is just Nc times the tree color factor, and the subleading-color structures are given by 

Gr„y(l,2, . . . ,n) = Tr(T"i . ..T'^^-i) Tr(r''^ . ..T"") , j>l. (9.6) 
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Again Sn is the set of all permutations of n objects, and Sn;c is the subset leaving 
GiCn-c invariant. In the trace-based color decomposition, fundamental-representation quark 
loops will contribute only to An-i, with a relative factor of Nf/Nc compared to adjoint- 
representation gluon loops. 

At one loop, the subleading-color amplitudes A^-c are completely determined in terms 
of the leading-color ones An-i (see eq. (7.2) of ref. The collinear behavior of the 

leading-color amplitudes was given in eq. ( p. 101 ), where the notation ^An ' was used instead 
of '^n;i'- The collinear behavior of the subleading-color amplitudes can be determined from 
the previously-mentioned relation, or from the observation that the two-particle collinear 
limit cannot reduce a pair of traces to a single trace. Therefore the tree-level (n — 1)- 
point amplitude A^^}_^ cannot enter into the limit. This implies that the one-loop splitting 
amplitude cannot appear either, and so the limit is (for c > 1), 

(. . . , , . . .) ^ Split-1(^; 0^^" > 4'-i;c-i(- • • , ^\ • • •) , a, 6 < c , 



4l)(...,a^^6^^...) ^SplitL"i(z;a^%6"'')4\,(...,P\...), a,6>c. (9.7) 



,^,,^(0)...„A„ ,A,, ,(1) 



Legs a and b must be cyclicly adjacent within the same color trace, otherwise the collinear 
limit is finite. In particular, the collinear limit is finite if the two legs lie in different traces. 

We can reassemble these properties of A^}c into a description of the collinear behavior 
of the full color-dressed amplitude, 

spiitL°l(z; lAi^ 2^2) + 

+ N, SplitL'l(2; 1^^ , 2^2) 



^(1) aia2a3...a„ j^jaia2ap 



(9.8) 

where Split ^]^(z; V^'^ , 2'^^) ^re precisely the color-stripped splitting amplitudes (|2.15|) - (|2.17|) . 



Note that the c > 1 contributions in eq. ( |9.7D all go into assembling J^^-i ™ ^Tsi term 
in brackets. We can also derive eq. (|9.8| ) directly from argumentation using the unitarity 
sewing rules (or by using light-cone gauge, ignoring prescription issues). The diagrams 
in figure § correspond to a factorization of the process including color indices. The first 
term in the brackets corresponds to the left diagram, and the second term to the right 
diagram. The color factor in the first term is precisely the tree-level one. The color factor 
in the second term is that of a one-loop vertex, f"-^"-2"-p in the pure-glue case. (Prom 
the diagrammatic point of view, there are three types of color contributions obtained from 
color dressing the diagrams in figure |^. However, all these color factors are proportional to 
Nc f^i^^ap _ Including fermions in the loop gives the Nf/Nc terms in eq. ( 2.15| ).) 



If we examine the cuts of a color-stripped 1 — > 2 splitting amplitude, each contribution 
has the product of a {j + 2)-point scattering amplitude on the left and a 1 — > j lower-loop 
splitting amplitude on the right, where j is the number of particles crossing the cut. (See, 
for example, figure ^.) We can exchange the two final-state particles by reflecting through 
the horizontal axis. This produces a factor of (— l)-^"*"^ on the left and (— l)-^"*"^ on the right 
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because of the properties of the amphtudes under reflection [37, ^. If the two final-state 
particles are gluons we get the same number of Fermi minus signs on the left and right 
from interchanging cut internal fermions, if there are any. Overall we always get a factor 
of —1; that is, the color-stripped splitting amplitude is anti-symmetric when we exchange 
the two arguments, including their helicities and z (1 — z). Bose symmetry then implies 
that the color factor must also be antisymmetric. The color-dressed splitting amplitude 
must be proportional to y^i^aip^ because there is no other antisymmetric invariant tensor. 
This argument holds to all loop orders. 



Note from eqs. (|2ll| ), (|t|) and that there are no subleading-color corrections 



to the pure-glue (A^/ = 0) contributions to the one-loop splitting amplitude Split*-^-*, which 
appears in the full-color collinear limit ([9.81) . To understand this fact, and the factor of 
in front of the splitting amplitude, consider the color factors in a diagrammatic calculation. 
One can use the Jacobi identity to show that the result of the color algebra in any diagram 
is a linear combination of the results for its various possible parent diagrams. It is therefore 
sufficient to consider the parent diagrams. At one loop, for pure-glue contributions, there 
is only one kind of parent diagram, the triangle diagram show in figure ^(a). Now, as far as 
the color algebra is concerned, any triangle subdiagram can replaced by a factor of —^Ca 
times a three-point color vertex. Similarly, any bubble subdiagram can be replaced by a 
factor of —Ca- Thus at one loop we just obtain an overall coefficient of Ca = up to 
A^c-independent factors. 

Now we turn to the two-loop case. The color decomposition generalizes in an obvious 

way, 

Ln/3J+1 L(i+ci+l)/2J 

^(2)a.a2...a„ ^ EE E Gr„;,,;,, Kl) , ^(2) , . . . , a(n)) X 

Cl=l C2>2Cl-l (TG5„/S„;ci,C2 

X ^„;ei;c2(fT(l),^(2),...,a(n)), (9.9) 

where 

Gr„;ci;c2(l,2,... ,n) = Tr(r"i ...T^'^i-i) Tr(r"'=i ...T^'^a-i) Tr(r"'=2 ...T""), (9.10) 

and a runs over the set of permutations 5„, modulo those in 5'„;ci;c2 which leave Gr„;ci;c2 
invariant. We identify Tr(l) = so, for example, 

Gr„;i;i(l,2,...,n) = A^^ ^r(r'^^ . . T'^") , (9.11) 
Gr„;iy(l,2, ... ,n) = A^cTr(r"i . . . T'^^-i) Tr(r''^ ...T""), j > 1. (9.12) 

The color-ordered amplitudes contain terms which depend on the number of quark fla- 
vors, for example, with factors Nf/Nc and {Nj/Nc)'^. In addition, at two loops An-i-i also 
contains terms of order arising from both planar and non-planar diagrams. In con- 

trast, while the two-loop gluon splitting amplitude does contain contributions of 0{Nf /Nc), 
0{Nf'^/N^), and 0{Nf/N^), as we shall see below, it contains no contributions of 0{l/N^). 

The collinear behavior of the leading-color amplitudes ^n;i;i was given in eq. (|2.19| ), 
where the notation ^A^ was used instead of ''An-i-i. At two loops the subleading-color 
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amplitudes cannot be determined solely from the leading-color ones. However, just as at 
one loop, consideration of the trace structure leads to the following collinear limits for the 
subleading-color amplitudes. The double-trace partial amplitudes behave as, 

..y^y^...)^ (spiitL°i(.; , b^^) ..,p\...) + 



+ SplitL'l(z; a^" , b^'^) Ai:i,,,_,{. . . , P^ . . .) ) , a,b < c, 



..y^y^...)^ J2( spiitL°i(.; a'^y) A^:i,y. ..,p\...) + 



(0). .^A„ ,Ai,x .(2) 



+spiyiiz;a^^y)A[yj 



while the triple-trace partial amplitudes behave as, 

4(2) ( „A„ ,Ai, N 

SplitL°i(z; a'^ , b'^) 4'li;c,-i;c.-i(- --^P' 



a,b > c. 
(9.13) 



A=± 



Am 

^n;ci;C2 V' • • J 



(...,a^",5^\...) 



allfe 



Y SplitL°i(z; , 6^^) 4'2l;c.;C2-l(- • • > ' 
A=± 

/l(2) ( „Aa tAfc \ 



a||b 

5]SplitL°i(z;a^",6^^)4^1i^,^^,^(...,P\.. 



a,b < ci 



ci < a,b < C2 



a,b > C2 ■ 



(9.14) 



A=± 



As at one loop, legs a and b must be cyclicly adjacent within the same color trace, otherwise 
the collinear limit is finite. 

At two loops, the generalization of the splitting amplitude from the leading-color struc- 
ture to the fully color-dressed version would seem to be more complicated, because of the 
appearance of non-planar contributions, such as the ones depicted in figure 28. However, 
as already mentioned in section ^, all such non-planar contributions vanish for the g ^ gg 
case. The fully color-dressed splitting behavior is given simply in terms of the color-stripped 
splitting amplitudes of section as, 



j^{2) aia2a3...an ^W^^ ^ Jaia2ap 



Y Ni spiit^\(z; 1^^ , 2yA^y^ 



(9.15) 



The proportionality of the color dressing to J^i^aap follows from the diagrammatic 
argument given after eq. (|9.8| ), generalized to two loops. To see why the pure-glue contri- 
butions to the splitting amplitude just have an overall factor of A^^ in front (as evidenced 
by eqs. ( 7.38 ) and ( 7.43| )), consider again the parent diagrams. We have already seen in 
section H that the non-planar diagrams do not contribute at two loops to the splitting 
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Figure 32: A diagram contributing to the three-loop g ^ gg splitting amplitude which has 
subleading-color, suppressed terms in its color factor. 

amplitude to gluons, so we focus on the planar diagrams. At two loops there are four 
kinds of planar parent diagrams. Each of these has either a bubble or triangle subdiagram 
(the two parents with triangle subdiagrams are shown in figure |2^ . Upon replacing the 
subdiagram factor of —Ca or —^Ca times a vertex, respectively, we are left with a triangle 
diagram, which generates a second factor of Ca, as at one loop. All pure-glue contribu- 
tions are therefore homogeneous in Nc of degree two. This latter argument does not extend 
to higher loops. Beginning at three loops there are pure-glue diagrams, both planar and 
non-planar, which have no surviving triangle or bubble subdiagrams at some stage of this 



reduction. The simplest example is depicted in figure 32. Such diagrams generally give 
rise to suppressed terms. They correspond to the existence of other SU(iVc) group 

invariants, such as {d"'^'^)'^, where d"'^'^ is the fully symmetrized trace of three generator ma- 
trices, which are not homogeneous in Nc- Of course, the fermion-loop terms in the g ^ gg 
splitting amplitudes have l/N^ suppressed contributions already at two loops, due to the 
existence of both Ca and Cp Casimir invariants. 

10. Conclusions and Outlook 

In this paper, we computed the behavior of a general two-loop amplitude in massless 
gauge theory as the momenta of two external gluons become collinear. (The behavior 
involving collinear fermions will be presented elsewhere |^^.) This behavior is universal, 
and accordingly governed by a set of 1 ^ 2 splitting amplitudes. For the g ^ gg splitting 
there are three independent helicity configurations (the remainder are related by symmetry 
or parity). One of these configurations vanishes at tree level. In an = 4 supersymmetric 
theory, the ratio of a non- vanishing two-loop splitting amplitude to the corresponding tree- 
level amplitude is helicity-independent, and given by eq. ( 7.16| ). In pure = 1 super- Yang- 



Mills theory, the helicity configuration that vanishes at tree level also vanishes to all orders 
in perturbation theory. The ratios for the other two cases are given byeqs. ( ^^Sj) and (|7^). 



In QCD, all three helicity amplitudes are non-vanishing; the ratios to the tree for two of 



them are given by eq. (7.38|), and the splitting amplitude for the remaining configuration 



by eq. ( 7.43| ). The finite remainders of amplitudes, after subtraction of poles predicted 



by Catani's formula, also have simple behavior in collinear limits in the 'color-trivial' 
case. These limits are given at one and two loops by eqs. (I8.21D and (|^), respectively. 
The functions entering these limits are given by eqs. ( S.31| )- (|8.34| ), ( ^.35|) , ( 8.37 ), ( 8.3^ ), 
and (pOl)^(p^. 
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These splitting amplitudes can be used as a check on calculations of higher-point two- 
loop QCD amplitudes; the collinear limits of these amplitudes must satisfy eqs. ( 7.3^ ) 
and (7.43). One can also apply the checks to the finite terms after subtraction of the poles 
predicted by Catani's formula, using eq. ( 8.2S| ). As mentioned previously, the compatibility 



of the divergent parts of our splitting amplitudes with Catani's formula provides a inductive 
proof of the latter, modulo some assumptions about the analytic behavior of the e-singular 
terms in amplitudes. 

The two-loop splitting amplitudes are also one of the ingredients required for computing 
the NNLO corrections to the Altarelli-Parisi kernel in an infrared approach. At NLO, one 
would add two terms |43|: the one- loop 1 — > 2 splitting amplitude interfered with its tree 
counterpart; and the tree-level 1 — > 3 splitting amplitude squared, integrated over the phase 
space of the unobserved partons. At NNLO, there are three ingredients: (a) the 1^2 
splitting amplitude computed here, interfered with its tree counterpart; (b) the interference 



of the one-loop |£5|| and tree-level ^ 1^3 splitting amplitudes, integrated over the 
unobserved phase space; and (c) the 1 —>■ A splitting amplitude squared [Q, integrated 
over the four-particle unobserved phase space. 

In AA = 4 supersymmetric gauge theories, the splitting amplitude has a remarkable 
property: it can be written as a polynomial in the one-loop and tree-level splitting ampli- 
tudes. This led to the conjecture that a similar property holds for two-loop amplitudes p3|| . 
Explicit calculations showed that this conjecture is true for the four-point amplitude. This 
simplicity suggests that substantial parts of the theory may be solvable. We confirm a 
leading-transcendentality relation between the = 4 supersymmetric and QCD splitting 
amplitudes at two loops; this posited relation was already used |5^] to extract the AA = 4 
leading-twist anomalous dimensions from those in QCD |48| |. 

We performed this computation using the unitarity-based sewing method. The method 
has demonstrated many advantages over conventional diagrammatic techniques over the 
years. It has made possible, for example, the computation of series of amplitudes for arbi- 



trarily many partons 11]. The present computation demonstrates another advantage: 



it offers a pathway featuring the physical-projection advantages of light-cone gauge, while 
avoiding the need for cumbersome prescriptions for dealing with the ill-defined integrals of 
the latter. Because it effectively combines many diagrams at an early stage of the calcula- 
tion, it also simplifies integrands considerably, and so minimizes the complexity of reducing 
loop integrals to master integrals. 
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A. Color-Space Collinear Limit of Catani's Formula 

In section ^ we demonstrated the consistency of the divergent part of the sphtting amph- 
tudes with the Catani formula at leading order in Nc, or more generally, for the color-trivial 
single-trace terms in the formula. In this appendix we demonstrate that the consistency 
also holds for the terms containing non-trivial color correlations. Because of these correla- 
tions, the collinear limits of the full Catani formula are a bit intricate. In this appendix, 
all quantities are understood to be renormalized. 

A.l Color-Space Notation 

In the color-space language of refs. |^, Q an amplitude is expressed in terms of an abstract 
vector in color space, {aH^^). To convert to the more standard color notation one may use 
an orthogonal basis of unit vectors |ai, 02, . . . , a„) with the property that 

^(L)aia2...a„(^^^ ^^2, A2; . . . ; K, A„) = (ai, ^3, . . . , a„|^(,^)) . (A.l) 

Color interactions are represented by associating a color charge Tj with the emission of 
a gluon from each parton i. The color charge Tj = {T^-} is a vector with respect to the 
generator label a, and an SU(A'^c) color matrix in the representation R of the outgoing 
parton i. For external gluons, R is the adjoint representation A, and T^^,^ = if^""^. (Note 
that the normalization of fundamental representation charge matrices is different in ref. [^8| 
from that used elsewhere in the paper. In this appendix we normalize the Tj according to 
ref. m.) 

In this notation each vector \An) is a color singlet, so color conservation is simply 

n 

Y,T^\An)=^, (A.2) 

i=l 

independent of the color representation of each leg. This identity incorporates the Jacobi 
identity and its generalizations very simply. Typical operators encountered in the discussion 
involve the combination 

{Tif{Tjf = Ti-Tj. (A.3) 

For i = j, eq. ( [A.3| ) reduces to a Casimir operator, Tf = Ci = Ca = Nc if leg i is a gluon, 
and Tj = d = Cp = (N^ - l)/(2iVc) if leg i is a quark or anti-quark (with Tr = 1/2). A 
useful property is 

Ti ■ Tj = Tj ■ T, , (A.4) 

which holds because charge matrices act on different index spaces. 

In the collinear limit where ki — > zkp, k2 ^ (1 — z)kp, a tree-level color-space ampli- 
tude satisfies 

14°)) - E SplitL°i \Ai'Ux)), (A.5) 

A=± 

where A denotes the helicity of P and the splitting functions Split^| are now operators 
acting on the color space. The operator Split links the color space with n — 1 legs to 
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Figure 33: The action of fp on tlie color index of leg P of an amplitude. The vertex represents 
a factor oi ip^"^"'' . 

that of n legs. That is, we define 

(ai,a2,a3, . . . ,a„| Split |^^^2i(A)) 

= SplitL°l 4-r"''""(^P> A; fcs, A3; ... ; fc„, A„), (A.6) 

where Spht^]^ is the color-ordered splitting function defined in eq. (|2.3| ). We define an 
operator 

where for the pure-glue case we have 



SplitL°l = /pSplitL°l, (A.7) 



(ai,a2,a3,...,an|/p|ap,a3,...,a„) = . (A.8) 

The operator fp is distinct from the operator Tp because it links color spaces of different 
dimensions. The action of the operator fp on an amplitude is illustrated in figure 
Color conservation on the three-vertex implies that 

fpTp = {Ti + T2)fp, (A.9) 



which is a special case of eq. (|A.2|) (up to a sign due to a swap of incoming and outgoing 
indices). This equation is equivalent to the Jacobi identity. 

From eq. (|9.8|), it is apparent that the behavior of the one-loop amplitude in color 



space in the collinear limit is very similar to that of the tree amplitude, 

A=± A=± 

where 



Split 2^2) = NcSplit^}l{l^\2^^) fp , (A.ll) 



with Split^l the color-stripped splitting amplitude. For helicity contributions where the 
tree splitting amplitude does not vanish, it is convenient to re-express this as 

Split^}lil^\2^^) = Ncf'g^'^'^^ {z, s) SplitL°l(l^\2^2)^ (A.12) 



where is given as in eqs. ( 2.15| )~( p.l7| ), and is renormalized according to eq. (^). 



The L-loop generalization is the obvious one, 

L 

l-^i^^) -EE Split^l \^t1\x)) , (A.13) 

1=0 A=± 
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where 

Split (5 = Split fp . (A.14) 

In the two- loop case, Split is the color-trivial splitting amplitude defined in eq. ( |2.19| ). 

Next we shall demonstrate that the full color-space dependence of Catani's formula 
for the divergences of one- and two-loop amplitudes is completely compatible with the 
divergent parts of our color-dressed splitting amplitudes. 

A. 2 One-Loop Warm Up 

At one loop, the infrared divergences of renormalized n-point amplitudes can be written 
compactly in color-space notation as |38], 



The color operator In'^ is 



1 e-^V'Ci) 
2r(l-e) 



1 



l^ 1 



£2 + e 



(A.15) 



(A.16) 



i=l j^i 

The sum runs over pairs of external legs. For external gluons i, we set 7^ equal to 

llCA-4.TRNf 



79 



6 



(A.17) 



(For external fermions, the ratio Jq/Tf = 3/2 is independent of the representation.) The 
subscript as. indicates that a quantity depends on the regularization and renormalization 
scheme. 

Now consider the limit of the In^ operator as the momenta of legs 1 and 2 become 
collinear. Inserting the decomposition of the amplitude into divergent and finite parts, 
eq. ( |A.15 ), into both sides of eq. ( A.1C| ) and taking the collinear limit, we find that the 
divergent parts must satisfy 



/W(6)SplitLll^ri(A)), 



I /l(O) 



Split«|4°l,(A)),s + 



+ Split(°1^2i(e)l4-i(A))R. + finite, (A.18) 



where collinear kinematics are implicit. On the left-hand side of the equation In\e) acts 
on an n-dimensional color space, while on the right-hand side I^}_^{e) acts on an (n — 1)- 
dimensional color space. A compact way to write the same equation is 



spiitL'i 



/«(e)SpiitL°i 



SplitL°llJ2i(e) + finite. 



(A.19) 



To make this more explicit, consider an n-gluon amplitude, for simplicity. Inserting the 
explicit forms of the operators, we have that the infrared-divergent parts of the renormalized 
one-loop splitting amplitude must be 



SplitL'l 



r(i 



1 . 
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J=3 

2 



J =3 



+ Ti • T2 



-S12 



-(1 - Z)spj 



Split(°l- 



+ 



i^j=3 



spiitL^l 



.J=3 



-Spj 



+ finite, 



(A.20) 



where we have separated the terms into those involving the colhnear legs and those not 
involving them. The terms that do not involve the collinear legs simply cancel, since 
Split ^1 commutes with any Tj where j > 3. After separating terms that involve spj from 
those that do not (allowing for shifts in the finite terms), we find. 



SplitL'l 



r(i - e) 



e2 ^ e 



eilnz^Ti-Tj + ln(l - ^ T2 • + 
V j=3 j=3 J 

2 



+ Ti-T2 



-Sl2 



spiitL°i + 



+ E ((Ti + T2) • T, SplitL^l - SplitL^l Tp . T,) (^) ^ } 



finite . 



(A.21) 



In the above, we have been careful in the placement of Split since it links (n — l)-point 
color space to n-point color space. To simplify eq. ( |A.21 ) we use eq. ( |A.9| ). Furthermore, 
using color conservation, eq. ( [A. 21 ), we have 



(0) 



Ti ■ Tj Split = -Ti • (Ti + T2) Split_^ 

J=3 



1 



.(0) 



2^,Ti + T2)' + Tf-T^)Split^;, 
ic-A Split 



(A.22) 



with a similar equation with leg 1 replaced by 2. Also, Ti-T2 \{{Ti+T2f -TI-TI) 
— i^Ca- Inserting these relations into the divergent parts of the renormalized splitting 
amplitude gives us 



SplitL'i 



-Ca' 



1 e-^'/'(i) 



2r(l-e) 



CAe 



A* 



-Z{1 - Z)si2 



spiitL°i + 



+ finite . 



(A.23) 



This result is in complete agreement with the divergent parts of the one-loop color-space 
splitting amplitude, as given in eqs. (2.14) and ( A.12| ), after renormalization according to 
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eq. (8/7). (Renormalization produces the term proportional to 7^ = 60 iii eq. ( |A.23D .) A 
similar analysis using mixed quark-gluon amplitudes yields the same result for the g ^ gg 
splitting amplitudes. This demonstrates that the expression ( A.15| ) for the divergences of 
one-loop amplitudes is fully compatible with the splitting amplitudes. 



A. 3 Review of Catani's Two-Loop Divergence Formula 



In a beautiful yet mysterious paper, Catani expressed the infrared-divergent parts of two- 
loop amplitudes as 



l4'^)Ra = I^nUie) +I«(e) |4'^)aa + |-4('^'^'^>k 

where |^i^^'^'^)Ra is the finite remainder and the operator /i^Ra(e) is 



a ! 



(A.24) 



+ r(l-e) 

An argument for this general structure has been given in ref. |^2[. The quantity -ftTna 
depends on the variant of dimensional regularization through the parameter 5^ |15]. It is 
given by [||, ^] 



K 



K 



FDH 

■M=l 
FDH 



-C2 Ca , 



T^QCD 



3-C2 

67 



18 



C2 



Ca , 



Ca - ^TnNf . 



(A.26) 
(A.27) 

(A.28) 



Note that in passing to the = 4 case we have switched away from the conventions of 
ref. [15 1, with regard to assigning e-dependent terms to i^PDH- (We did not want to destroy 

the uniform transcendental weight which all functions have in the AA = 4 case.) 

(2) 

The function aa contains only single poles. It splits into two types of terms, 



4er(l -e) 



i=l j^i 



3 rj-i2 



+ H 



(2) 



(A.29) 



(2) (2) (2) 

where H^^ is either Hg or Hg , depending on whether particle i is a gluon or a quark. 



The constants and H^"' are given by ||, |2|, |15|, |3|] 



r(2) 



+ (^-IcaCf + ^Cl + ^CpTnNf^ (1 - 5r) , 



(A.30) 
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and (including also the supersymmetric cases here [15, |3 



<)'^=^ = |C1, (A.31) 
i^r'^=^ = (| + ^C2-0Cl, (A.32) 
= (I + + - (I + '^y^TuNf + CTnNf + |r|iV/ + 

+ (-^^1 + IcaTrNj^ (1 - 5^). (A.33) 



There is no real proof that eq. ( A.29| ) is the precise form of the n-point divergences, but 
previous calculations strongly indicate that it is right. There are several conventions in 
the literature for 'dressing' the terms proportional to H^"^^ with factors of {fJ-'^ /{—Sij)Y to 
make them dimensionally consistent. Different conventions agree at 0{l/e), but generate 
different finite remainders, which would have different collinear behavior, affecting the 
results of section H. Here we adopt a convention motivated by Catani's treatment |pq] of 



the 1/e poles proportional to ji at one loop in eq. (A. 16). Note that the 1/e pole part of 
the terms containing -ffj-^^ in eq. (A. 29) is actually proportional to the color identity matrix 
1, since - Ej^i^i ■ Tj = Tj, due to eq. Q). 



Less is known about the function (e), which contains genuinely non-trivial color 
structure at 0{l/e). For the four-point amplitudes it is given by |15, p^] 



(e) = -4 In 



-Sl2 
-S23 



In 



-523 
-Sl3 



In 



-■S13 
-Sl2 



T^-To,To- T, 



(A.34) 



with ln((— si2)/(— S23)) lnsi2 — ln(— S23) — «vr in the s-channel, etc. A simple ansatz, 
generalizing eq. ( |A.34| ) to arbitrary n-parton amplitudes, is 



(il,*2,«3) 



a2a3rpairpa2rpa3 j 



In 



'«2«3 



'n«3 



In 



'n«3 



(A.35) 



where the sum is over distinct triplets of external legs, with 7^ «2 7^ ^3- For n = 4, there 
are four such triplets (omit any one of the four partons). Using [T^,TI'] = if"'^'^T^, color 



conservation ( A. 2 ), and antisymmetry of /"^^, it is easy to see that each triplet gives an 
equal contribution, and eq. (|A.34|) is recovered. 

A. 4 Collinear Compatibility with Catani's Two-loop Divergence Formula 

Next let us verify that the e poles in the collinear splitting amplitudes match those predicted 
by Catani for the full color-space two-loop amplitude. In order to analyze these divergent 
terms, it is convenient to define the one-loop amplitude as being a color operator acting 
on the tree amplitude (assume it does not vanish), 



= (/«(e)+M«fi-(6))|^W)H.s, 



(A.36) 
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since we want to combine divergent and finite parts into a single entity. We do not actually 
need the explicit form of the Mn^^^ operator, although it is straightforward to construct 
this operator in one's favorite color basis. 

The full color-space two-loop amplitude is 



2 " ^ ' V " ^ ' e J r(l-e) 
|^W)r.s+ finite. 



/«(2e) + 
(A.37) 



After adding and subtracting finite terms to combine the term with finite pieces, in 

much the same way as was done for the color-trivial parts in eq. (|8.16|), we obtain. 



n /R-S 



1 

2 

+ 

+ 



r(l) 



1 



r(2), 



4e 

e+^V'(i)r(l - 2e) 



rfi 



4er(l 



(2) 



e 



r(i), 



2el 



2e 



finite. (A.38) 



First consider the terms on the penultimate line in eq. ( |A.38| ) containing the operator 
In^ . Since the same operator appears as in the one-loop case, the collinear limit of these 
terms may be determined in the same way as the one-loop case, using eqs. ( A.19| ) and ( A.12| ) 
in particular. Equation ( [A.19[) says nothing about the 0{€'^) terms in /^^^ . Fortunately 
these terms are identical for In\&) and /n^(2e); hence their contribution to the singular 
terms in eq. ( A.38| ), when they are multiplied by bo/e, cancels. Following the one-loop 
discussion, we obtain contributions to the two-loop splitting amplitude of the form 



-ry{e) + 



e+^^(i)r(l - 2e ) fbo 

r(i - e) 



(2e) ) Split + finite . (A.39) 



(Recall that r^^\ introduced in section ^, denotes the renormalized splitting ratio.) 

Now consider the terms containing H^-^^ on the last line of eq. ( A.38| ). These terms are 
written in the same form as ill \ but with one less power of l/c, so their collinear behavior 
again follows from the one- loop discussion. Their contribution to the two- loop splitting 
amplitude is just 



(2) 



4e 



Split + finite , 



(A.40) 



where Hg'^^ is given in eq. ( A.33| ) 
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Next consider the {Mn\(-))'^ term in eq. ( |A.3^ ). In the collinear limits for this term 
we have, 



l(M«(e))'splitL°il^il) 



(A.41) 



In order to evaluate this we use the collinear relation, 



SplitL°iM«, + SplitL^ 



(A.42) 



which is a rewriting of eq. ( |A.10 ) after acting on |w4^^2]^) and using eq. (lAlS^ ). Using 
eq. ( |A.42 ), we find that in the collinear limit, 



l(M«)'splitL°l ^ iM«(SplitL°lMi^l, + Split(^) 



To simplify the Split ^| term we may multiply eq. ( [A.42 ) by Ncr^g\e) to obtain 



Split ^ Split M^^l, + Split N,rf{e) 



(1) 71^(1) 



(A.43) 



(A.44) 



where we used eq. (|A.12| ). Then applying eqs. (|A.42| ) and (|A.44D to eq. ( |A.43D yields 

l(M«)'split(°i ^ i(SplitL°iM«,M«, + 2Split« + Split« iV,r«(e) 

= SplitL°i I M^:lMl. + Split« + liN4\e)f SplitL^l . 

(A.45) 

The first term on the last line may be identified as a contribution to Split the 
second as a contribution to Split^| while the third is the contribution we include 

mSplitLll4°2i). 

Finally we examine the second line of eq. ([A.38| ), containing the commutator term and 



. After some manipulations involving eqs. ( A. 42]) and (A.44) we may write the part 
of the commutator term that does not contribute to Split^| j^A^^^) as 



-ImW (i« SplitL^l - SplitL°l li^l, - SplitW) + 



+ l(lWSplitL°i 



Split(°i 1^2, - SplitLl ) {MI:L, + N^f','' ) 



r(l) 



(A.46) 



where the finite combination appearing is, 



/«spiitL"l-spiitL"l/i^2i- Split 

n 

\nz Ti +ln(l - z) T2] • ^ In(-spj) Split + 0{e) . (A.47) 



J=3 



Only the universal 1/e terms in Mn contribute in eq. ( |A.46 ) to the order we need, and 
we may obtain these terms from ^. Inserting these terms, and moving them to the right 
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with the help of color commutators, the total contribution to Split '^^^ |„4^^2]^) 
commutator term is 



from the 



I 

Ye 



In 



n 

f-'-T^T^rr H-spi) (in(-spi) + Inz + ln(l - z) - ln{-su)) + 



j=3 



i=3 j=3,jy^i 



Split . (A.48) 



This result turns out to be precisely the negative of the corresponding contribution from 

^ (2) 

the iJ„ term. Thus the second line of eq. ( A.38 ) only contributes to the tree- level splitting 
amplitude term, Split^j \A^li{X)), in eq. ( A. 13 ). 

Combining the contributions to Split^^l in eqs. ( [A39D and ([Ool ) and ( [OsD , we 
obtain. 



Split 2^2 ^ 



SplitL'l(l^\2^2)/p 



bo Ui).. , e+--^Wr(l-2e) 



.(1) 



(26) + 



+ 



1 H. 



(2) 



+ finite 



SplitL°l(l^\2^2)/, 



4e 

N^r'§\e) SplitL°l(l^S2^2)y.^ + finite. 



(A.49) 



(2) A K 

where rg [e) is given in eq. (|8.18D . This in turn agrees with the divergent parts of the 
splitting amplitudes for AA = 4 and AA = 1 supersymmetric theories, and for QCD, as given 
in section after they are renormalized according to eq. ( |8.8D . Thus our splitting ampli- 



tudes are fully compatible with Catani's color-space formula for the infrared divergences 
of scattering amplitudes. 

As noted already, this agreement may be turned around to prove the validity of Catani's 
formula inductively, up to reasonable assumptions about the analytic structure of the e- 
singular parts, i.e. such that they do not have vanishing collinear limits in all channels. 



B. Relabeling Algorithm 

For completeness, in this appendix we review a simple procedure for mechanically per- 
forming relabelings. For the splitting amplitude problem discussed in this paper, it is not 
difficult to identify the diagrammatic structure of a given term and map the momentum 
labels used in the cuts to those used in the integration. However, for more general problems 
it may be useful to perform these steps algorithmically. To do so, first one chooses a set 
of parent diagrams and momentum labels used in the integration. To sort the terms into 
diagrams it is useful to collect the terms on each propagator or light-cone denominator 
type. The propagators and light-cone denominators encode the diagram type, but with 
labels not matching the ones used in integration. Assuming the external legs are ordered 
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the same way in the cut and integration labels (if not, one should permute the external 
legs), the problem amounts to finding a mapping between a given term in the cuts and the 
diagram to which it belongs. For an n-point L loop diagram there are L independent loop 
momenta and n — 1 independent external momenta, under momentum conservation. If the 
loop propagators carry momenta pi in one set of labelings and qi in another labeling, then 
a change of variables should exist relating the two of the form, 



n-l 



= ^ aijPj + ^ hijkj , (B.l) 
i=i i=i 

where the aij and bij are in the set {—1, 0, 1}. This form assumes that for both the original 
and final momentum labels the momentum of each internal leg is a sum or difference of the 
independent momenta, as naturally arises either in cut or Feynman diagram momentum 
routings. One sweeps over all changes of variables in eq. ( [B.l| ) until a match is found. If 
no change of variables matches then the subsequent diagrams should be checked until a 
match is found. An efficient way to rule out putative changes of variables is to first set 
all external momenta to zero, then re-introduce them one by one, eliminating candidate 
diagrams and solving for the hij. 

As a simple example, consider the double triangle diagram with a light-cone denomina- 
tor occurring in our calculation of the two-loop splitting amplitude. From its origin in the 
three-particle cuts shown in figure ^6|(b), where the cut momenta are {gi, —qi — q2 — ki — 
k2,Q2}, this diagram is described by the set of propagators and light-cone denominators 

{ql,ql {q2 + fe2)^ + h + kif, +qi + ki + /C2)^ {qi + ki + kif, qi-n}. (B.2) 

The denominator variables used in the integration routines, shown in figure ^9|(a), are 

{Pi, {Pi +ki+ k2f,pl, {p2 -ki- k2f, {p2 - kif, {pi +P2f,Pi • n} (B.3) 

where p2 is labeled as p^ in the figure. After sweeping through the change of variables in 



eq. (B.l) one finds that the two sets of momenta of momenta are related by 

qi=Pi, q2=P2- ki- k2. (B.4) 



With this change of variable the two sets of propagators and denominators in eqs. (|B^ 
and (|B.3D are identical, except for the order in which they appear. 
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